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HAUPTAUFSÄTZE 


Die Strömung um einen ebenen Tragflügel mit Querruder. 


Von €. Schmieden in Rostock !\. 


nläßlich theoretischer Untersuchungen des Flugtechnischen Instituts der Technischen 

Hochschule Berlin ergab sich das Problem der konformen Abbildung des Außenbereichs 
eines gewissen Flügelprofils auf eine Halbebene. Dieser Aufgabe hat sich in dankenswerter 
Weise Herr Prof. ©. Schmieden, Universität Rostock, unterzogen; dabei ergab sich außerdenı 
das Strömungspotential für den unendlichen Parallelstrom um den ruhenden Flügel und die 
kraftwirkung der Strömung auf den Flügel. Wir bringen hier die Entwicklungen von 
. Schmieden im wesentlichen wörtlich wieder (nach einer brieflichen Mitteilung). 


Wir idealisieren das Flügelprofl des ebenen Tragflügels mit Querruder durch zwei 
segeneinander abgeknickte Strecken BE und ED (Abb. 1). Wir nennen die Tiefe des Quer- 


ruders ED=tp, die des Flügels ohne Querruder BE tr und die Gesamttiefe des Trag- 


’ R : r A BE . 
Nügels I, tp+tp; den Ausdruck } } setzen wir gleich 4, wobei O<i<1; y sei der 
Ft Ir 

tuderausschlag (zunächst y >60). Das Profil des Tragflügels ist durch den Streckenzug BUCDE 


lest veleet. 


Il. Abbildung der z-Ebene auf die obere !-Ebene (Abbildung des Streckenzuges BUÜDE 

er Abb. 2 auf die reelle Achse der tEbene von Abb. 4). Es ıst bekannt, daß man die 

bene Strömung eines unendlichen Parallelstromes um ein beliebig gestaltetes Flügelprofil 

geben kann, wenn die konforme Abbildung des von dem Profil bestimmten Außenbereichs 

r 2-Ebene auf das Äußere eines Kreises bzw. auf eine Halbebene gegeben ist. Hier wird 

gewisser Weise die Sache umgekehrt, indem von einer speziellen Strömung um das Profil 
(U die konforme Abbildung geschlossen wird. 


Der Tragflügel werde mit der Anströmgeschwindigkeit V,, parallel zur “-Achse, von 
ks her angeströmt. Dabei möge der Flügel den Anstellwinkel a besitzen (Abb. 1). Die 
'kulation um den Flügel wird in der üblichen Weise durch die Forderung „zlatter Abfluß 
der Hinterkante* festgelegt. | 


!) Für den Druck fertiggestellt von Dr. G. Ellen berger, Flugtechnisches Institut der Teehnisehen Hochschule 
ın (Prof. Dr. Ing. H. Wagner). 
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Im allgemeinen wird der vordere Staupunkt A irgendwo auf der Unterseite bzw. a 
der Oberseite des Tragflügels liegen, so daß die Vorderkante mit unendlich großer G 
schwindigkeit umströmt wird. Man sieht aber sofort, daß bei jedem fest gegebenen A und 
ein gewisses «a, existieren muß, für welches der vordere Staupunkt in die Vorderkante 
hineinfällt (Abb. 2). Da in dieser Lage die Geschwindigkeit auf der Unterseite durehschnii 
lich größer sein muß als auf der Oberseite, muß sich für diesen Fall ein Abtrieb ergebe: 
Von dieser Lage, wobei «a, natürlich eine Funktion von / und y ist, gehen wir aus: Betracht: 
wir die Geschwindigkeit entlang des Tragflügels, so ist sie überall endlich, nur an d 
konvexen Eeke E wird sie unendlich (an der konkaven Eeke Ü wird sie Null). Durch di 
Funktion 
lz 1% 


( 

e In—” +i-0 

, In | f d = er 

(V Betrag, © Richtung der wahren Geschwindigkeit) 


wird somit in diesem Falle die z-Ebene abgebildet auf einen Parallelstreifen mit dem Null. 





punkt im Innern. Der obere Rand hat die Gleichung 9=y-— a,, der untere ©: Gy. (Es 
ist klar, daß für den betrachteten Strömungszustand y > «a, sein muß Abb. 3 —.) 
z-£bene Y z- bene Yy 
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Abb. 1 bis 4. Der Zusammenhang zwischen der z-, <- und £-Ebene. 
Bilden wir den Streifen so auf die obere #f-Halbebene ab, wie in Abb. 4 angegeben, so 
erhalten wır 


® ,„ da Y 
Ä nV.) end. a N Fe Vase De N, 


Dem Unendlichfernen der z-Ebene entspricht CO, also 


.ao 
ı 7 


f e ER a a Ne BT De a DE ee et ee 2). 


Iım Punkte #, haben wir jetzt einen Dipol und einen Wirbel anzubringen und das komplexe 
Potential so zu bestimmen, daß die reelle f- Achse Stromlinie wird und die Staupunkte naclı 


t=z°" und = — v* fallen. Die erste Bedingung erfüllt der folgende Ansatz: 
2 ei? ei? pr ” N eiun (b reelle freie (3) 
m db - z ! a 3). 
Ferien eier 22 t-—e-i«r| Konstante) | 


Kine weitere Bedingung, der wir die Konstanten zu unterwerfen haben, ist die Schließungs- 
dz 


- 


bedingung: Integrieren wir N ın der f-Ebene entlang der reellen Achse von — oo bis +x. 
{ Cu 


so müssen wir in der z2-Ebene an den Ausgangspunkt zurückkommen, also 


(az, (dedw I | dm N 
\ur FT v.\ea‘ 


ar 
Da sich für >» der Integrand verhält wie f” , können wir das Integral über den Hal 
kreis mit unendlich großem Radius in der oberen Halbebene hinzufügen, ohne den Integralw‘ 
zu verändern. Es muß also das folgende Konturintegral verschwinden: 


P N 1 » ! iß8 , i?3 IR / { N 
- dt: ) € — ( -—. ’ } At. 
\ Id ein? (eier TIalt-— ein _e-iun]f 


I « 


Kontur: ID 
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‘ dem Rande liegt keine Singularität, die einzige Singularität des Integranden im Innern 
" Kontur liegt bei £=t,. Also muß das Residuum bei #, verschwinden. Wegen 


f 7 (eier Lt eiur) a - ei li + - E iur(g eier) >... ) 


ot 
we y / 
Rele!)=ef-#M) — +, 0, 
nı 2m 
ler | 
: 0 untn (4) 
Pia WEN ee at 2 ae Da ee 


s ist 2'>0 für y>0, d.h. wir erhalten tatsächlich einen Abtrieb, wie oben behauptet 
vurde. Wir erhalten jetzt 


ih 


d mw | e° jun P iun y | ' 
d l dt e iu 7T )? (t u ei un j? pr r . 


ri / +2 cos un f ı ") (»). 


„T \ JT 


(’—- 2tcos un + 1) 


2b.sin un 





Die Staupunkte fallen in die Wurzeln von 


| ‘ 
yv cos ur 7 T 
+2 - ——t- 0, 
IT y y 
| | 
„T „T 


also wird 
ar 2 
„? E.. | —| sin? ı \ | 
' 1 | P% COS Ani 


N (6). 
v” = ” | | en sın? u a+Z CoS u " 





x" und »* liegen beide dicht bei Eins (für kleine y), sind aber größer Eins für y>V®. 


ür z erhalten wir nunmehr das Integral 


% V it - e iu 1)? (t eine Tr 


il e iag( ’ e-iun ein 7 | | 1 \ 
2) - t st J 4 ‘ \; R j ) 
je e) 


Nun ist 


| t7 eiunm t7 eiun y | ein r f7 
I 
(t 


„dt r — At, 
BrrFy t—ettn mn )tlt—eiun) ° 
f7 e- iu f7 pP tun u P iunm A 
| wert — . dt. 
. (f ed / p ! E tw? 7ı k tit E ur) 
\\eiter 
id pn pn n 
y eiurf: Ä | En 
e ..: 4 “ 
er ame Zirtar, 


Umemer: Fr | & | AR 
| . dt. 


- Be Be Ta] 41 = 7 


«/ 





! kn Ze N rin e \ Ztschr. f. ang 
196 Schmieden, Die Strömungz um einen ebenen Tragtlügel mit Querruder Math. und Mr 





Also heben sıch die übriebleibenden Integrale fort und es folet 


ibe-iw | eiun e-iun 

2 - / 7 u “ . im = Pr . . . . . . . (7 

V, (per —e-iun] 
oder zusammeneereehnet 
1+ 

I, sın ur-f 

2 — gta — > m) er Armee en (S 
® f° 2tcos un + | 


wobei der Koordinatenanfang in dem Anschlußpunkt des Ruders an den Flügel liegt (U bzw. / 


Für die Gesamttiefe des Flügels erhalten wir 


a 4 | zig 3 Me \ 
2b sın un (#°) ad Miu 
ER ER ‘ 
7 r =: Bu, | 3 20 (.) 
V. „"—Zxtcosun+1l »*—2v’cosun- | 
während man für 4 erhält 
„1i+2 
(vr) 
‚ vv! —-2v’cosun-tl ( 
= — a er a a a a  ° 
F \ I 
2b sın un 
Werden Quadrate von — vernachlässigt, so erhält man aus (6), (9) und (10) 
7 
u ' | COS UT 
h = y.smun 4 = aa ea Ka ee SACHE, 


Wir machen uns jetzt von der Strömung, die zur Auffindung der Abbildung führte, 
frei, indem wir setzen 


7 - . j . 
T, z | e sn ur. 


wo f nach (9), (10) und (11) eine positive Funktion von y und 4 ist, die gleich Eins ist, 


‚ei statt 2 (d.h 


[0 


wenn (Juadrate von —- vernachlässigt werden. Ferner schreiben wir 
T 


wir drehen das Profil um den Winkel a,) und erhalten so 


a 


5 (Sa 
2t cos MIT | 


z f,-sın’un-f j: 


aıst mach (10) bestimmt als Funktion von y und A. Die Abb. (Sa) enthält keinerlei Größen 


mehr, die sich auf eine Strömung beziehen. Sıe gilt für „,<1, also auch für negatives 


>> 
(die „-Achse fällt jetzt in die Riehtung von BC). Die gefundene konforme Abbildung ist : 
Beispiel für den Schwarz-Uhristoffelschen Abbildungssatz. 


Il. Auftrieb und Moment bei beliebigen Anstellwinkeln.. Nachdem wir die Abbildung 
funktion (S) kennen, können wir Auftrieb und Moment bei beliebigen Anstellwinkeln berechne 
Dazu brauehen wir nur den Dipol ın f, um den Winkel a zu drehen, wenn in der z-Ebe:i 
die Anströmgeschwindigkeit mit der .-Achse den Winkel a einschließen soll (vgl. S. 1" 
Anm. 2). Der Anstellwinkel des Flügels ist dann a,+«a. Jetzt wird natürlich auch © 
/\rkulation eine andere: sie ist wieder so zu bestimmen, daß der eine Staupunkt der Strömu! 
in der Ebene ın das Bild der Hinterkante, also ın f v” fällt. Wır setzen also 


fe’ ıfa (2) P il T () & ! eiunm \ 
i:b —— 
een tn ra tn] 
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. d m u. 2) ® ry* ® 
J erhalten aus der Bedingung t 0 für # v” die Zirkulation 
(i 
| (- sin? ua 
I y ; zz } 
— COS a — sın a . En Ger VE Ua 
en 7 sın ur 
Pa N Das ® 
or andere Staupunkt fällt dann nach 
w.J iu \8 
| | | ; sıIn? u 7 cos u 7 
er 1 ' 
| = 
r 7 
(14), 
() ar \ 5 Y 2 , z 
1 sin utcosa- sin a[eos UN— | | sın? u a) 
. ’ „T : ur 
. | <zx° für a>O 
1 - sin UTCOSsa— sin alcos ur | | sın? u 1) 
7T e „T 3 








lieet also für «>0O auf der Unterseite des Flügels dieht bei der Vorderkante. Setzen wir 


























(10). / ee!" gr, so erhalten wir bis auf Terme in 7? ?), 
/ 7 2» 1 / I ' 9’ 
(d w\? ae rt - ae a y\ 
| \= Vz e?iut " e-2ie‘1+2re iun|, e!® + - +7 E ien|, Br" | 
dz \ In I, EA; 7 r 
(iv), 
P yınm 7T !e it I 2 ’ ‘) , ’ | 
.- . _ (e: I } er - P a - ‘ t 4’ P Y;j IT | e! dt t ‘ us 
sin? u sın un \2n 1 7 an K; ; 
(11). 
(dl: l I, P i@y Ä ö l | ar . | E Jun / ve iu | 
N ir etun — I — T—e-tr "7? — A— Peer 0,0, 
ührte, et V. | 7 | Isın’un 2asin un ) 
dz Me-:iom 7 h y | 8-3 pn iye-tun| | 
; - I7T e?2iun | er iun L7?| zu / Pe 
dt Vz rs 7 14 sin’un 2nsinun|. e 
(le), 
ie Yyınnc R e-?iun | 
1S Ist, . | +7 3: "a Sr nee 
| 2sın un ksın? un | 
(d.h 6 
! T P ion f T u e ’ un + ... 
| 7 
(Sa ür Auftrieb und Moment haben wir nunmehr 
} ,‚o [(dw\? ‚o [/(dw\dz 
rößen P.—-iP,„=i- | dz ! | At, 
2 „ 2 dz 2 dz} dt 
ves 
Oo dw\? Oo dw.” dz 
M = N & z| dz S R $ | 2 At"), 
st. dz Br dz At 
Die Auswertung ergibt 
lung 
We 6) „\2 u 
hneı IS PP > P -. ‘(5 1 «) r2 2b e / ER ) Im . = 
« ttiP,=e \? eiiz sınay 1 sin? ua eosasın ut, 
ebeı - @\V„sinun We Br 
1!) 
he 27 
mu! °) Für — 0 wird! du eta0, also "„e *"rdadureh ist der obige Ansatz für «= nachgeprüft. Der 
d 2 
I 
ind, warum f ” nieht nach T entwiekelt wird, siehe Anm. 3. 
(1: 2; 


”») Beim Einsetzen von (15) und (17) hebt sich ” herans. 
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1 / AR 2 | »\? | 
M zoV; — sin acosa[1-+2cos u | | | sin? u n\ 
E | \ „Sn ur) | ’ 77 j 


> 
“rs - 


er I 

’ j ie ’ ) Zune | 
sinuncosun 2 sın’ asın u ıt- | | | sın? un 
‚sin, Ä = Ä - Pe 





Es ist bei dem Moment zu beachten, daß der Bezugspunkt in € liegt, also in dem Ansat 
punkt des Ruders an den Flügel. Für $ folgt aus (13) und (15) der Kutta-Joukowskisel 
Satz p 


Wir erkennen, daß die Zirkulation um den Flügel in der z-Ebene bis auf den Faktor b mi 
der Zirkulation, die wir mit Hilfe der t-Ebene erhielten, übereinstimmt: /’ ist dimensionslos. 
/’b hat die Dimension einer Zirkulation und gibt zusammen mit oV, die Dimension einer 
Kraft/cem, wie es sein muß. 


Die Berechnung von ® läßt sich in etwas anderer Form einfacher leisten wie folgt: 


Es ıst 
en ie ‚o ((dw\dz ‚o ((dw\”dt 
Pr=FP,—iP, 5415] ge re 


oder nach (S) und (12) mit ei" —f 


- 


Ban .o h u ei ff Er» 
Re ı#?| — 1) 21 —-cosurn +— +1] — tl — t,)? 
| 77 7 T | 
| R ‚„.t-Tsin un P 
je-tunr a) —t” —etur—e)F —L)”-+ = (t —t,)(t iz 


Man erkennt, daß die beim Ausquadrieren der eckigen Klammer erscheinenden Glieder nur 
zum Teil einen Beitrag zum Residuum an der Stelle #, liefern. Einen Beitrag liefert nur 
feier ce)(t - #)Y und das in /' lineare Glied. Nachdem aber die Kraft nach dem Kutta- 
Joukowskischen Satz sich proportional 7" ergeben muß, kann nur das ın /' lineare Glied 
einen Beitrag liefern (man überzeugt sich leicht, daß das andere Glied tatsächlich das Resi- 
duum Null ergibt). Es ıst dann 


n obV,eiweiwz-Jf fg 1 FE 
ER 2 d | dt. 
a / « ! Fr | Y l ME; ’ . RER ’ | l 
z “ii COS u T 7 ar 


Hier ist das Residuum aber sofort abzulesen. Für die Berechnung des Momentes ergıbt 
sieh nach dieser Methode kein Vorteil gegenüber der Reihenentwieklung nach (15) und (16). 


Zur Kontrolle berechnen wir noch den Wert der Ausdrücke (18) für y=0, also für 
verschwindenden Ruderausschlag. Man erhält 
? . 1 eosun 


A=n oV; t,sına: M noV: (2cos um +1l)sinacosa: 4 2 
Die Luftkraft wird gleich der des ebenen Flügels und greift in ein Viertel Flügeltiefe von 
der Vorderkante an, wie es sein muß (?2cos urn +13 —44). Vernachlässigen wir wieder 


Quadrate von -—- , so bekommen wir 
JT 





v2 |» Y . | : | COS UT (Kky 

I=roV,t,|sına cos asın u T|; / - U} 
- L s - l | .) ‘ nr 
N) .T | / 
(l 
.o. A | . f .) Y *- . .) Y . 62) er 

M zoV, sin a cos all +2c0s u) SINUTCOSs ur — 2 sın’asın un 

Nur 1 7ı ‘ ' be ‘ 


Der dritte Term in der Klammer für das Moment ist im Bereich normaler Anstellwinkel 


klein. daß man ıhn ruhize fortlassen kann. a) 
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Bestimmung der Luftkräfte auf einen ebenen Tragflügel 
mit (Juerruder. 
Von @. Ellenberger in Berlin'). 


ruders bei Flügelschwingungen kennenzulernen. Dazu ist vor allem die Kenntnis der 

ıftkräfte notwendig, die bei einer Schwingung des Flügels auf diesen wirken. Zur Be- 

mmung der Luftkräfte auf einen Körper unter Berücksichtigung der sich ablösenden 
Virbel bei nicht stationärer ebener Strömung sind zwei Methoden anwendbar. 


)‘ vorliegende Arbeit ist aus dem Bedürfnis heraus entstanden, den Einfluß des Quer- 


[. Methode von H. Wagner?) Dieser bestimmt unter Zuhilfenahme der konformen 
Abbildung die sich ablösenden Wirbel aus der Bedingung elatten Abflusses an der 
Hinterkante des Flügels. Es ergeben sich so die Kräfte auf einen ebenen Tragflügel, 
der die allgemeinste Bewegung ausführt. 


3, Methode von W. Ackermann und W.Birnbaum’). Diese ersetzen den Flügel 
durch eine tragende Wirbelfläche und gelangen unter Verzicht auf die Hilfsmittel 
der konformen Abbildung zu gewissen Ansätzen für die Luftkräfte. In einer zweiten 
Arbeit wurden diese Ansätze von W. Birnba um’) auf einen schwingenden Trag- 
flügel angewendet. 


Im folgenden wird die erste Methode nach H. Wagner benutzt. Es werden die Luft- 
kräfte auf einen ebenen Tragflügel mit Querruder berechnet, für den Fall, daß der Flügel die 
alleemeinste Bewegung ausführt. In Abschnitt I wird das Problem und die Berechnungs- 
methode auseinandergesetzt. Die Abschnitte II bis VI enthalten die mathematische Durch- 
führung der Berechnung der Luftkräfte. In Abschnitt VII ıst das Ergebnis gebracht mit 
einer Zusammenstellung der Voraussetzungen und mit dem für die praktische Verwendung der 
luftkräfte notwendigen Kurven- und Zahlenmaterial (das Ergebnis kann natürlich auch als 
(‚rundlage dienen zur Berechnung der Luftkräfte auf Höhen- und Seitenleitwerke). Die ge- 
Iundenen Ansätze für die Luftkräfte werden dann noch angewendet: 


l. Auf die allgemeinste Bewegung eines Flügels ohne (uerruder, 


2. auf die stationäre Bewegung eines Flügels mit Querruder, 
3. auf einen harmonisch schwingenden Flügel mit Querruder. 


I. Das Problem und die Berechnungsmethode. 


Das Profil des Flügels mit Querruder ersetzen wir durch zwei gegeneinander abge- 
winkelte Strecken (Länge 1, und I,, Abb. 1). Der Flügel ist senkrecht zur Ebene der Abb. 1 
unendlieh ausgedehnt zu denken. Er ist geometrisch durch folgende Daten bestimmt: 








' Y 
2-Ebene 
ke ———— l y — e] 
L er — 
K ß \ N ’ 





Abb. 1. Angaben für die geometrische Form und für die Bewegung des Flügels. 
Gesamtflügeltiefe 1: 1=1,-+1,, Verhältnis A von Querrudertiefe zu Gesamtflügeltiefe: 
'uerruderausschlag: y. 


ı Mitteilung aus dem Flugteehnisehen Institut der Teehnisehen Hocehsehule Berlin (Prof. Dr.-Ing. H. Wagner). 
°, H. Wagner: Dynamischer Auftrieb von Tragfllügeln, ZAMM 1925, S. 17 bis 35. 


9 W.Birnbaum: Die tragende Wirbelfläche als Hilfsmittel zur Behandlung des ebenen Problems der Trag- 
zeltheorie, ZAMM 1933, 8. 290 bis 297. 


°), W. Birnbaum: ZAMM 1994, S. 378 bis 202. 
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Wir legen also für unsere Betrachtungen ein symmetrisches Profil zugrunde, das w 
dureh sein Profilskelett ersetzen. Für Profile mit gekrümmtem Profilskelett läßt sieh d 
Erweiterung der Ergebnisse leieht durchführen (siehe später Abschnitt VI. 


Der alleemeinste Bewerungszustand des Flügels seı dargestellt durch Überlagerung vo 


l. einer Translation des Flügels, gegeben durch die Geschwindigkeit U (Komp: 
nenten V,, V,) des Punktes O, 

2. einer Drehung des ganzen Flügels um O0 mit der Winkelgeschwindigkeit », (eni 
gegen dem Uhrzeigersinn positiv gerechnet), 

3. einer Drehung des Querruders (1,) um O relativ zum Hauptflügel (1,) mit der Winke 
geschwindigkeit @, (entgegen dem Uhrzeigersinn positiv gerechnet). 


Das .r, „System sei ein mit dem Hauptflügel (2,) fest verbundenes Koordinatensystem. 


Die Flüssiekeit (Luft), in der sich der Flügel bewegt, denken wir uns reibungsfrei. 
imkompressibel und von konstanter Dichte o. Für die Bestimmung der Strömung können 
wir die allgemeinen Methoden der zweidimensionalen Potentialtheorie (speziell der Funktionen 
theorie) anwenden. Es gilt dann für die Strömungsgeschwindigkeit d mit dem Strömungs- 
potential 


vb grady, divv0, 


Wir setzen ferner voraus, daß die Flüssigkeit im Unendliehen ruht und daß für sehr große 
Abstände d vom Kniekpunkt 9 gilt 


BT... | 
2] E55 er EEE EEE 


wo ce eine Konstante ist. Die Bedingung (1) garantiert uns bis auf eine unwesentliche Kon- 
stante die Eindeutigkeit des Potentials 9°). 


Wir lassen uns nun von folgendem Gedankengang leiten: Nach H. Waener (Anm. ?) 


ist die Kraft PB auf einen bewegten Körper unter den zuletzt gemachten Voraussetzungen 
eereben dureh den Ausdruck 


\ d £ 
IS 07,19 dt ne Te der ee a: Sie (2), 
F L Fr 


wobei das Integral über die Oberfläche F des Körpers und über Flächen F’ zu erstrecken ist, 
welehe den Raum (im ebenen Falle die Ebene) einfach zusammenhängend machen. Der 
mehrfach zusammenhängende Raum (bzw. die mehrfach zusammenhängende Ebene) ist dabei 
zebildet dureh den Außenraum, der vom Körper und den vorhandenen Wirbelfäden bestimnit 
wird. Der Vektor df des Flächenelementes ist vom betrachteten Flüssiekeitsraum nach außen 
gerichtet: 4 bedeutet den Differentialquotienten nach der Zeit r. Für das Moment WM gilt 


{ 


ein entsprechender Ausdruck 


\ d £& +.) 
m 0 4,\9 rt, Ar] Be ee re Re SP Er > () 
F r Fr 


v ist der Radiusvektor von einem ortsfesten Raumpunkt nach dem Flächenelement df. Wı 
werden im folgenden als Momentenpunkt denjenigen raumfesten Punkt wählen, der in den 
betrachteten Zeitpunkt mit dem Kniekpunkt 9 des Flügels zusammenfällt. Für unser ebene: 
Problem ergeben sieh P und WE pro Zentimeter Flügelbreite (senkrecht zur Ebene der Abb. | 


Für die Bereehnung von S und W® muß das Potential y bekannt sein. Wir zerlegeı 
es in vier Teilpotentiale 9,. 95. 9, und g,. Dabei bedeutet: y, das Potential der Strömung 
die einer Translation des Flügels mit der Geschwindigkeit B entspricht; 9, das Potential d« 
Strömung, die einer Drehung des Flügels um 0 mit der Winkelgeschwindigkeit &, entspricht 
y, das Potential der Strömung, die einer Drehung des Querruders um O0 relativ zum Haup! 


flügel (2,) mit der Winkelgeschwindigkeit @, 2, entspricht: 9, das Potential der Strömung 


Obire Voraussetzungen sind im folgenden immer erfüllt. Es ist dann die Zirkulation längs einer genlge! 
eroßen den Flizrel (und die Wirbel) umgebenden Kurve Null. 
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hei ruhendem Flügel dem Vorhandensein der Wirbel hinter dem Flügel und der Zirkulation 
den Flügel entsprieht (vel. H. Wagner, S. 35, siehe Anm. 2). Die Zusammensetzung 
vier Teilströmungen gibt gerade die Strömung, die dem allzemeinsten Bewegzungszustand 
Flügels entspricht. Unsere Aufgabe besteht darin, das Potential 9 zu bestimmen und 
ııs nach (2) und (3)°) die Kräfte auf den Flügel zu berechnen. 


II. Die grundlegende konforme Abbildung. 


Bevor wir an die Berechnung der Teilpotentiale 9,, %s. 95; 9, gehen, brauchen wir die 
forme Abbildung des von dem Flügelprofil (Abb. 2) bestimmten Außenbereichs der z-Ebene 
»+iy) auf eine Halbebene der komplexen Veränderlichen #. Diese Abbildung wurde 


F/ 
ZLDEene 


EN 


NIS 





- 


t- Ebene 


u —_ — 
— - 


7) 








Abb. 2. Der Zusammenhang zwischen der z-, f- und z*-Ebene. 


vor einiger Zeit von U. Schmieden, Rostock, gefunden. Die Einzelheiten der konformen 
Abbildung sind in der vorhergehenden Arbeit von ©. Schmieden (S. 193) mitgeteilt. Wir 
brinzen hier zusammenfassend das Ergebnis. Die konforme Abbildung des von dem Flügel- 


nr 


profil bestimmten Außenbereichs der z-Ebene auf die obere f-Halbebene lautet für . 3: 


l 


} 


7 


u. 1 - 
sın? un-t 


:=—fl, A, 


2tcosun-+| 


u - 4 (In t +iaretb) E Mr. i ‚ es r : i 
wober # e" mit D<aretnr; f=fli,y) und «== u(/,y) sind bestimmte reelle 
unktionen von A und y. Für die praktische Verwendung der Formeln für die Kräfte 


, die klein sind, also - I. Für diesen Fall der 


interessieren nur Querruderausschläge ; z 


7 
ar \«- 


„ersten Näherung* wird, wenn wir (- und höhere Potenzen vernachlässigen, 
| 


e. . 1 ecosun ’ 
fi, y)=1 und 4 : a et 


» 
_ 


mit VD u—|€ 


I — 


D) 
ap \*- 


er Ausdruck „erste Näherung“ soll im folgenden immer bedeuten, daß wir f und höhere 
i z JT 


’otenzen von p vernachlässigen). Wir brauchen später noch die Abbildung der oberen 


| * 7 > .*. . . . . ’ . ’ 
'lalbebene auf das Äußere des Einheitskreises einer z’-Ebene (2’=."-+iy', Abb. 2) 


Bi e IiUun 
Me t ein 7 . R ‘ R R R . R . 6 . ; . ’ . A (6). 


und e=i“” sind die Nullstellen der Nennerfunktion von (4). Wir haben dann folgendes 
-prechen (in Abb. 2 sind die entsprechenden #-Werte auch in der z- und z’-Ebene an- 
hrieben): 


6) Die im Taxt vorkommenden in Kiammern ersetzte Zahlen sind Hinweise auf Gleichungen und Beziehungen, 
reits entwickelt wurden. 
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dz : | 
“ dt I 
; = 
f, l,e! 0 "; 
[IB (0 0 e ?iun 
! I, 0 # 
X 0 v0 | 
e! c 3 X 
ei = X v0 
I | EN us y 
/ | | sıIn” ua cos ur 
1 | 7 / 7 | | 
7 i . . o ‘ . . (4 
1 | a» \3 ‚ ar 
R | | | - sSINn’ u 7 — Cos ur 
- Y \ 7ı D zT N 
| 
T 
sı IL IT ® 117 
d 2 1 ! | 4 117 > ” | ve ' 
dz’ / | in da 2 4’ r a’? 
7 io 2 r | 
ee 57 / | ’ | 
1» f #2 | | 2t—- cosun-+ +1 
(2 fl in? | TI T 7 
sın ur r z 
dt (f? 2tcos ur + 1) 





III. Das Potential Y der Strömung. 


Ws ıst q Y+92 +93 +94 91, 92, 9, werden für e <] berechnet, und zwar so, 
daß die Zirkulation um den Flügel jeweils Null ist. Durch 9, berücksichtigen wir die tat- 
sächlieh um den Flügel vorhandene Zirkulation (y, wird nur für y=0 ermittelt). 


1. Berechnung von y,. (Gesucht ıst das Potential für die Strömung um den Flügel, der 
sich mit der Translationszeschwindigekeit 3 (V, Veospß; V, -Vsin 5: V=|8]|, siehe 
Abb. 1) bewegt. Wir gehen von dem komplexen Strömungspotential w,=gy,-+tiy, aus und 
bestimmen zunächst y, aus der hier für die Strömungsfunktion y, bekannten Grenzbedingung. 
Ks gilt nämlich längs der Kontur AK des Flügels für die bekannte Normalgeschwindigkeit 


7, 
(Abb. 1) 
On 


(ds 


m, F (“r .YV e Oyık 
on ed ds) a. K 


\ 08 
Daraus folet durch Integration längs s, wenn wir eine belanglose Konstante unterdrücken: 
y’ K p „ V, LK . . . . . . . . . . . . . . . . . D (9). 


Es ergibt sich Jetzt folgende Randwertaufgabe. Gesucht ist eine harmonische Funktion yr, (a, ') 
der z-Ebene, die auf dem Rand (Flügelprofil) die vorgegebenen Werte yıx (8) annimnıit. 
Diese Funktion ist eindeutig bestimmt. y, ist vermöge der konformen Abbildungen (4) und (6) 
auch eine harmonische Funktion von =’ und „’ (Abb. 2). Die zugehörigen Randwerte au! 
dem Einheitskreis der z’-Ebene (Bild des Randes der z-Ebene) sind durch die konformen 
Abbildungen (4) und (6) bekannt. Wir können daher mit Hilfe des Poissonschen Integrals ı 
als Funktion von .«’ und „° ermitteln. Nachdem durch (4) und (6) der durch das Flügelpro!\ 
bestimmte Außenbereich der z-Ebene auf das Äußere des Einheitskreises der z’-Ebene a 
ebildet wird, brauchen wir das Poissonsche Integral für das Äußere des Einheitskreis: 
ls ıst dann 


| \ (r’ Nyıkdd 
| 


r 2reosiq MM An HE Sn ug 


‘) Vgl. z.B. Lamb: Lehrbuch der Hydrodynamik 1931, S. 91. 
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Bedeutung der neu auftretenden Buchstaben ergibt sich aus Abb. 2. Es ist ferner nach (8) 
VIK Vamflze-'N 9. 


Integral (9) nach 9 schreiben wir vermöge der Abbildung (6) um in ein Integral nach t 
eell!). Es ergibt sich nach einigen Umformungen 


fiV (r? — 1)tEIm tr e-iitdt . 
A )asinun +) Hr) HP) —2rcosp-(t?- N) -Arsing-t] wa) 
+o 
ı) t— cos un 


obei zur Abkürzung FF =ctg, gesetzt wurde. Wir betrachten zunächst das 


teeral P, 


sin un 


fiVv (r? — 1)t-tTe-iß dt 
P,= u. s » y} 2) [4 E ro ".. ’ . (10). 
„zasınar (+) Hr) (IHR?) 2rcosp-(f? 1) Arsingp-t] 
+ 2 
ls ist dann 
en 2 ed ch ee ir er rl 


j ' | 
Für >»=(>x) verhält sich der Integrand von (10) wie „. d.h. wir können das 


DE 


<Y. Wir erhalten so 


Integral P, durch ein Umlaufintegral ersetzen, da - 


P,=-2riRd(!=1)+Redl=t)) :- © © :» 2 2.0.0. RR). 


2rsinpg-+itr —1) 
l+r" - 2rcosg 
gelegenen Nullstellen der Nennerfunktion von (10). Die praktische Durchrechnung ergibt für 


ii wir, = sind die beiden in der oberen ’-Halbebene (f- Halbebene) 


ı 


w, unter Beachtung von (11) und (12) 
fIV e-ikr+ru—P) n 4 
mw, —«e | . „’ pilur+yu _L. \ sn ıiß» 
‚— Pf, rtıYy, 4 Zu: ' Er „’ 2 ct : 
oder nach (6) 
iflVsin unlei«"+re-9 e-iwa+ru—P) i PER ii ' 
u Ye ) / eiun t p-iun Vze-t’ r ) cos (u 7 f- Y P) 
(13). 
Vs; sın(yua pP) tsin(lunmyu pP) en PER; A h ’ 
fIVsın un BE Dfeo 1 Vze Hg, eoslun + yu—P) 
: _ 2tcosun 2 





Deutung von (13): Der erste Bestandteil stellt das Potential einer Strömung um den ruhenden 

lügel dar, wobei die Geschwindigkeit im Unendlichen V ei? ist®). Der zweite Bestandteil 

hebt die Parallelströmung im Unendlichen auf und gibt zusammen mit dem ersten Bestandteil, 
. dw 

5 eos(um-Fyu— pP), unser gesuchtes w,. Für v, Ep 

st die Bedingung (1) erfüllt (Beweis am besten aus der Integraldarstellung (9) für y, in 

Verbindung mit (7)). 





bis auf die belanglose Konstante da 


2. Berechnung von y,. Gesucht ist das Potential für die Strömung um den Flügel, der 
sıch mit der Winkelgeschwindigkeit », um 0 dreht. Wir gehen wie bei der Berechnung 
on q, vor. Es gilt jetzt für die Strömungsfunktion y, des komplexen Strömungspotentials 

7» +iy, längs der Berandung K (Abb. 1) 


OY, | O2 K 


11 
v),'2 : 
on ) 


Os 
"araus folgt, wenn wir eine belanglose Konstante weglassen, 
Dr 


yn K > 9) \ 


_ 


n 

_ 

n 
—_— 


(14). 


o 


v) 


sn bedeutet imaginärer Teil, Ne reeller Teil. 

», Vgl. die vorhergehende Arbeit von C, Sehmieden S. 196 Gl. 12 mit !=V0, 

10), Querstrich bedeutet Übergang zum konjugiert komplexen Wert. 

1, Vel.z.B. Lamb: Lehrbuch der Hydrodynamik 1931, S. 95. 

1?) Das Minuszeichen erseheint hier, weil wir b— + grady setzen, im Gegensatz zul,.amb, der mit v= grady 
ınet, 
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y, als Funktion von .r’ und „y’ ergibt sich wie bei (W) als Poissonsches Integral 


| \ I 1) UK dA 1; 
: 2Zn)il+r’—2rcos(g ) 
0 
Es ist ferner nach (14) und (4) 
.. 9 . 1 a 1 = 
m, || oo,’ sin! un u Wh . (# reellf 
tl’ & ” 0) ‚ De ve .!» 
ds 2 2 (?— 2tcosun +1) 


Das Inteeral (15) nach 9 schreiben wir wieder vermöge der Abbildung (6) um in ein Integral 
nach # (f reell!\. Es ergibt sich nach einigen Umformungen 


FL 


+ 
o,f I? (r? 1) (#?) 7 .cdt 15 
U. . . r Io , r ıo PN . ’ m . ( di) 
: | 2asinuon + PA Hr) HR’) —2rcosg- (fh) Arsing-t] | 
(vel. dazu (9a)). Dieses reelle Integral läßt sich durch komplexe Integration lösen. Zu dem 
Zwecke betrachten wir das folgende Integral P’, über eine in der oberen t-Halbebene ana- 
Ivtische Funktion 


m \ o,ft —l)E Tat a 
, — ® x N [267) ” [2] ro . ’ } 
-— 2asinaoa (HF?) + Fr) +) 2rcosg- (fl) Arsing-f] 
a 2 — 2 ! In!t +iaret) £ e 2 E : 
mit / E ‚VB=aretmn: für t>o(">x) verhält sich der Integrand wie 
| ar 


d.h. wir können das Integral wieder dureh ein Umlaufintegral ersetzen, da “ <1I. 
‚„T 


Wir erhalten 
P.=—Ianinot SE. I) FB FERN: +» * #5 a 


2 1 
, ‚. 2rsnp-+t(r—]) . 
H, i und #, ‚ae. sind die beiden in der oberen f’- Halbebene (f-Halbebene) ge- 


I+r" -2rcosg 


lerxenen Nullstellen der Nennerfunktion von (16) (die Residuen sind für den in (16) auftretenden 
Integranden zu bereehnen). Der Zusammenhang zwischen P, und y, wird durch die Ver- 


i } : b BEN Ed 3+2- . 
schiedenheit der Funktionen (f) und # " bedingt: 
.. I. ( 
Wir zerlegen das Integral (15a) in zwei Teilintegrale J, = \... dt und Jy=\... dt. Es wird 
m L 
dann auf Grund der Festsetzung über are f 
J, e2ii T Jyı Pr . i - f f 5 ö R R s s ; . : j (17). 


Da ./, und ./,, reell sind, lassen sich beide aus (17) bereehnen. Es ist 
Jtd1=v,=RetPı+tgy-ISm{P} . -» : 2 2 200000. 8 


Die Berechnunz von (16a) ereibt 


(!) g* I" . { Il +2 “| : : nY j 
5 rer ch. ie — I)eter er e:9 (1 m Ire-tt ur’ +1 
a A 7 
Re er 
(4 e-tut2re (l+—) . (16l 
7 
} P 21 re 24 7 
e! (} { ! e-!1)( en | ; 
! e ‚ 





Die praktische Aufspaltung von P, in reellen und imaginären Teil ist nicht notwendig. Ma 
kann direkt das komplexe Potential =, 9, +-iy, aus (15) und (165) ablesen: 


13) Die Berechnung von Ps hat Herr von Borbely ausgeführt. 
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1 tgyti)leteieuatere H4isinaz(l Hr - \z’eit H2yru)_42’eile PER 
32 ni \ T 
je i(2u T { 2yu) en pP ilum ayu P iur t 2) u 
‚(tg y—1) u kisinunli+] — (19). 
’ \ 2° IT) z 2 J 
u 
„’ pin _ e-iumz pe 
. r Yzyı a” 3; o ud 
a) er ee 








u dw, 


Gl. (19) gibt w, für - <19) Free Pr ist ddie Bedingung (1) erfüllt. 


3. Berechnung von y,. Gesucht ist das Potential der Strömung um den Flügel, wenn 

Iı das Querruder um 0 mit der Winkelgeschwindigkeit ©, dreht. Die Ermittlung von q, 

bereits durch die Berechnung von 9, mitgeleistet worden. Wir haben jetzt als Rand- 
dingung für die Strömungsfunktion y, des komplexen Potentials mw, = g,+tiy, 


u für 20 
7 a Ds r ” 
| =-19P für t<V®. 


Das der Gl. (15a) entsprechende Poissonsche Integral wird jetzt 


SF. 


. F n 
| a, 7° (r?-- 1)(F) "dt 
Ws‘ , 2rsin un +) Hr) (I +’) 2rcosg -(F? 1) Irsing- |’ 


d.h. y, = Jr, wenn in J; statt », sinngemäß @, geschrieben wird. Aus (17) berechnet 
sich J, zu 
map, 


J, U’s - . . . . . . D . . . . . . . . (20). 


sın 2y 


ar 


Ks ıst dann nach (16b) und (20) für —- </| 











JT 
DR f’ I? | j Be Y 1: 
2 0, =. z?’et@utt2i) -4isinaali+ 1: etlun +27) —_4z’eilzun+2yu) 
FR es 3» sin 2y | f 7 
P-iQur+2yu) u | „\e-i(ur+2yru) e-i@ua+2yu) 
r Fr 4isiın un | 1 - > | “ (21). 
2 x 
2" pin 7 E ium = „' 1 2 
| e?inn (z’ P- Yıu ri m 
2 | 2 
dr, . : i - 
"ür v, “ ist die Bedingung (1) erfüllt. 


dz 


4. Berechnung von y,. 9, ist das Potential, das die von der Hinterkante des Flügels 
sich ablösenden Wirbel erzeugen. Jeder abgelöste Wirbel mit der Zirkulation /" bewirkt eine 


entzerengesetzt gleiche Zirkulation !’' um den Flügel, da die Gesamtzirkulation nach (1) 
tür jeden Zeitpunkt Null sein muß. Wir setzen “- 1 voraus; ferner kleinen Anstell- 
oz s i : £ z , AV. AV 
ıkel 5 des Flügels, d.h. V,, soll klein sein gegenüber V,; ferner seien ,”, !, ds 

dı dr 
E do, . . Bi v . Pr r. . . . 
Be 2 klein. Die Stärke (/') der sich ablösenden Wirbel wird dann auch klein sein. 
v„ AV, dV de do, ö 
dukte der kleinen Größen —-, & *,0,, -, “, 7 untereinander werden 


dt’ dr’ dr’ * dr 
an geeieneter Stelle immer vernachlässigen, da sie klein von zweiter Ordnung sind (wir 


r 
14), Der Fall y „ läßt sich dureh Grenzübergang erledigen. Er interessiert uns im folgenden nicht, da wir 


l voraussetzen. 


y 


r . ... T .. . o. 
15), Die Fälle „= 0 und = „ erhält man dureh Grenzübergang. 
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verweisen später des öfteren auf diese Voraussetzungen unter der Bezeichnung Abschn. II, 
Für die Bestimmung von y, kann nun auf Grund dieser Vernachlässigung y=0 gesetzt werd: 
denn für y=F20 würde sich ein Zusatzglied ergeben, das proportional 7'-y ist. Außerd: 
können wir für den weiteren Verlauf des Flügels von der Eigenbewegung der abgelöst 
Wirbel absehen; denn bei Berücksichtigung der Eigenbewegung der Wirbel würden wir fi 
das Potential y9, ein Zusatzglied erhalten, das von zweiter Ordnung klein ist. Bei der B: 
wegung des Flügels mit der Geschwindigkeit ® bleiben dann die Wirbel im Raume fes 
Nachdem wir kleinen Anstellwinkel 5 vorausgesetzt haben, benutzen wir für die Bereehnuno 
von 9, für die Riehtung von ® die V,-Richtung, indem wir die Wirbel auf der x-Achs« 
annehmen (der Fehler, der dadurch bei der Bestimmung von 9, entsteht, ist von zweit: 
Ordnung klein, also zu vernachlässigen). Durch 9, haben wir schon die Translation d: 
Flügels mit der Geschwindigkeit U berücksichtigt. Zur Ermittlung von y, müssen wir daher 
den Flügel als ruhend voraussetzen und das Potential für eine sich auf der »- Achse von (ı 
Hinterkante aus erstreckende Wirbelschicht berechnen. 

Wir betrachten zunächst einen einzigen Wirbel mit der Zirkulation /' (entgegen dem 
Uhrzeigersinn positiv gezählt) im Abstande x von der Mitte 0’ des Flügels (Abb. 3). Wir 
bilden den Flügel dureh die Funktion 


I 
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Abb. 3. Der Zusammenhang zwischen der Z- und z-Ebene zur Bestimmung von %3. 


auf einen Kreis um 0” mit dem Radius b ab. Bezüglich der Bezeichnungen siehe Abb. > 
(für diese Untersuchung ist vorübergehend 1=2b gesetzt, ferner geht die y- Achse nicht wie 
bisher durch 0, sondern durch 0’). Mit Hilfe der konformen Abbildung (22) läßt sich das 
gesuchte Potential in bekannter Weise angeben. Es ist mit den Bezeichnungen der Abb. > 





!' I’ | 
95 ,\€ &)=5 ev 
oder Be. Een #2 1° GE Des. en Sen Pre 
Li in 2X 
TR ( - iy, =— 17 = 
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mit 


24,m=d Au -P., 2. 8 ra 


(Dabei ist von der Tatsache Gebrauch gemacht, daß die Größe von I’ durch die konform 
Abbildung (22) nieht geändert wird.) Haben wir nieht nur einen einzigen Wirbelfaden, sondern 
eine sich von „eo — b bis «= .e, auf der .r- Achse erstreckende Wirbelschieht mit der Zirkulation 
’—=u(r) de, so wird 





Xo Xo 
i(, Z—-X a 
7 vule)de und m \In umdR . » (255 
e FAR 2ın) Z—-X, 
2 h 
X 
2 ar aaa Er dw, . 
Die Zirkulation um den Flügel ist jetzt -\u(ae)dxe. Für v,- gi, Ist die Bedingung (1) « 
b ie 


füllt, weil die Gesamtzirkulation der Wirbel (die alle im Endlichen liegen) Null ist. 








„ange 
Id Mi 


— 


III, 
'erdı 
‚erd« 
]öst: 
ir für 
er B:. 
e [e: 
hnuno 
-Ach: 
‚weiten 
N dı S 


daher 
on der 


1 den) 


Wir 


(22) 


Abh. 3 
ht wie 
h das 
Abb. 3 


ge 

- 
.. 

—— 


»)».) 


(z2a) 


form 
‚ndern 
ılatıon 


(Di 
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". Die Bedingung für endliche Strömungsgeschwindigkeit an der Hinterkante 
des Flügels. 


Der allgemein vorgegebene Bewegungszustand des Flügels ist bestimmt durch die 
dy 


‘ 


uponenten B(Vx,V,), ®,, 2,1, Die Strömung um den Flügel soll nun so bestimmt 


len, daß ihre Geschwindigkeit an der Hinterkante des Flügels (Hinterkante des (uer- 
ors) endlich ist (die Hinterkante ist auf Grund der konformen Abbildung (4) gegeben durch 
/, vgl. auch (7)). Diese Bedingung muß gestellt werden, damit wir glatten Abfluß von 
Hinterkante erhalten. Unsere Bedingung muß also lauten 





dw dm dw, , dw 
5.3 2 3 A . 
Droot | r + + endlich , 
dz dz dz ° dz j t, 

a ®  ! ; 
> die Variable der Abb. 1 bedeutet. Nachdem aber für ft, sich N © ergibt, können 

(2 ’ 

;r die Bedingung auch schreiben 
dw, , dw, , dw, , dw, 

; I a a ee 


P r 7 
de: de’ d’ dz’h=ı, 


dw, dw, dw, 
de’ de’ dz’ 
,. w,, w, sind nach (13), (19) und (21) als Funktionen von 2’ bekannt. Ferner ist nach (4), 


(7). (22) und (25) für y=0 


& “) (dw, d 

delt, Fr dz)y=ı 
(dabei ist zu beachten, daß das z der Gl. (4) sich von dem z der Gl. (22) durch die additive 
Konstante bcos ur unterscheidet). Es ergibt sich so für (24), wenn wir noch durch den 


\ach Abschnitt III, 4 darf y=0 gesetzt werden bei der Berechnung von 


l od n; 5... 
u lt ze 


) 


helanglosen von Null verschiedenen Faktor © et"” dividieren, in „erster Näherung“ 


y or q,| | | r 1 i 
Vy 2 V.y y +9 (ds +3 x fl == 3 . . . . R (24a). 


ür eine Wirbelschicht ergibt sich entsprechend 


No 


‚Bi ‚9% ‚IfıjecHt 
arg at 9" zı\) vl 


bh 


al)de=0 ... 0. ..(24b). 
J 


(os 4, 9, Sind Funktionen von «ur und sind in Abschnitt VII angegeben. Wir führen zur 
\bkürzunge ein 


.. | ch . 
Bun j .)” 
\ U 7 l | x l / . . . . . . . . . . . . . . . . (2), 
hzw, 
Xo 
Mai 
Kr / » » Ns 
Westin... 2.2.2200 ie) 
bh 


"ine elementare Betrachtung zeigt, daß V,* diejenige Geschwindigkeit ist, welehe die Wirbel 
I bzw. (ae) de) für die auf die Z-Ebene abgebildete Strömung (Abb. 3) an der Stelle Z—b 
\ervorrufen. Die Bedingung (24a) und (24b) besagt, daß die Strömungsgeschwindiekeit der 
I die Z- Ebene transformierten Strömung für denjenigen Punkt Null ist, welcher der Hinter- 
ıte des Flügels entspricht. Die Bedingung endlicher Geschwindigkeit läßt sich durch ent- 
vrechende Wahl von u(x) immer erfüllen. 

(Für den Fall, daß der Wirbel 7’ im Unendlichen liegt. erhält man aus (24a) 


Go +, q,l „I I’ 
y \ xy J od, —+ I: nd, + ==), 


Yy 4 t : | “N 


ıell für stationäre Strömung (9, = w, 0) 


| V, Io Vx Yıdı i= E. 


t 


‘ 
/ 
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Dieselbe Bedingung erhält man aber auch, wenn man für die stationäre Strömung um « 
Flürel dureh Einführung einer Zirkulation / um den Flügel die Bedingung endlie) 
Strömuneseesechwindiekeit an der Hinterkante erfüllt; vgl. dazu die vorhergehende Arbeit ı 
(, Schmieden Gl. (13). Für die entstehende Strömung ıst dann (1) nicht mehr erfüllt, 
der Gegenwirbel (+ 7) ins Unendliche gerückt ist.) 

V. Berechnung der Kraft auf den ganzen Flügel. , 
(Für kleine Ruderausschläge y.) 


(‚esucht ist nach (2) und (3) 


dı z dı\ } 
FR FF 


Entsprechend der Zerlegung des Potentials in y=94,+9»+ 9, + 9, werden wir je vier Teil. 
inteerale erhalten. Wir berechnen zunächst 
Bein U Heidi: + a a 5 Fi eier, 
F+P F-+F' 
F’ ist irgendeine Fläche, hier beim ebenen Problem eine Gerade, welche die Ebene einfach 
zusammenhängend macht. Zu jedem abgelösten Wirbel (Wirbelfaden) gehört eine solehe 
Fläche F’. Wir wählen sie speziell so, daß durch sie der betreffende Wirbel (/") vom be: 
trachteten Flüssiekeitsraum ausgeschieden wird (Abb. 4). 








14 
| z- £bene 
ke l} 1. l y 
| ıdf r O | Z Bi 
Be EN N, -—Wibelschicht — f 
N nn 
Se 


Abb. 4. Zur Bereehnung der Teilintegrale für die Kräfte und Momente. Der Flügel mit Wirbelschieht und Fläche F’ 


7.9, und y, sind längs F’ überall eindeutig; daher ıst 


o,di=0; Iml,411=0 P#=133: .. nr 0 0 M 
F’ F’ 


Nach (26) und (27) ıst dann noch zu berechnen 


N, Yy #7: Mu \q m At] =1,2,3) und t,=ig,df: M,= \y,|t. At] re.) 
F F F+F' F+F' 
Die vektoriellen Integrale (25) wollen wir mit Hilfe der komplexen Bezeichnungsweise 


umschreiben. Es ist (Abb. 4) 


dt i-dz-1: Ir, df]: 2-.dz-1 
(@- Achse, - Achse und positives 2 dz bilden ein Rechtssystem). Ferner nach (15), und wenn 
wir z2= z-e!* schreiben: 
\y,idz Viiz-sin(<— P)idz. 
F F 
Das Inteeral der rechten Seite ist aber Null, auf Grund der einfachen Form unseres Flüge!s. 
Wir erhalten daher 


Y e nn . . n 
N, \y,tdz w,idz. 
F F 
Nach analogen Überlerungeen ist 
M =ipy,2dz=|w, 2daz 
] F 
Kbenso für y, und g,: also 
“,=iw,idzs; M,=|iw,#dz 6=123. 
F F 


Diese Integrale lassen sieh nun dureh Integration im Komplexen auswerten. Dieses soll 
zwei Beispielen (St, und M,) gezeigt werden. Die Berechnung von 8, und M, ist bereits \ 
NH. Wagner (Anm. 2) ausgeführt worden, dessen Ergebnis wir einfach übernehmen wol! 


e \ j d d 
Wir werden jedoch in Abschnitt VI VrBaLLT, und 0 N) nach anderer Methode n 


einmal bereehnen. 





it, da 


r Teil- 


(6), 


infach 
solehe 


m be: 


iche F’ 


(27). 


(28). 


sweise 


Wenn 


.. hi 
lüge|s. 
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1. Berechnung von $t,. Nach Abb. 2 können wir schreiben 


SU Be Be}. 
St, : \ w,idz= \ m, 1 tt $ wÄz dz. 
F + 2’ l 


dz 
dz 
als Funktionen von z’ bekannt. Der letzte Summand in », (21) ergibt nach (4) und (6) 
Funktion von t bzw. z 


ist am zwecekmäßigsten, mit der z’-Variablen zu rechnen. ı«, und ‚ sind nach (21) und 


) + 2 _° 

2 172 sin! ee „2 
»,PFsin’un t 2, 2 (29) 
.) „. 9) Y Di .) . . D) m. ‚) F .) . . . . . . . —r . 

2sın 2} (?— 2tcosun+|1) sın2y 2 
.ıs zugehörige Integral über die Oberfläche F des Flügels liefert demnach keinen Beitrag, 
ıs für die praktische Berechnung von Vorteil ist. Den Rest des Potentials w,, der für das 
Inteeral einen Beitrag liefert, nennen wir w,*. Wir haben dann 


dz ,, 
st, n- d ,“ N 2’ d u = Ir ) Rei (z’ un 0) ! 
TE 


gr l 


2. Berechnung von M,. Nach Abb. 2 ist 


u \ a4 $ zen 2. 
IH. == w.Zzqad2z Ww. 2 dü2. 
3 3 3 d 2’ 

F 2’ l 


Im M, wieder durch ein Residuum zu erhalten, betrachten wir zunächst das Integral 7;: 


Bi a rn r ’ 
P,=\m,2da: Puehde. 
F 2 — | 


Der letzte Summand in dem Ausdruck für w, (21) liefert wıe oben keinen Beitrag für das 
Integral. Der Integrand läßt sich nach (21), (4), (6) und (7) wieder als Funktion von z’ 
schreiben und wir erhalten Ä 


dz . j 
P=Pme, ‚dz: Braun . . HD 


„"' 


| =] 


Der Zusammenhang zwischen P, und M, ergibt sich (analog wie früher zwischen P, und y,, 


vgl. (18)), wenn wir berücksichtigen, daß bei P, unter dem Integral der Faktor t ” [infolge 


dz i F A ! h : 
von 2), bei M, dagegen nicht auftritt, und wenn wir die Festsetzung O<areti=n 


beachten: 
B :Jı eti7 + Jı Be a ee a (30a), 


i M,=J + Jı . . . . . . . . . . . . . . . . . . . . (31), 
wobei 
mt Br #—e-?iun 0 
dz " dz dz N dz 
Jı- f w,2->de = w,2 ,,‚dt Ju=T w,2-—- de = \ w,2 ,,dt; 
3° dz’ a TE ı “de 7 ia 
et Yiun 0 nl 1 + 
2’ = l „’ po 1 


el. dazu (7) und Abb. 2. J; und .J], sind reell; beide lassen sich daher aus (30) und (30a) 
»erechnen, und damit ist M, nach (31) bekannt. 


Nachdem nun die Methode der Berechnung an den zwei Beispielen 8, und M, vor- 
«führt wurde, bringen wir das Ergebnis der praktischen Rechnung von (28) in „erster 
(herung“ und mit Vernachlässigung (Abschn. III, 4) der Produkte kleiner Größen. 


AP... ö zul? ; ö 
st, ne 4 ( | u 9 \ x7) M, 3 4 | Ir | Y + I; \ x Y) 
al’ al’ 


,—=—i y h® M, 4 I®ı rn 6 NE 
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Die M, (vr —=1,2,5) sind auf O bezogen. Die Beiwerte 9, bis q, sind Funktionen von 
d.h. nach (5) Funktionen von /: sie sind in Abschnitt VII zusammeneestellt. Die Werte 


st, und M, übernehmen wir von H. Wagner (vgl. dort S. 30 bis 34, siehe Anm. 2): 
St, \ JE dt | | h? M, = \ VER [t, Ar] > xcyVx" h? 


FF rip 


wobei JM, auf 0° bezogen ist. 
Auf Grund der Ausdrücke (32) und (33) lassen sich jetzt (2) und (3) berechnen. Es ist 


t 
d # d \ v 
1S 01,9 di= 0 er N, I» 
FF 
dı\ ’ Ad | . 
u. 041? [r. dj] = 01, \vydil . . . ee 
F+F' FF 


Wie früher erwähnt, bezieht sich vr auf einen raumfesten Punkt (Kniekpunkt O mit z=0) 
Wir wählen für den Momentenpunkt vorübergehend einen außerhalb des Flügels gelegenen 
Punkt 9 mit dem Radıusvektor R (Abb. 4), und führen die Differentiation in (35) aus: 


-@ | AN | 
" ni \ " d T 9 an] . \ | d T i 4 a) (394 
FF rip 
Nun ist aber: 
N=-TW+T | 
AR dr ze Ya nn. Fe u GR a GE a = 
co, £ 23 +4 ’), 
dr dı ) ,T | 
mit 
if Bi .. 
E I», für den Hauptflügel 
ı!) w ni . 
19,70 für das (uerruder. 


> 

«»,xt hat stets die Riehtune von df. so daß der zweite Bestandteil von (36) keinen Beitrag 
zum zweiten Integral der rechten Seite von (35a) liefert. Lassen wir 9 nach O rücken 
(t,=0V), so wird mit (28), (36) und (35a) 


dl ” v | 
Nur ) \ MW F ’ Sl, 


| Ip 
nn dr y  UIO, y\ - . . . . . . ü ; ; r , , 5 (35 bl 
| | 
Die Zeitabhängiekeit von st, und JM, in (33) ıst durch „we gegeben, wenn wir beachten, dab 
die Wirbel im Raume ruhen und das .r, „Koordinatensystem der Abb. 3 mit dem Flügel fest 


A BE nn . RR 
verbunden ist. Es eilt (= Relativgeschwindiekeit der Wirbel bezüglich des mit 
( T a “ 
dem Flügel fest verbundenen ., »-Systems). Wir erhalten so speziell für die Kraft ®,, die 
von y, herrührt nach (33) 
d hi ir 


v n dı 7 .e e | | ge8 I? 


Dafür können wir unter Benutzung der Identität 


r + h 
| ‚e? I? A | ge® h? 


und nach (25) schreiben 


r 


SS, 4 ‚ p ‚Znb V, + ! 0 V. b | r” h? ' 


Für das zugehörige Moment haben wir entsprechend nach (33) und (35b) 


dd 


AM, en 
- dt 


Vol. I x — 6°) 


o Vz$; (bezogen auf O’N). 








. 
‘ 
Pie) ST 


”);* 
(„hl 


eitrag 
iicken 


(35 b) 


}; dab 
] fest 


s mit 


\,, die 





16, Heft 4 


net 1936 Ellenbereer. Bestimmune der Luftkräfte auf einen ebenen Traetlügel 1] 





eine Wirbelschiecht, die sich mit der Stärke != u (w)dae von zb bis 2 = x 
«hse erstreekt, muß noch nach . integriert werden. Dies liefert 


„ auf der 


X 


la iurtig arız. Tr | vla)dx 
Y, —ı0 V.2nbV,"+io VY.b \ | ge? h? 
bh 
0 
\ o\ x Ih? ulr) dr 
M, \ = , 
’ ) ] ge® h? 
h 


Die praktische Berechnung von (34) und (35) ergibt unter Berücksichtigung der in Ab- 
sitt III. 4 angegebenen Größenordnungen mit Vernachlässieune der Produkte kleiner 
sen, wobei für V,* nach (24b) und (25) 





en y,l 4; | 
Far ’ V,, ” V.y 2 o, «ie . m, 
wesetzt ı1st: 
oral! d V, E m i de ö - 
25 7 | I ar V.„V„-G Vey+6G,[1 Ir G,IV.o, 
x a. SA 
‚do, , 4 ulr)dal 
PT +GI Vo, V.\, he 
h 
oa®|[ 1aV, do, , 1 „de, | 
N) - | G.V”v+G.P um ie r.1? ’+-G a 4 
IN 3 wrr Fa Va y- s ä: r7 Vo, (| dr ‘nf >| re 5; 


Dabei ist das Moment M auf den in drei Viertel der Flügeltiefe T vor der Hinterkante ge- 
evenen Punkt P„ bezogen (s. Abb. 5). Die Beiwerte @, bis @, sind Funktionen von um 
bzw. 7) und sind in Absehnitt VII zusammengestellt. Die Kraft B erscheint als reine Auf- 


triebskraft (Faktor i!). 


VI. Berechnung der Kraft auf das (Juerruder. 
(Für kleine Ruderausschläge y.) 

"ür die praktische Verwendung interessiert nicht nur Kraft und Moment auf den ganzen 
"lügel. sondern auch Kraft BB, und Moment W, für das Querruder R. Für Bz und Ir 
können wir nicht die Beziehungen (2) und (3) benutzen, da diese nur für einen ganzen Körper 
selten. Zur Aufstellung einer Formel für BB, und Wr gehen wır daher direkt von der 
Bernoullischen Gleichung für instationäre Potentialströmung aus: 


p - bD°’ Oo -C (rt). 


(Gesucht ist, wenn FF, die Oberfläche des Querruders bedeutet: 


’ ’ ä 
N e 0 ’ £ ( ‘ e. y m 
157; \pdi= > \v dt e\ 3-46; \cımaj 0). 
FR FR Fi FR | 
Wir führen die zeitliche Änderung des Potentials in bezug auf das mit dem Flügel fest ver- 
| . (Jor\ 
indene «, y-System eım | \ und erhalten 
ee \or/r 
PP a DA ug 
- (Val } / : (Up, D), 
or Alor|r \orfr 
hei 
z \ —» —m 
BEFF Y 4 o, Tot 


(eschwindiekeit eines Punktes von R ist, dessen Lage vom Punkt O aus durch den Radius- 
tor x bestimmt ist. Es wird dann nach einer weiteren Umformung 


. (0 -. i i " 7) E n N k 

DLyR l - Bn,,d)d —\v?df-+ı ! 140). 

1: 0 ri, o\i R yet > \ Aat-q o\4 lör N, ( } 
FR FR Fr Fr 
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Die beiden letzten Integrale der rechten Seite machen für die Berechnung Schwierigkeit 

N \ 

” * . . ( e za - v [2 ” .. * [2 

da sie in g quadratisch sind za, Di, af). Um hier mit möglichst wenig Rech 

( ? u s 
aufwand weiterzukommen, benutzen wir die Voraussetzungen, die auch bei der Aufstellung 
der Formeln (35) und (39) gemacht wurden (vgl. Abschnitt III, 4). Es ist ersichtlich, 4.5 
dann die beiden in g quadratischen Integrale zu vernachlässigen sind, weil in ihnen die Pr». 
dukte der als klein vorausgesetzten Bewegungskomponenten vorkommen’. Das Syml»l 


IO\ 


ö s ( em 
ist eleichwertie dem in den Formeln (2) und (3) auftretenden . Wir erhalten so aııs 


Orr dr 
(40), wenn wir noch für (97,d) bis auf Vernachlässigung der Produkte kleiner Bewegungs 
\ 
.„0c 
komponenten (Abschnitt III, 4) Vx ' schreiben, 
0 x 
dı\ i 09 ,.. 
17, ydi+oV, 2 Fr 
Pa 0 4,\9 Tri *\5« T (41), 
Fr FR 
Analog zeigt man für das Moment 
d r "oO q 
Mr vogodt oV,\ r,d ae ar Go Ger ae eG Ge a ee p4 
h 4 AN! 7 I+e x\) —_ I (42). 
Fr v=konst. FR 


Es ıst wie früher df=i-dz-1; [v,df]=z-dz-1. Für den ersten Summanden von (41) bzw. 
(42) erhalten wir wieder (zunächst ohne Differentiation nach der Zeit) entsprechend der Auf. 
teilung von g je vier Teilintegrale: 

K,r=\mdf bzw Mr=\o,lWdjl =1,233,49). 22222208). 


. « 


FR FR 


Für den zweiten Summanden von (41) bzw. (42) erhalten wir 


N N 

n OO, .r . a Ü IV, r r 

r=oV.\,, df bzw. Mir=oV;\ 3 [nf] v=1234)... 0. (4). 
FR FR 


\tipr können wir noeh umformen. Es ist nämlich, wenn w,=g,-+iy, eines der vier kom- 
plexen Teilpotentiale bedeutet, 


Eh 
\OYy 15 
! dt=V, 
\ 0x j 
FR 
also 
N : / 1° t — @ 
; U4y . , a,, r e | \ 
Eu d ap idz=io V,\mw, A N I 
VO Vx \ Nır ! 0 Vx \ dz 0 Vx \ ’(,_o | 
FR FR 
(wenn f von VO bis x läuft, dann wird die Oberfläche F7, des Querruders durchlaufen, 
vel. (4). Für Mir in (44) ergibt sich entsprechend 
de Fe Date 
.\0o, ’ . wi, I | 4 
ol u. di]=o V, zdz—oV,imw,2 oV,\w,dze ..... (4b) 
S \ 51 ( x\ dz % | ) r M Ü “\ ) 
FR FR FR 


\e,.dz ist die „Komponente von \yp,„df. Die übrigen Ausdrücke in (44a) und (44b) sind 
Fn FR 

bekannt. Für die Berechnung von BP, und DW, zerlegen wir diese nach den vier Teil- 
potentialen unter Verwendung von (41), (42), (43) und (44): 


Üün r 
Por e dr Nu,at sr 
IS EnHn . .;. 2: a Hr 5 a 
, d i * 
% y R 0 d 1 M, R "Tr MM: R 


Die Methode der Berechnung wird für q, (Berechnung St3r und Mar) und 9, (Berechnu 
von Yır und Wur) näher ausgeführt. 


!#) Die an der Vorderkante des Flügels (infolge der dort herrschenden unendlich großen Strömungsgeschwin« 
keit) wirkende Saugkraft ist auch zu vernachlässigen, da sie von zweiter Ordnung klein ist. 








ını$ 


(41), 


(42). 


hzw. 


A u [- 


(44). 


kom: 


(4t a) 


rufen. 


(44). 


Inu! 


wine 
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{. Berechnung von St3r. Nach Abb. 3 können wir schreiben 


2 


; Bu, | | 
ers \r: dr: \midz min Ban 
FR FR 0 aret=n 


» letzte Summand von #, liefert keinen Beitrag zu dem Integral (vgl. (29)). Den Rest von 
der zu dem Integral einen Beitrag liefert, nennen wir wieder m,*. Wir wollen (46) durch 


mplexe Integration erhalten und betrachten deshalb zunächst das Integral P, 


0 — ıo 
„.dz u! [ „.da m 
P=\m, FRAU: \ Fri, ra ee ee 1, an 
| — » arel 7 0 aret +7 


> 


„.dz. u \ a; 
or Integrand w,“ Fr ist durch die Festsetzung — a = are = + zeine in der ganzen t-Ebene 


nalytische Funktion von ft (mit Ausnahme der negativen reellen Achse der #-Ebene). Man 
‚ojet leieht unter Benutzung von (21), (6) und (7), daß der Integrand für >» sich verhält 








wie .  P, läßt sich daher durch ein Umlaufintegral ersetzen. Der Zusammenhang 
41? L 
wisehen P, und St3r ist durch die Definition von arcf gegeben. Es ist 
P =#ıre!!  Rırne-!!"), .» 2 2 a 3 a a a arte +» RR). 
Andererseits Ist 
P,=-2niRllt=ei")+Redlt=e-t""). . .». . 2 2 2.2.47). 


Am einfachsten wird die Berechnung, wenn man mit der z’-Variablen rechnet. Es wird dann 
nach (47a) und (47b) 


P,=S83r2isiny= —27i[Re (2! =) + Ref (2’=0)]. 


Fr 
Fe, 
[7 


‚ zu bereehnen, der nach (7) und (21) als 


? 


’ R a 
Hier sind die Residuen für den Integranden ıv,* iv] 


«\ 


"unktion von 2’ bekannt ist. Die Zusammenhangsverhältnisse ergeben sich aus Abb. 2; dort 
ist der Integrationsweg für P, in der z’-Ebene und in der #-Ebene durch Pfeile angedeutet. 
Die beiden Residuen sind konjugiert imaginär, so daß sich für P, ergibt 
P,=—2ni:2Re {Ref (2’= o)N. 
Rei (z2= x) ist nach dem früheren (Abschnitt V, 1) bekannt. Nun wird schließlich 
> .) 
P, an 


gıR= 3, -- =" Ren). 
2isın Y sın y ' 


2. Berechnung von Mar. Es ist zu ermitteln 


M;rk=\w,*zdz. 
FR 
Die Berechnungsmethode geht analog wie früher (vgl. Abschn. VI, 1). Wir betrachten zu- 
nächst das Integral P, 


. . 
P,=\n,'e 7,dt mit nzarctsHtn. 


er Zusammenhang mit Mar ist 


P,= Mar e?it! — Mare-?ir, 





| P, läßt sich wieder als komplexes Integral auswerten am besten mit Benutzung der 
\artablen 
P,: 2 ri [Ref (2? =) + Ref (z’ = ON] 


dz ’ a 
dem Integranden w,* z da Auch hier sind die beiden Residuen konjugiert imaginär, so 


17, Dieser Kunstgriff der Bereehnung von N. pr aus P3 wurde von Herrn von Borbely vorgeschlagen, 
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daß der Wert von P, nach unseren früheren Berechnungen (Absehnitt V, 2) bekannt ji 
nämlich 


Pe 2 ni-2 Re Re) (2’= ©)) 
und 
Pr EA; ar fe 
Mar —— Ne Rei (z >0)%, 
2isın2y sın 2y ai 4 
3. Berechnung von Yır. Gesucht ist zunächst 
KM i R - ‘ 4 d f | y \ N dz . . . ° . . . . . . . . ( | 
FR FR 


sei hier die -Variable der Abb. 3. Unter Benutzune der unter Abschnitt III. 4 erhalten 
Beziehungen ergibt sieh speziell für die Berandung (Oberfläche F7,) des Querruders 


YA leid: 2 heos ): dz ,sındd»W: VISV 


ou *® 


Für einen einzelnen Wirbelfaden im Abstand .- von der Mitte 0’ des Flügels haben wir daı 
für (48) 


Zr IF ı 
/ I 
\ır id3 \ vsınddu ch \y sin Jdı ib, \2 ısındd), 
.n T an 
I ( Io 
wobeı 
D we cos U) 

ı/) Art cos BP 
4 r l, cos ı) 


(Bezüglich der neu auftretenden Bezeichnungen siehe Abb. 3.) Es ist praktischer, nach 


/ ZU 
Inteerieren: 
iTEIiIE—- LPT vsinırdy I 
Mın —— ; ib cos d,) mit & 
ui en N\ \ı nn 2 k? 2 cos yr? \ 
Nach Inteeration und weiterer Vereinfachune unter Benutzune von (22a) ereibt sich 
Nun ! (77 VEN AG I cos ı),) Zu L bEi . . . . . . ö (48a). 


Nun ist OO = beos d, oder nach (H) OO = beos u a, also #4, = ur. Ferner ist y, gegeben dureh 


I 
h -zcosun 
vos y, (49), 
2 —beosun 
so daß (48a) ereibt 
! & l, COS MIT P „| 
NMıR ! BE Aare cos (ur heosum) — un ] 7 ya Du . ..(48b) 
TU d eos u 7 ; | 
Weiter bereehnet man nach (44a) 
m Ey ie a 
O Dr ri" A 0 oV,J 
’ 4 .Q „X - 
ol d ioV, T, ! (y, T) Br Ay (U 
u -\9z T J x ra m IK 
Fr 
Jetzt ergibt sich, wenn nach (45) Yır der Bestandteil von B7 ist, der von 9, herrührt, 
dl 5 oO En - = 
NS =: #) N ent! . d R } . ; ; . . . ; j (rl 
in 0 7, ARtE x\3 j 
FR 
RS DE, EEE dx 
Die Zeitabhängigkeit von Mı4r Ist dureh = (und yr (@w)) gegeben, wobeı wie früher 
| dr 


zu beachten ıst. Es ereibt sich somit für (51) aus (4Sb), (49) und (50) 


‚oV,I xzua-+bsinun 
Ber R= = ’ 5 , 
Jı ] a h? 


(Für #1 ist das Querruder gleich dem ganzen Flügel, und wir erhalten für (52) in U! 
einstimmung mit H. Wagner (vgl. dort S. 30, siehe Anm. 2) 


x). sy. d / & . / ’ 2 
ri R r - dd T \ a dT - ! Y r ] x" I)? . 








Ne 


Mi 


en 


daı 


IS al. 


lureh 


(4), 


4Sh). 
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Kraft Yır läßt sich mit der Identität (37) und nach (25) umformen: 


. 0 h V. | | 


\ ‚., fur. sinumT 
Juri 1 aun Vy+ 
zt | 


| 
] ze? I? | 


‚en wir nieht einen einzelnen Wirbelfaden, sondern eine von @=b bis «= .r, Sich er- 
“kende Wirbelsehieht mit der Stärke [= u(#)dx, so haben wir noch nach = zu inte- 
ren. Wir erhalten so 


No 
(8 ‚oO h| - >= Jo nlar)d r) 2 
Pın ! ) m VrV, +5,\ | r? DR? (.»>) 
j, 
X 
r % | + h 
V, FI ptw)dw. 
b 


und ,, sind Funktionen von «zz und sind in Abschnitt VII angegeben. 


4. Berechnung von Wur. a) Es ist zuerst zu berechnen 


Maas dil mit [nafl-&de, 
FR 


wenn zunächst 0’ (Abb. 3) Momentenpunkt ist. Wir erhalten dann für einen einzelnen 


Wirbelfaden 


ua () 
re rw. | 
Mır \ vysındcosV dd : \ >asındcosddd) ee 5° 
a” en. 
0" it 7A 


Nach Einführung von y als Integrationsveränderliche, Integration und weiterer Vereinfachung 


unter Benutzung von (22a) 





I’ A _ a > 
M4ır= 5 [(r — y,) (&? — 6? cos? un) -(e un +bsin un)yx° |. . .(ö4a). 
sa 
_ ae d dd 
Nun wird nach (45) mit (54a) und f 
dr dr 
d h oV.f' i iur sinurscos un) -2a(e ua + bsin um) 5) 
> = - F ® ’ nr iss - ; . i +.) . 
% rn >. ki 17T pP.) Var? — 1? (9) 


hb) Es ist weiter zu bereehnen nach (45) 


le. dee 
FR FR 


\uchdem der Ursprung des Koordinatensystems für diese Berechnung 0" ist, erhalten wir 


r \' y . % l Vx l' u 6 
0 \ x \? 34? ug = cos utly, T) x ee‘ eo u A, 
ner ıst nach (48Sb) 
oV.T — La 
oVx \ y,dx : ; Ir—y,)le -besum)—unya® WW .. 0.0. PN. 


FR 


Die Zusammenfassung von a) und b) ergibt so für Wur (45), bezogen auf 0’ aus (99), 
und (57) 


obV,.IT 2#rsin un +b(ua+ sin un cos ua 


)) > . 9 
lep’ I 4 Va: — 0? 


(DI). 
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Y . . Fr 1. . % . ° 
(Für «—=1 erhalten wir wieder das Moment für den ganzen Flügel und in Übereinstimmu 
mit H. Wagner (vgl. dort S. 33 und 34, siehe Anm. 2) 


d\ i no I 
Wır I, — O y: \ vr It, df] 0 \ x I? 9 2 h: . 
F+F' 


Für eine Wirbelschieht, die sieh mit der Zirkulation !=u(rz)dxz von z=bbis 2x, ı 


streckt, ergibt sich mit (37), wenn wir noch unter Benutzung von OO =bcos ur und (53) das 
Moment auf O beziehen: 


Xo 
“L Jıa V, \ ula)da 


% —onb?|e V.\ o P) 
ıR Ü JıyYad!y A Ir | 22 — h? 





. (58a. 
b 
Y,, und 9,, (Funktionen von ur) sind in Abschnitt VII angegeben. 


Wir geben jetzt zusammenfassend das Ergebnis der praktischen Berechnung von (43) 
und (44) in „erster Näherung* und mit Vernachlässigungen nach Abschnitt III, 4: 


Sır= er [9, VYy— 91x 7]; tiri En Ir Vz Vu + 9 Va’ rl; 

Mır 7 [9; V, — 91 Vx }]; Mir="® 7 E [— 9; Vx«Vy+ 9.» rl]; 

Mon i z 9: ori" 1 R Ya @,; 

M»R = Js ®,: M>3R s 7 E YVıo,;; 

Narmi m Js 95} tır=—i n 7 E JoVx ®s; Zr |; )} 
MsRr= T 0; M'r . ; E I VYx Ds; 

Nur di £ kt y,)(a - beos un) — unyar— DR]; Nir=i e £ E rl; 

Mır: 2 In y) le? bF cos’ um) (ve ua -+-bsin un) ya’ — b?| . 

Mir © d. Ela (vr) +unrye? b?| 





Die Mr und Mir (vr —=1,2,3) sind auf O0, Mır und Mir auf 0’ bezogen. Die Beiwerte q, 
bis y,, sind Funktionen von «u und nach (5) also Funktionen von 4. Sie sind in Abschnitt VII 
zusammengestellt. 


Mit (59) ergibt sich jetzt für (41), wobei für V,* nach (24b) und (25) 





( q,l 4,1 

Wer -Ayy - Eu + 

R d | A, f f 2 

gesetzt Ist: 
OT I d V, N ’ r r ‚ r»- ) ) da ' r 
Sr =i. I | Go! dt enValy : GV,’ + Gis 1’ Er G,V.o, 
| HERE 60) 
‚do, GG. lule)de 
h 





und für (42) bezogen auf O 


oT I? z d V, ’ y y y gr Y 5) de, | ‘ r 
Mn n { (1,1 ep r G,,VxVy GV + ER 2 dı +G,1Vx®, 
( Er \d u WM 
„do, j zer \ ula)dr 
b 





Die Beiwerte @,, bis @,, sind in Abschnitt VII zusammengestellt, 








ng 
Mi 


nu 


Be 


\ das 


DSal,. 


(9), 


te 9, 


tVIl 


(60) 





ınd 16, Heft 4 
\ngust 1936 


Ellenberger, Bestimmung der Luftkräfte auf einen ebenen Tragtlügrel 217 





VII. Ergebnis. 


Bevor wir die endgültigen Formeln bringen, wollen wir für die Ausbildung der Unstetig- 
it hinter dem Flügel andere Koordinaten einführen (Abb. 5) (vgl. auch H. Wagner S. 21, 
ahe Anm. 2). Der Flügel habe sich von der Ruhe aus von s, bis s bewegt. Die Wirbel, 





u 3 Pe e | = 
2 ’ L — & "273 
5 MM I. 
2? 0 ı ula)d | 
Fe Se EN 
u ie L_, T2 a uam | 
4587,5] L;— ES; nl ’ | 
nr 2 


Abb. 5. Einführung der Koordinate « für die Ausbildung der Unstetigkeit # (x) hinter dem Flügel. 


lie gebildet wurden und im Raume ruhen, seien durch die Koordinate « in ihrer Lage fest- 
velegt. Es ist dann, wenn wir beachten, daß der Anstellwinkel 5 klein vorausgesetzt wurde 
vcosp=a!)2—=s- -a—+b(s und a nach links positiv gezählt). Unsere Bedingzung (25a) für 
die Stärke der gebildeten Wirbel wird dann mit der Flügeltiefe = 2b 

r 
1 \ vI+s a 


al J Vs-a 
So 


Vy: wlalda .... N 3 


(u (a) entgegen dem Uhrzeigersinn positiv gezählt). (a) ist durch die Integralgleichung (62) 
bestimmt; V,* und damit nach (24b) und (25a) V,,y,@,@, müssen als Funktionen von s 
vereben sein. 


Das Endergebnis wird jetzt nach (38), (39), (60), (61) (mit V, Vv, V„=-Vß,V=|®3|, 
>= Anstellwinkel): 























‚.onl|, d(V p) ui ..: „do > u „do, 
15 4 I dr +4V?Pp G,V”’y+G6G,? dtv, G,F a; 
2 ula)da 
G.IV eo —V 6.) 
r„ Il On pr \ \ y(s a® + l(s me; (be). 
so 
onl’|l d(V pP) de I den, 
ur = 3 a E FI. ı 9 1 =. “ ' ( P) > (4 . 
IK n 4 zZ, G,V’y+G,l ey ‚Yo | ri „IVo, 
We ist hier auf P,„ bezogen; Abb. 5. 
u ol dıV P) a er ‚de 
PR { Gl dr +6,V?’P— 6,V’7+6,®P > 
de, G ; ula)da 
a7 G,.! 2 2 a “y 
GulVo, Fl dr GuIVo, In \ Yıs a’+lis- a) 
. (64). 
en oalP| , ,d(VP) , oa ‚do, j 
iR 4 G,1 dr +G,,V?B G,V’y+G.l dr GV, 
da (i.- ula)da 
G, 2 e.V 
iz dı G„1Vo, In 1: (s — a” + lis - a) 
in ist auf O bezogen. 
Die Stärke u(a) der sich ablösenden Wirbel ist durch (62) bestimmt: 
0 rR I Fi; q,l g;1 1 / I u- 
V"=-Vß+ 4 v7 4 + j @2° 1\] r+z - v(a)da . . . (6P). 
0 


und ®, sind reine Auftriebskräfte. 
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Zur Abkürzung setzen wir u 


Beiwerte 9, DIS 0; 





n dann 


Io (vr, , Hı (7,3 (fs tr, 2 70 
! 1+-Scos’ 
= Go; 9,=0:; g,=05; % 3 
| n er 5 ie 
4, 3. 1) \ | r SC08s° „IH JIsın cos -sın 7 COS» 
>) .T 
) v 
(Ja 2 7 7 Jo SR 
-) 
PR r [sin 7 1) COS 3 ia Ep 
TEENS E > 
4 3 [7 + sin? ı, 27 sin) cos]: Yıa sın 7; 
| i 
Yin T“ Y sın 3, Js 45 RE 
r |sın? 7 7 sın 1) COS |: Iso (1,, 


sın 7 cos, 


sin y] [1 


Die Beiwerte G, bis @,, sind Funktionen von «ur und daher nach (5) auch Funktion: 
7. Wir habe 


’ sin IR 


.) . x ’ 
=n<osı]|; 


COS n]; 


IS 1, sın cos n]|: 


’ | eos 
/. .) 
’ | . cos R 
(1, ly+sın n]; (1, 2: G,=2008n —1; 
T e & 
(: | ’ sın? „| ( 2 ; 
1 > Sin 3, ycosn a f, > 7 — sın TR cosn]|; 
’ Zr eos 1 
(i SINN 3 1 r cos 2 1; Y.= . 
6 - / I ‘ I). 
ZI Ep ya 
| | r . 2 og > | 
=, FL isın ] 1, COS ’) rn sın 7 COS n — > sın n\2-r cos |; 
- | 
| 2 
‘ | ® R : l Mn 3 . Y .q . N . : 
G=5,.|% SINnNncoSsy+% Sin“ |; ff. I) - sın 1) cos 17]: (,, 
un] +) zı 
| | | sor BE u PR 
5 ® 7 (i, oe sINn?9 TCcos) 75 sn 7 UVu> DB 7SsINn)Ccosı 
ad +) „T 
N .) 
| I | 
(1, — I2(sın 7 1) cos ı)) ) m sin 2] COS n]: 
' „T 
G;. ;» 27 + sin y)(sin ) = eos y) -n(y — sinncosn)]; 
n“” 
Ga=3) ‚u S cos’ y) (sinn eosn —- n)-+6 sin? n cos]; 
| I i sin’ | 
Eros COS N (MACOS sın 9) a-71;5 
& 7 / j e, | 
’ 9 a n Be Rn N cin? . . KA sınZ2 
ef. on el” (1 SCHOSS" N) sin 7 ) ) SIN 7 C0S 177 Isın 1] 
| En = | | 
In ) u Isın 7 1) COS 12 (sın 3, cos) sın 2) COS ]. 
-. 2.38” 








die ım Verlaufe der Untersuchung vorgekommen 
auf die Gleiehuneen. in denen sie zuerst auftreten: 


siehe (24a): 


siehe (32 


siehe (59 


siehe (58a 


siehe (9° 
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IN. Beiwerte sind für den Bereich 0 = u = 1 bestimmt, d. h. nach (5) für den Bereich 0-41. 
Vorzeichen der @ und g sind so gewählt, daß für kleine ur=n G>V0 bzw. g>V gilt. 
Verlauf der Beiwerte @ ist abhängig von „= «7 (und abhängig von 4) in den Abb. 6 
’ 7 : u 
4 für das praktisch interessierende Intervall „=0bis „=, (d.h. 4-0 bis =") dar- 
f Zee - T nenn. una EN .. 
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Gy 
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9]; / 

65 
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h) I ” 2 23 “3 23 a 90 0 9 BB 27 3% 4 5 63 7 81 30 














0 005 00 006 01 0% 02 035 03 035 0* O6 05 0 0008 002 006 01 0% 02 083 03 08 Gr GW 4 
AsaTE 15877 y 
\bb. 6. Das Tiefenverhältnis 4 abhängig von n, Die Bei Abb. 7. Die Beiwerte G,, G3, G>, O19, Gıs abhängig von 
rte G3, Ge, G7, Go, G1o, Gr, Gıa abhängig von n bzw. 4. ' n bzw. 4. 
02 
“| 
gu | 
| ß 
| | 1/7 
os) 
Up 
03) | 
| | 
| | 
| 015} 
| | 
002} | On 
men | 
O125\ 
Do | 
4a) 6 
f ( . or 
| & 
O\ | 
E | Op 
| 
0075) 2 
(32 i | 
. | 
| 
” 
(2) | 0025| 
1 | 
| Q4 
| Fo 
IS | 7 | On 0} ar VEEERERESEN, VEESEIEEEEE, 
j / | 7 
[ g MM I KEHS.H&C„HUOCU 9 ee u u u I ED SE EZ 
EEE nn aan A a SE A 
00 0a 006 01.06 02 0285 03 035 0% 0% 05 0 005 0a 00 Gr 08 2 08 43 035 0% Qws 05 
AEBTE | 
8. Die Beiwerte Gs, G13, G19, Gaı abhängig von Abb. 9. Die Beiwerte Ga, Gıs, G17, O1, Gs0, G55 abhängig 
u n bzw. 4. von n bzw. 4. 
(): 


tellt. Außerdem sind die Zahlenwerte für @, bis @G,, in einer Tabelle I angegeben, wobei 


Intervall 7=V bis = = in 20 gleiche Teile geteilt ist. 
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Wir stellen die Voraussetzuneen. unter denen die zeefundenen Formeln für die Krä! 
und Momente eelten. noch einmal zusammen: Gereben Ist (Abb. 1) der Flügel mit der Gesan 


op Be . ’ | cos n Dr 
tiefe I, mit der Querrudertiefe 1, | Ber = |und mit dem Querruderausschlag | | 
- - 70 
Der Bewegungszustand des Flügels ist bestimmt durch die Winkelgeschwindigkeit o, 
2. > A er dy . 
Hauptflügels ({,) um ® und die Komponenten V, V,‚V, Vsinf, ,=7-, die auf d 
| ar 


mit dem Hauptflügel (1,) festverbundene «=, „System bezogen sind. Der Anstellwinkel 7 ı 
klein vorausgesetzt, d.h. V, ist klein gegenüber V (5 ist der Winkel zwischen dem Hau; 
flügel (Z,) und der Richtung von B; 5 wird vom Hauptflügel aus gemessen, und zwar posit 
dV, dV, do, do 
dt’ de’ "tt de’"e ga 
auf Grund der Inteeralgeleiehung (64) ist dann auch (a) klein. Produkte der Größen, d 


enteeeen dem Uhrzeieersinn): ferner sind klein vorausgesetzt 


ar 


wir klein anzenommen haben |mit Einschluß von - ), sind vernachlässigt. Mist auf P.. WW, 
JT n 


m 


ist auf O bezogen (entgegen dem Uhrzeigersinn positiv); die Kräfte sind nach oben positi 
verechnet. 

Kine zusätzliche (schwache) Krümmung des Profils gibt Zusatzkräfte, die nach Lage un: 
Größe konstant sind, da sie nur von der (zeitlich) konstanten Krümmung abhängen. Für den 
Fall, daß das Protilskelett ein Kreisbogen oder eine Kurve dritter Ordnung ist, erhält man 
die Zusatzkräfte durch Kombination der Formeln (63) und (64 mit den Ergebnissen für 
eekrümmte Profile in stationärer Strömung, wie sie z. B. in Hütte I, 1931, S. 401, angegeben 
sind. Die von der Krümmung herrührenden konstanten Zusatzkräfte spielen bei Flügel- 
schwineuneen und im alleemeinen auch bei Stabilitätsbetrachtungen keine Rolle. 


Anwendungen. 


1. Luftkräfte auf den Flügel ohne Querruder. Die Luftkräfte auf den Flügel ohne Quer- 
' de, 

ruder erhält man aus (695), wenn man y=0 und oe ra 0) setzt. Den Drehpunkt 0 
können wir noch beliebige wählen. z. B. so, daß die Ausbildung der Unstetiekeit hinter dem 
"lüzel nur von V,, abhängt. Nach (65) ist dieses der Fall, wenn q, = 0 ist; das ergibt u" 
Ks ist dann A "Ja, d.h. der Drehpunkt 0 fällt nach dem Punkt 7P’,, der in ein Viertel der 
Flüreltiefe vor der Hinterkante lieert (Abb. 5). Wir erhalten so in Übereinstimmung mit 
HH. Waener (vel. dort S. 55, siehe Anm. 2) 


orl| d(V 9) , ! do 2 Ri nla)da | 
B—=*— 1 = +47°B . v\ ge |, 
t | dr I dr a !Yles—a”’—+l-(s—a)) 

m. erFlt an, Fam 1. | 

I |4 dr 32 dı | | 


Für die Ausbildune der Unstetiekeit u (a) ereibt sich die Inteeraleleiechung 


Ve=;;\] | - u(a)da. 


Su) 


2. Luftkräfte bei stationärer Bewegung des Flügels. Der Sonderfall der stationären DB 
h d(V pP) r fi 
weeung ergibt sich, wenn wir setzen @, =, / 0: ferner müssen die Wirbel s: 
= (ET 
weit entfernt vom Flügel angenommen werden, d.h. =» (siehe Abb. 5). Wir erhalten >o 
aus (65) und (6b), bzw. aus (385). (39). (60) und (61) 


rl ER 
= VUB-G)]: M=77 VG); 
onl RI. 
Bey? 7 V’[@G,,P ER E pn = N V?[@G,;P Gr] 


Die Zirkulation /' um den Flügel ist dabeı 


Partelj 


) 
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Ye ) 16 SEE Ellenberzer, Bestimmung der Luftkräfte auf einen ebenen Tragflügel 22: 
räl 3, Luftkräfte bei harmonischen Schwingungen des Flügels. Die gefundenen Formeln (63) 
an | (64) für die Luftkräfte wollen wir jetzt noch auf den Fall anwenden, daß der Flügel beı 


| ‚stanter Vorwärtsegeschwindiekeit harmonische Schwingungen ausführt, d.h. daß die Be- 
| sungskomponenten V,,, @,, @, periodische Funktionen der Zeit sind. Für konstante Vor- 


d | »tseeschwindigkeit heißt das aber, daß V,, ®,, ®, auch periodische Funktionen des zurück- 
voten Weges s sind. Unsere Aufgabe besteht jetzt darin, u (a) auf Grund der Integral- 
q chung (65) zu bestimmen und daraus das in den Formeln (65) und (66) auftretende 
hı | (eeral .,/ zu berechnen: 
un] i 
sit u (a)da u 
lc . -\ ' (s a)? er | (s a‘ (66). 
d so 


» übrigen Ausdrücke in den zuletzt genannten Formeln sind bei gegebenen V,,y, ®,, @, 
Sy) } . nr . . > 
zur ekt gegeben. Entsprechend der Annahme harmonischer Schwingungen machen wir folgen- 
»n Ansatz für die Bewegung des Flügels (Abb. 10): 














sıtı 
und 1Y, 
den 
man 
lür 
eben 
ügel- 
[A587.10] X 
Juer- Abb. 10. Festlegung der Schwingung des Flügels dureh die Komponenten 1, 9 und y. 
kt 0 
ls seı 
dem 
= 1 Aei’’ Ausschlag y, des Punktes O0 gegenüber einem raumfesten Koordinatensystem 
der (#,.9,)» 
mit ’ B ' , \ 
v=g,+ 9, Prehwinkel 9 des Hauptflügels (l,) bezogen auf das raumfeste Koordinaten- 
Ya = Bei’ a a dq ‚ 
system (#,,9,), wobei gq, zeitlich konstant und N , Ist, 
u RE «T 
y„-y, DPrehwinkel y des Querruders bezogen auf den Flügel (@*, „-Koordinatensystem), 
(eirt dy ’ 
wobei y, zeitlich konstant und L- 0, Ist. 
(T ? 
Die «dimensionslosen, im allgemeinen komplexen, Konstanten A, B,C geben Amplitude und 
’hase, wobei diese Konstanten hier noch unbestimmt bleiben. » ist die Kreisfrequenz der 
Schwingung, r die Zeit. Die physikalische Lösung sei gegeben durch den reellen Teil von 
/.y und y. Unser Ansatz ergibt für ,, eine rein harmonische Schwingung (ein zeitlich 
I Er Y e d „ » . 
konstanter Ausschlag „, hat für das Folgende keine Bedeutung, da nur N maßeebend wird), 
(T 
bi 


y bzw. y eine Schwingung um den zeitlich konstanten Winkel 9, bzw. y,. Zu unseren 
seht 'nüheren Voraussetzungen machen wir noch die Annahme, daß 9 klein ist. Es ıst dann 
siehe Abb. 10; 9 und 5 sind vom Flügel aus gezählt, und zwar entgegen dem Uhrzeigersinn 
ıtiv gerechnet): 


Ve, =— Veos(g pP): V V. V eos p = V 
j dy, ö (67). 
Yy 2 V(®—P) V, Vp | 


dt 


“ erhalten somit aus (67) für den Anstellwinkel P 


’ 1 dy, 
P: 4 V 


I dy, 
Yo Fr Yı V : 


dt dt 
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Der Anstellwinkel setzt sich also zusammen aus einem zeitlich konstanten Anteil (yp,) und ein« 


| I dy ’ >. 
rein harmonischen Anteil |y, a 7 . Der Zerlegung von /p entspricht nach (67) eiı 
Zerlegung für V, 
day, 
Y, V’n Yo, dr’ 


"ür die Untersuchung der Schwingungen können wir von dem zeitlich konstanten Antı 
’y, und ebenso von der Wirkung des zeitlich konstanten Ruderausschlages y, absehen. 

denn diese erzeben nach (63) (und (64)) nur einen zeitlich konstanten Auftrieb ®,, hervor- 

eerufen durch eine zeitlich konstante Zirkulation 7, um den Flügel, die unabhängig von de: 


Ausbildung der Wirbelschieht « (a) ist. Es eilt für diesen Auftrieb 


S, . 3 V N, 


nit 


| u; 
N=alVl-n+yr)- 


(Wenn das Flügelprofil gekrümmt ist, dann kommt zu ®, infolge der Krümmung noch ein 
weiterer zeitlich konstanter Auftrieb hinzu, der für die Schwingungen auch unwesentlich ist.) 
Bei horizontalem Fluge ist ®, gleich dem Gewicht des Flügels (bzw. Flugzeuges). Unsere 
Integralgleichung (65) wird jetzt mit den rein harmonischen Anteilen von %,, y und y 


In hand Iı/ I Z 
Vor dr 4 Yr | de 4 dı ı\) Er „"la)da = 8). 


Infolge der rein harmonischen Schwingungen des Flügels wird sich hinter demselben eine 
(im Raum ruhende) zeitlich unveränderliche Unstetigkeitsfläche ergeben, die periodisch von 
dem zurückgelegten Weg s abhängt (Abb. 11). Eine rein periodische Abhängigkeit der Un- 
stetiekeit mit konstanter Amplitude längs s können wir nur erwarten, wenn wir von dem 
Anfahrvorgange (für s=s, sei der Flügel in Ruhe) absehen. 








Abb. 11. Die Ausbildung der Unstetigkeit x. («) hinter dem Flügel bei harmonischen Schwingungen des Flügels. 


Wir machen deshalb für u (@) den Ansatz 


2) 
ul) U.e EEE RITTER EEE EEE 7° 
wobeı 
vi 2nl = 
-PE 5 ) EUR ı 3 2 I See |} 
\ l 


Als physikalische Lösung wählen wir Ne tu (a)y, entsprechend dem reellen Teil unseres An- 
satzes für 9,, 9, und y,. U ist so zu bestimmen, daß die Integralgleichung (68) erfüllt i=' 
Diese ergibt jetzt mit unseren Ansätzen für 9.9. Yu (a): 





as na NE 20‘ 
vıA-Sas B.(1 YIomıi )+0-(® -2 mi le 
" 3 Pr te 
i20 
- > I 
-1\] eu en 





\Wır führen die dimensionslose Veränderliche £ ein durch die Substitution 


el=s—a 
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ine erhalten dann für (65a) 
. S So N So 
e 
1 l I 
IS r N ui Y N E 
i2w I i2o | > P U i20 BULDES , 
u Pe \ | IH —1le"as-—e "tler?as. . (68b), 
a | £ a 
N 0 


»j zur Abkürzung 











nt« 4 Ge; ‚Ps 
D=A:.2wi Bi 2 wi 2) cf +2mi°) 
hen, | Lt, 4 | 
vor. ” u Er h 
d tzt ist. Wir betrachten „stationären Schwingungszustand“, d.h. der Flügel ist sehr weit 
w zung 
* y ® S nd. » \ [Y@‘ 
'srnt von der Anfahrstelle ss, zu denken. Es ist also N " sehr groß, und Gl. (6Sb) 
hf 
VD 207 1! 20, (f | li 207 nn 
e oe I ı\/, (wo) + e et ar 4 
T fe ) Is‘ in rn In 
FI I 
| | 02 Iı/ | EN - 
em «)=\|/ I+  — 1leos2w&ds (so) \ \ I+.- —1|sin2w&de. (72). 
Ist.) S S 
0 0 
ISere 
(|. (71) muß zumindest für den reellen Teil für alle s erfüllt sein, d.h. es muß gelten: 
U l 20°] b 
Sn 2 l — Fi 
” e Er ER) 2 FR 
(69). Re 27 | (8) 
(Gehen wir jetzt wieder zur komplexen Schreibweise über, so können wir den zweiten 
eine Summanden der rechten Seite von (71) weglassen und erhalten für U 
von DV 
, h 1) .. 
n- I uns ; | (4 1). 
dem - (Fr | 
fe lfe+ = 
Damit ist U der Größe und Phase nach bekannt. Ne fu (a)} (69) befriedigt die Integral- 
oleiehung (68), wenn U durch (74) gegeben ist. s, muß dabei so gewählt werden, daß die 
Bedineune (73) erfüllt ist; dieses ist immer möglich. 
Wir sind jetzt in der Lage, das in (63) und (64) auftretende Integral ./ (66) zu berechnen. 
\lit der dimensionslosen Veränderlicehen £ ergibt sich unter Benutzung von (69) 
I 
ls. R s 5 de iD oo . & x S 
‚ i2w u FE - ,„ i20 i Fan 
J le | = U8e lpe is]. RE N 
> = [« Is : 
‚ ver? 
0 
mt 
JS SF 
(69 cos?2mwm£Ed: sn? on.Sd:s (76) 
Ir), Ve — — q =| == Fa ib). 
' ver? I) VRr+E 
0 0 
(70). Schließlich ergibt sich aus (74) und (75) ” 
Ye is 120 
\ J ıD\ 4 e I (“). 
All rn N Ehe 
Ist I. Ir | 2 „| 
zefundene Wert von J ist in die Formeln (63) und (64) einzuführen, womit die Be- 
nung der Kräfte auf einen harmonisch schwingenden Flügel grundsätzlich erledigt ist. 
Luftkräfte sind rein harmonische Funktionen der Zeit bzw. des zurückgelegten Weges s. 
Bei gegebenen A, B,C ist Phase und Größe von J nach (77) bestimmt durch den Faktor 
Su 
-_ Pet 
JT D\ : ' 4 
f. IE T 2m, 


Beate B[t 20: +c.("420i) 
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4, B und C ergeben sieh aus der Lösung der linearen Differentialzleiehunzen, welch: 
Schwingungen des mit Masse belegten Flügelprofils unter der Einwirkung von elastise 
Kräften und der Luftkräfte beschreiben. Die Untersuchung dieser Flügelsehwingungen w 
hier nicht gebracht, da sie über den Rahmen dieser Arbeit hinausgehen würde. 


Die Integrale (72) und (76) lassen sich in geschlossener Form durch Besselsche Fn 
tionen ausdrücken. Es ergibt sieh für sie nach von Borbely'*) mit den Bezeichnungen 
Jahnke-Emde'’) 


: | 
7 j . \ ur = 
Gele) a sın m +.J, (m) cos + N, (en) . ys(®) sin © N, (wm) +08 w..J, kw) (if 
| 


») 





i JT ; ® j 
le. (w) | Isın m (.J,, (en) N, «))) cos ın (J, (m) -N,(o))|. 
(12a 
Is (m) ug Fr [sin m (J, o»)—+ N,(m)) ! cos (J, (on) N, (o))] 
mit 
> | >| | > 119 
Ne) Y.io) J (o)In a. N (en) Yo) J, (wo) In | 
T 1 | ri 


J, und Y, sind die Besselschen Funktionen nullter Ordnung erster bzw. zweiter Art. Ent- 
sprechend sind J, und Y, die Besselschen Funktionen erster Ordnung erster bzw. zweite 
Art. Der Kurvenverlauf von (72a) und (76a) ist in Abb. 12 aufeezeichnet. Die Zahlenwerte 
sınd in einer Tabelle II zusammengestellt. 








107] 


\bb. 12. Die Funktionen fe, fs» fe und 9 in Abhängigkeit von w. 


Ta be] | e 11. 











PER 0.10 ).>0 0.35 1.50 1.0 1.5 >U0 252 3.0 1.0) AR 
A 1.7805 1.1288 1.1035 1.0467 07979 0 0ETIS 05917 05348 0.4918 0.4299 0.3» 
? 7106 06572 OEINT 05850 05054 04547 0.4181 0.3901 0,3674 | 0.3324 | 0,8 
Te 2.5543 1.9732 1.5556 1.3194 0.9366 0,7594 | 0.6535 05813 05283 04546 0.4 
1.3184 1.1867 1.0535 0050] 1,76 16552 0.5826 0.5292 4880 0.1279 0) 

>, von Borbely: Mathematischer Beitrag zur Theorie der Flügelsehwingungen, ZAMM 1936, S. 1. 


18, Jahnke - Emde: Funktionentafeln mit Formeln und Kurven, 1928. 
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um Problem der Anlaufströmung einer Flüssigkeit im geraden 
Rohr mit Kreisring- und Kreisquerschnitt. 
Von Wilhelm Miller in Aachen. 


1. Einleitung. Wenn eine reibende Flüssigkeit unter der Wirkung eines Druckes durch 

N Kapillarrohr gepreßt wird, so stellt sich bekanntlich die dem stationären Zustand ent- 

‚erhende Geschwindigkeitsverteilung über dem (Querschnitt (im Falle des Kreisrohrs die 

oiseuilleseche Verteilung) erst nach einer gewissen Zeit, bzw. nach Durchlaufung einer 

‚n der Rohrmündung an gerechneten Anlaufstrecke ein, während beim Eintritt in das Rohr 

nähernd die ideale gleichmäßige Verteilung herrscht. Die genaue experimentelle Unter- 

chune dieses Anlaufvorganges, insbesondere der Druckverhältnisse, hat sich vor allem für 

;» Zähigkeitsmessung als notwendig erwiesen, die nur dann auf Genauigkeit Anspruch erheben 
ınn. wenn die Meßstellen des Rohres einen gewissen, von der Reynoldsschen Zahl ab- 

ineieen Abstand von der Einlauf- und Ablaufstelle besitzen, oder wenn durch ein der 

sundlegenden Relation zwischen Geschwindigkeit und Druck beizufügendes Korrekturglied 
»r Tatsache Rechnung getragen wird, daß die Druckkraft die anfänglich erhebliche Grenz- 
-hiehtreibung überwinden und die Beschleunigung auf der Anlaufstreeke erzeugen muß'). 
Die vorliegende Untersuchung, die theoretischer Art ist, bezieht sich nicht eigentlich 

uf diesen wirklichen Anlaufvorgang; sie beschränkt sich auf die Gewinnung einer mathematisch 
<treneen Darstellung der analogen laminaren Überganesströmung als einer exakten Lösung 


ler Stokes-Navierschen Differentialgleichungen, unter der Annahme, daß die reinaxial 


seriehtete Strömung das Rohr in seiner ganzen unendlichen Länge ausfüllt und der Querschnitt 
von zwei konzentrischen Kreisen begrenzt ist. Wir werden dabei das Druckgefälle zunächst 
ls konstant voraussetzen uud in einem Schlußabsehnitt den Ansatz auf den Fall eines zeitlich 
veränderlichen Druckes erweitern. In der mathematischen Durehführung, die insbesondere 
von der Theorie der Fourier-Besselschen Reihen Gebrauch macht, kann ich mich auf 
mehrere andere Arbeiten berufen, die ähnliche Methoden verwenden, insbesondere auf die 
Abhandlungen von P. Szymanski und G. Vogelpohl über die instationäre Strömung im 
reisrohr, die hier eine gewisse Verallgemeinerung erfahren sollen. 


2. Die Differentialgleichung und der stationäre Strömungszustand in einem Hohlzylinder. 
Wenn wir rotationssvmmetrische Begrenzung und axial gerichtete Strömung voraussetzen 
und den Abstand von der Achse mit r, den Reibungskoeffizienten mit zu und die kinematische 


/ähiekeit mit » bezeichnen, so gilt bei einem Druckgradienten für die Geschwindiekeit 


02 
2 " die Gleichung 
0° v Ilo®v Iop \ Ior N 
or roör udz vdL 
2 .. * Y .. > | % pP . ‚ . . 
\Wır wollen uns zunächst auf den Fall beschränken, daß „ einer Konstanten gleich ist, 


ug 
die wir gleich — te setzen werden. Als Querschnitt soll ein Kreisring mit dem inneren 
kadıus b und dem äußeren Radius a angenommen werden. Das Problem möge dann so 
estellt sein, daß zu Beginn der Bewegung die Geschwindigkeit über den ganzen Querschnitt 
ten konstanten Wert U besitzt und daß im Laufe der Zeit der stationäre Zustand sich 
erausbilden soll, der durch das Gleiehgewieht zwischen Druck- und Reibungskraft gekenn- 
sichnet ist. Die allgemeine, von der Zeit f unabhängize Lösung unserer Gleichung 


RE de ee a BE re 6 5777 


diesem stationären Zustand entspricht, hat die Form 
Kern ie aaa ae a ra a + 


ıı die Haftbedingung # = 0 am inneren und äußeren Rand zu erfüllen, haben wir zu setzen 


I’ lo "lot j 2 In 
Ri. 1. 12 / Et a A 
” A lo dA 
= IR > h 


', Vgl. L. Sehiller: Die Entwieklung der laminaren Gesehwindigkeitsverteilunge und ihre Bedeutung für die 
vkeitsmessung, ZAMM 2 (1922), S. 96. 

Ferner Wilh. Müller: Einführung in die Theorie der zähen Flüssigkeiten, Leipzig 10939, 8. 58T. 

Vgel.auech H.Scehliehting: Laminare Kanaleinlaufströmung, Vorträge der Hauptversammlung in Bad Pyrmont 
esellsch. f. a. M. u. M., ZAMM Nr. 6 (1934), S. 20. 
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Man sieht dann, daß der Ausdruck mit b)>V0 in das Poiseuillesche Gesetz für das Kre; 
u ,=c(0’ r”) 


übergeht. Die mittlere Geschwindigkeit « bestimmt sich aus der Beziehung 








dl 
ı 5° .) / .) ER ? r? | 1 2 1 \j« 
rl 9a erdr=2nelgdı 4 +, fr (igr 5 )|, 
h | | 
5 os a 
Nach Winsetzunz der Grenzen und Einführung der «dimensionslosen Größen zeig: 
dl (dl 
ereibt sieh mit 
| Ti 
€ l+ 4 I 
| 
lg 
‘4 
für « der Wert 
ea 
u " 
und damit für die Gesehwindigkeitsfunktion der Ausdruck 
| | | 
(1 o’)le (1-q’)lg ji ie — 
2 ® ‘ 0 a k & ” 
| 0° (1 4”) > 3 
u - | R € ° | 
(1 +")le (1 4") le 
4 








Die ceharakteristische Größe e, die in der Tabelle berechnet ist, wächst vom Werte Null 
für 4 1 bis zum Werte Eins für den Fall g=0 des Kreiszylinders. Wir haben in der Abh. | 
die Geschwindiekeitsverteilung für einige Werte von g bei konstant bleibendem a zur Dar- 
stellung gebracht. Man erkennt, daß die Unsymmetrie über die Ringbreite a — b um so größer 
ist, je kleiner das Verhältnis y wird, bis sich im Falle b=0 die beiden einem Meridianschnitt 
entsprechenden Zweige der Verteilungskurve in der Achse zur Poiseuilleschen Parabel 
zusammenschließen. 


Tabelle für die Größe: 








7 € 
() 1 
0] ).DN 
(1,2 0.443 
0.4 0,24 
0.9 0.17 
6 0.114 
ON 02 
1.0 0) 





Abb. 1. 


3. Die vollständige instationäre Lösung für den Hohlzylinder. Als Lösung der homogen“ 
Gleichung 
l 0 " 


Ir re ee et ae 
vr of 


nehmen wir eine Reihendarstellung von der Form 


v A:2yt, n 
\ A;e ' Z, (4; r) 





Null 
‚bb. 1 
Dar- 
rößer 
-hnitt 
wrabel 
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der jedes Glied in einen Zeitfaktor und einen Ortsfaktor zerlegt ist und Z,(r) eine noch 


her zu bestimmende Zylinderfunktion nullter Ordnung bedeutet. Dann ist jedenfalls 


” 0%, Tv, 

sine Lösung der nieht-homogenen Gl. (1), die unseren Bedingungen angepaßt werden kann. 
Von der Lösung ist zu verlangen, daß sie für alle Werte #=F0 für r=b, r=a verschwindet 
ind für #>0 in die Funktion 

(vv. >= U 2 a Pr a rt 5 Dr De dl DD... En ind Ze 
übergeht. Dann wird unter der Voraussetzung der Konvereenz die resultierende Geschwindie- 
keit in dem Zeitintervall 0 bis t=x einen Übergang von der Anfangsverteilung # = U 
in den stationären Endzustand » = r, darstellen. Als Zylinderfunktion benutzen wir nun eine 
lineare Verbindung der gewöhnlichen Besselschen Funktion .J,(/r) erster Art und der 
Funktion Y,(Ar) nullter Ordnung zweiter Art. Für die allgemeine Funktion Y,@r) gilt die 
Reihendarstellung 





.) | My nr | | 
i er en = r Hr I g J N | Br 
Ya) YA = 















n F. . . . . (6), 
n! 1 fair J,() E n—+2s 
) N | ! \ ( 1)" Fe 
- een n S “ S. ———— sin Ss) r J 
sy) N l 








« 
ww 





in der le — 0.11593 zu setzen Ist und die reehts stehenden Funktionen die gewöhnlichen 
‘ . vr. . . . . . . r 

Funktionen erster Art bezeichnen. Wie man sieht, besitzt die Funktion im Nullpunkt = =0 

eine lorarithmische Singularität, die aber in unserem Falle wegen == /r und Z>0, r>V 

auseeschlossen ist. Um die Haftbedingung zu erfüllen, setzen wir nun?) 







I (4 r) Y.(4 N) S,(Ar) 
J,(+b) Y.(Ab) J,(4b)-Y,(Ab) 





Z, (Ar) = (). 






Die Funktion verschwindet zunächst für r = 5b. Setzen wir aber weiter fest, daß die Para- 
meter / die Wurzeln der Gleichung 
LANE - IRB ERDERB.... 5 2 nn ee ee 
















sind, so ist auch die zweite Randbedingung erfüllt’). Übrigens kann man statt Z,(/r) auch 
die „komplementäre* Funktion 
J/.(ır) Y.(Ar) 


Z,(+r) J.(ie) Y,(Aa) 


















benutzen. Es bleibt dann nur noch übrig, die Koeffizienten A; der Bessel-Fouı ierschen 
Entwicklung der Forderung (5) oder der Beziehung 
x 
U ce [d® r”- flog r| > 4,7, (ir) DE u En U be Se E.V 
il 
anzupassen. Zu dem Zwecke benutzen wir die für jede Zylinderfunktion gültigen Integral- 
lormeln 


; r n r a r |. r n r a a r a r a | 
r Zulzir) Zu Ar r) d dh I Ja FE c Zn+ı(llir) Zulkkr) — Apr Zuldir) Zu (Ar r)| 
. RR k; 





| (10). 
i SL. 2 a .\ | 5 | n°” r Y a r ‚9 n \ 
l Zn (7; Y) dr= 5 7 | 7: s | Zw (4% r) + Zu (A, r) 
& . (4; r) ) 
Setzen wir die Gültigkeit von (8) voraus, so ergibt sich in unserem Falle wegen Z, («) L, (x) 


a 


I» Zur) Zar) dr=0 Mfür i+k) 
) 


a (il), 


Zur (r)dr= a’ Z’ (ira) = 


h 


h? Zr (/, 1) 

















°:) Vel. A.Grav and G. B. Mathews: A Treatise on Bessel Funetions, London 1927, S. 116. 
», Vgl. A. Kalähne: Über die Wurzeln einiger Zylinderfunktionen und gewisser aus ihnen gebildeter 
ıchungen, Zeitsehr. f. Math. u. Phys. 54 (1907), S. >>. H.S. Carlslaw, Conduetion of heat, London 1922, S. 128, 


® h 
Iindet sich der Beweis, daß die Wurzeln von (7) alle reell sind, wenn reell ıst, 
a 
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wo alleemein 
J.(Ar)Y,(Ab)— Y.(Ar)J,(Ab) Sir) 


Zular) J,6b) Y,(Ab) J,(Ab) Y,(Ab) 


einzuführen ist. Multiplizieren wir (9) mit r Z,(/,r) und integrieren über das Inter 
r—hb bis r - a, so erhalten wir die Gleichung 





ac | SEHR 
c(d - r?+flogr)| r 7, @, r) dr= 5 Ar | Z°’ (ira) —b? Z° (2, b)] a ; 


zur Bestimmung der Konstanten A,. Um die Integration auf der linken Seite auszuführ: 
benutzen wir die bekannten Formeln für Besselsche Funktionen?) 


Al” Zu ()] ee 5. ‚(®%), bzw. Ar" . Zu r)] ir" Zn (4 r) 
Damit ereibt sich 


dl 


" | in S Age 
Ir zn) dr; [aZ,Gra) 62) 
. 75 


h 


"erner erhält man dureh partielle Integration 


(dl (dl 

. | . f | | |“ ».) » 
re r)dr: [rra r Z, (4, r) 19" 7,(/rr)| [127,0 r) dr 
“ L.], . | | Al | ıh Ak . 
h) h 
lcı ) | a 

— If? Z, (/;r)\ 

7 Al; | e 


0} ı ) 


Ei u | 
- RN 
Al | | 





| h 


In ähnlicher Weise folgt durch Benutzung von (153) und der Randbedingung 
dl Ad 


I Ca Su 
| r lo: "Z,@;r)dr logr-r Z, (4, r)) | 7,&«r) dr 
r Al, | \h FR 
I 
| | . la | | |“ 
 I\rlogr-Z, (i+«r)| ++; 17, (Ar r) 
Ak | bh A} 


| ja 


ırloer-Z, (4, "| 
) 


+ | | 
1 ıhb 


[ 





In 
Die linke Seite der Gl. (12) geht daher über ın 


, u a Sur SE ee... 2 
(l ed’) +-[aZ, (Axa)—bZ, (A,b)| + =—- Z, (Ara) — 3. (4b —-; 
Al I]; Al ı 


ce 2 . 14 | 1} 
i,? Ks ZEN a) h" Zs (ty h)| 
1; a 


laloga-Z, (/,a) —blogb-Z, (4, b)]. 


Setzt man die unter (2a) angegebenen Werte von d’ und f ein, so heben sich die logarıt! 
mischen Glieder heraus, und es bleibt 


-) 


Ü h BR... un | de ns r 
. Ja 7,0; BZ, (4, b)] ; „ja? Z,0,a) — 6’ 2,(4,b)] 5 A: le ZA) WZ°’(Ar,b)| (dr) 
A], a . _ 


Benutzen wir ferner die bekannte beziehung 


2, (2) + 2Z,(@)=— Z, 


PR 
in Verbindung mit Z,(/,a)=0V und setzen für e den oben gefundenen Wert e= = 
erhalten wir für A, nach Kürzune mit aZ, (/,a) -bZ, (4, b) 
21... Su 
Ir | € (/, a)” 
az, (A, a)+bZ,(4, b) 


\ treatise on Dessel funetions, Cambr. 1922, S. 134 und 1395. 
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nit entsteht aus (4) für die resultierende Geschwindigkeit unter Voraussetzung der Kon- 
senz der Ausdruck 
A 
ad r Su 
a | { (4; a)” 
, / 


Zlr)e- tr 
/;a Z,(1;a)+4;b Z, (4; b) 


) 
| 
| 
| 
j 


; wieder dureh ö ersetzt wird und die Summe über alle Nullstellen der Funktion (8) zu 
° . . . a I u . vf 

treeken ist. Mit den dimensionslosen Größen ae ui 1; Ss, T=7 erhalten 
( dt 


r daraus 


4") lo 


S /,(0s;) e sit 
= (20). 
N j 


87/85 |4,68)+ 44, (as) 
‘ 

"ührt man jetzt statt Z, und Z, wieder die Funktionen S 
snchtet (8) in der Form 


le 
„und 8, bzw. J und Y ein und 

J,(s;) Y.(s;) 

I, (y s;) L, (4 s;) 
nd die weitere Beziehung’) 

J, (2) Y.(.r) J, (a) E. (‚r) 
so wird der Nenner in dem Reihenglied 
s;|S,(s;) + yS,(g s;)| 


Die Gesehwindiekeitsfunktion vereinfacht sich daher zu 


| 


ıd 


(1 “) lv 


| 
1 


or 
- 


T 3:8 
”s 


\+/.(8;) [/ to s;) Yu(qs;) Y.(os;) J,lqs;)l 
R e 
Js) + I lg si) 





\\enn man statt S(os, qs) die Funktion S=- S(os,s) einführt, so gelangt man auf Grund der 
einfach ableitbaren Beziehung 


J,(s;) S,(8;0) Jg Ss) S,(s;0) 


_ 


u «dem gleichen Ausdruck. Ferner ist es auch möglich, an Stelle von Y, die Hankelsche 
"unktion H,=,J,-+-iY, zu benutzen‘). Man sieht dann, daß die entsprechend gebildete 
Verbindung 


J,(s) TA (48) J.(q Ss) i IH, (S) 
oleichen Wurzeln hat wie die oben benutzte Funktion 
J,(s) Y.(g s) J(y s)Y,(s). 


Zur zahlenmäßigen Ausreehnung benötigt man vor allem die Nullstellen der Funktion 
48), (die für einige Werte von q von Kalähne bestimmt und in den Jahnke-Emde- 
en Tabellen wiedergegeben sind. Wir haben den Wert 4-05 zugrunde gelegt und die 
inktionswerte aus den Tabellenwerken von Jahnke-Emde, Havashi und Watson 
nommen, bzw. für größere Argumente (>25) nach den bekannten Näherungsformeln neu 
mmt. Es zeigt sich, daß für die kleineren Zeitwerte (r 0,001 bis 0,005) etwa sechs bis 
on Glieder der Reihe erforderlich sind, um einen bis zur zweiten oder dritten Dezimalen 
auen Wert der Geschwindigkeit zu erhalten, während etwa bei r-- 0,05 bereits ein Glied 
gt. Die Ergebnisse im Falle U=u sind in der beigezebenen Tabelle und in der Abb. ? 


'ergegeben. 


el. G. N. Watson:na.a.0. 8. 77. 
el. 


\ 
° V Il. Ss.Carlslaw: Conduetion of heat, London 1921, S. 17. 
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Tabelle für die Geschwindirkeitsverteilung 
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Abb. 3. 


085 006 
T 


Inn 


der Anlaufströmung in einem Hohlzylinder 








u 
l 1 ® 
n [ — U). 
| ! 7 2 
' Oo 0.59 0.6 0,7 0,8 0,9 0,95 
0) 1 1 1 1 1 1 
0.00 0.693 01,974 1.01 1,03 1,05 0,74 
0.0025 0.565 0.897 1.09 1.093 0,953 0.601 
0.005 0.51 0,845 1.174 1.16 0,868 0,515 
0.01 0.403 0.850 1.27 1.23 0,825 0,476 
0.025 0.52 0.809 1.35 1.29 0.852 0,485 
0.05 0.54 0.051 1,13 1.36 0.807 0.503 
X 0,56 0,970 1.16 1.39 1,0912 0.514 


Zur Berechnung möze noch bemerkt werden, daß für größere Argumente folgende asymptotische 
Darstellung heranzezoren werden kann: 


J.(2) » 


Daraus folet z. D. 


Da ferner 


#40) 2, (4) 


so ereibt sıch 


N (s) a“ 


1 


Fi (q s) m 


y 


cos [r N a) 
| 
| Ber 
(S, (8 J (s)Y,(qs) 


(s) 


FAR 


sın E 


J (a) n 


u TS] 4 


„T 


ng Y, (.r) 


(ZH. 


Daher kann man z.B. für den Nenner in der im Geschwindigkeitsausdruck auftretenden R: 


(12y 7) 


aneenähert setzen 


s[S,()+g S,(qs)] x 
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.) 
ınserem Falle haben wir mit 9=<y> abwechselnd die Werte 0,414 0,264 und 
ya JT 
»414 1,54. Man stellt fest, daß diese Werte bereits bei den ersten Wurzeln (s, = 6,245, 
12.547) mit eroßer Annäherung zutreffen. 
Die graphische Darstellung zeigt, daß sich im ersten Stadium der Bewegung ähnlich 
» beim wirklichen Anlaufvorgang zwei Gebiete ziemlich scharf voneinander abheben, das 
» Wand benachbarte Gebiet, in dem die Geschwindigkeit vom Randwerte Null steil ansteigt, 
— | das mittlere Gebiet, wo das Geschwindigkeitsprofil einen fast geradlinigen Verlauf zeigt, 
allerdings wegen der Ungleichheit der Reibung am inneren und äußeren Rand von der 
erriehtung etwas abweicht. Unter dem Einfluß der Reibung findet ein Ausgleich beider 
— ‚biete statt, derart, daß sieh die rineförmige Grenzschiecht mehr und mehr vergrößert und 
hliießlich über den ganzen (Juerschnitt ausbreitet, während die Geschwindigkeitsverteilung 
nen parabelartigen Charakter annimmt. Dieser Endzustand wird zwar erst nach theoretisch 
nendlieh langer Zeit hergestellt; die Berechnung zeigt aber, daß bereits für 70,05 etwa 
7" der Endgeschwindigkeit erreicht ist. Das würde bei einem Rohrradius @a=0,t em und 
nn R =. a”rt er 
ner Zähigkeitszahl »—= 0,01 em?/sek der Zeit t 0,5sek entsprechen. Die Wirkung der 
| ro 
Reibung spricht sich vor allem darin aus, daß die mittlere Geschwindigkeit nicht mehr 
konstant bleibt, wenn der Anfangeswert U dem Endwert « gleichgesetzt wird. Wir erhalten 
‚len Mittelwert nach der Formel 
a l 
In ’ 2 t 
 — - F [erar ‚[rodo. 
(a? —- 6°), Ki Bu 
h 4 
ine elementare Rechnung ergibt 
SL 
3, Y { N >», >», 
2 3 > . Rn) . . . . . . . —t . 
u LE es?’/ sy? |J,(ys) + Ju lsı)] 
==] 
\lit guter Annäherung läßt sich (28) zufolze dafür setzen 
I 
v 1 | l fi 5 | | 1%] 4 u 
2; > R | + \ led, 3 
u | 4 — € S;j sl 4 | 1’ | 4 
i-1 
Dieses Verhältnis ist in der Abb. 3 ın Abhängigkeit von r dargestellt. Man erkennt, daß die 
Te mittlere Geschwindigkeit im ersten Augenblick der Bewegung, wo die Reibune und daher 
die Verzögerung unendlich groß sind, stark abfällt, sehr bald ein Minimum etwa von der 
Größe O,SS erreicht, um dann allmählich auf den Grenzwert u anzusteigen. 
Für die am äußeren und inneren Zylinderrand wirkende, auf die Längeneinheit des 
| /xlinders bezogene Reibung erhält man mit Benutzung der Beziehungen 
(IH. j 
Q 2 @ 2./,(48) 
SQ»D,(S( 8,18) — 
1 ! T e: zı J,(s) 
he Ausdrücke 
SL. 
N | | a\ Y S J ( Bi > rs;2 
0? ( ( . 4Ss;)« ! 
on Ru =2nua Imu u) 2 . + \ 1 | Ai da 
7 ö 0 Y r—=ıu € : | EEE fi Si, J,(s;) + ./,14 s;) , 
Ss q TE 
n. (30), 
oe ı / (W r ”r Ss J (s;))e= "si 
k;=2nub a Ituu | I 2 \ 1 | —— 
0 Orb Pe N las 
\ } 4 / : 
« als Funktion der Zeit r in der Abb, 3 veranschaulicht sind. Die Reibungszahlen 
I R; o ” . : . r i 
Z, und £; werden für 7-0 unendlich groß, sinken sehr schnell bis auf 
| ruu 4 au 
| ni b ” Aare 1 / 1 —ıq? 
'e Minimalbeträge und erreichen für >» asymptotisch die Grenzwerte — [2 je 
’ n 
BE, 
| Ti Pr . y . | ® od [u ® 
ı FT) also in unserem Falle mit 9 =— die Werte 5,4 und 3,41. 
le — u 
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4. Anlaufströmung im Kreisrohr‘. Wenn der Innenradius des Rohres Null wii 
(b=0,14 0), so geht, da .J,(0) = 1 und Y,(0) logarithmisch unendlich wird, die Funktion ; 
in die gewöhnliche Besselsche Funktion ./, über, Da ferner mit : I der stationäre En: 


zustand durch die Geschwindigkeitsverteilung », =2u(l1- 0°) gegeben ist, so gewinnt ma 
für den ganzen Bewegungsverlauf den Ausdruck 


SI 
? u S J.(4;r) 
v=»>u } ‚+ \ | — le ze (31 
) TE u (1; a)’//,;aJ (1;a) 
ie} 
oder in dimensionsloser Schreibweise 
a ur we S J,(s;0) ’ | (31 
er o” = e ° ya 
Hu | gb u 2) Ni f (S;) ) 


wo die Summe über alle Nullstellen s; der Besselschen Funktion ‚, zu erstrecken ist. 


Tabelle für die Gesehwindigekeitsverteilung der Anlaufströmung im Kreisrohr. 








1 0 (0) 0.2 0.4 0.6 (0.8 0,0 
v 1 | 1 1 1 1 
0.005 1.058 1.058 1.057 1.046 1.00 0.72 
0,0125 1.102 1.101 1.10 1.068 0.861 0.57 
0,025 U 1.20 1.19 1.10 0,74 0,41 
0,05 1.39 1.37 1.206 1.09 0.662 0.375 
0,1 1.62 1.57 1.41 1.11 0.63 0.355 
(1.2 1.86 1.78 1.57 1.20 0.68 0.341 
0, 1.07 1,80 1.65 1.26 0,71 0.38 
x 2,00 1.92 1.68 ‚28 0.72 0,38 


Wir haben auch die diesem Ausdruck entsprechende Verteilung der Gesehwindigkeit in 
Abhängiekeit von r für den besonderen Fall U = « ausgerechnet und in Abb. 4 graphiseh 
veranschaulicht. Die Art des Übereanges in den Poisseuilleschen Endzustand stimmt im 
wesentlichen mit den Verhältnissen beim kreisringartigen (Querschnitt überein. Es ergibt 
sich dabei, daß für 7-05, also bei einem Kreisrohr von a= 0,2 em und einer Zähigkeit 
v—- 001 em?s, d.h. also nach zwei Sekunden die Geschwindigkeit bereits 98.5", des End- 


wertes erreicht hat. Für die mittlere Geschwindiekeit entsteht aus (29) der Ausdruck 


/ "ir S\ 1 
| | \ | .e si“ (32) 
H —— u N; Ss; 
N l 























0 08 00 05 „00 05 090 408 


Abh. 4. Abb. ». 


') Vel. P.Szymanski: Quelgues solutions exactes des eqnations de V’hydrodynamique du fluide visqus 
dans le eas d’un tube evlindrigne, Journal math. pures et appl., Paris 1992, S. 67. (G.Aresi: Sul moto di un fu 
visecoso, Ineompressibile, lungo un tubo eilindrieco, Rend. Senin. Face. Sci Univ. Cagliari, 4 (1954), S. 91-03. 
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rend die mittlere Beschleunigung den Wert hat 


0 dr (U 8 
ı lee ei 


«“ a u 0) 


für U u der Mittelwert für >V0 sich dem Werte « nähert, so erhält man als Neben- 
ıltat folgende Relation zwischen den Nullstellen 


HS L- 
nr: | \l | 
ss? 4’ s' 32 
ır die Abweichung des Mittelwertes gerenüber dem Anfangs- und Endwert U u ergibt sich 
F. 
y l S 
ee Aa 
7 Sj Sj 
ye 


\Wenn die anfängliche und mittlere Endgesehwindigekeit als gleich angenommen werden, so 
zeiet sich auch hier wieder, daß im ersten Stadium der Bewegung eine Verzögerung eintritt, 
‚lie als Äquivalent der Arbeit anzusehen ist, die erforderlieh ist, um die der Wand aufliegende 
(renzsehieht, bzw. die hier ansässige Wirbelenergie zu überwinden. Das Minimum der 
Geschwindigkeit, das etwa die Größe 7 0,53 « hat, entspricht dem Zeitwerte 70,035. Um 
zu erreiehen, daß die mittlere Geschwindigkeit von einem gewissen Augenbliek an ungefähr 


„A 
konstant bleibt, muß das Verhältnis z etwa gleich 1.34 bis 1,36 anzenommen werden. Der 
dem Verhältnis 1,34 entsprechende Verlauf ist in der Abb. 5 dargestellt worden. 


Aus dem örtlichen Mittelwert kann man ferner den mittleren von der Flüssigkeit zurück- 
seleeten Weg berechnen. Es wird 


m 0 
i i BE , 
“ühren wir die Revnoldssche Zahl RN = ein. so haben wir 
® y’ 
FT. 
| Y l = 
| ( 
> aNIT—4 \ | Je e- "8i*) (>34) 
| \g;? g.® 
l eo oı ! 
ie} 
Fa re a?rT E 2 ’ 
Dividiert man durch die Zeit = . so erhält man weiter den zeitlichen Mittelwert der 
v 


(Geschwindiekeit in der Form 


[A T 


"ner ergibt sich für die Stärke des Ringwirbels, der für die Reibung bestimmend ist, 


I 
Ior Il or " \ "J/ (os;) u 
w=, > 20- 5 ;) e . 
20r 2ado da “ "a Eh Ss; 
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Die auf die Längeneinheit des Zylinders bezogene Randreibung wird dann 


FL. 
| v S\ | 
R | Jı u [A E { 7 u H 20 In \ 1 | 2 7 5;2]| 
ı —— si | 
i—1 
y. h wa, ” Au | 
Die Wirbelverteilung e ıst für verschiedene Zeitpunkte in der Abb. 6 dargestellt. 
zeigt einen Übergang von der anfänglichen Verteilung m —=0 für o<1, bzw. w—=» für o 
. P . : . i r m da u R 
in die Poiseuillesche lineare Verteilung 20. Der Koeffizient { der Raı 
. ä Itum 


reibung ist anfanes unendlich groß, nimmt dann sehr schnell ab und erreicht bei T= 00; 
einen Kleinstwert und nähert sich dann asymptotisch dem Grenzwerte 2 


E —777 
| 
| 
1 
| 





02 0» 06 08 10 72 1% 16 18 20 22 du 26 28 30 32 | 
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Abb.,6. 


>. Rohrströmung bei zeitlich veränderlichem Druckgefälle. Zum Zwecke der Anwendung 

auf einige interessante Fragestellunzen wollen wir einen vorbereitenden Ansatz geben für den 

allgemeinen Fall des veränderlicehen Druckes. Indem wir das gesamte Druckgefälle in einen 

konstanten und einen von der Zeit abhängeizen Teil zerleren, soll die Differentialgeleichung ın 

der dimensionslosen Form 

0’ Io" 0% En 

ara ke Fn+,, a een 

zugrunde gelegt werden, wo e'=ca? zu setzen ist, und als (uerschnitt des Rohres zunächst 

wieder ein Kreisrine mit den Radien a und b angenommen werden. Um die Lösung der 

Gleichung zu finden, die sich von (55) dureh das Fehlen der Konstanten unterscheidet, benutzen 
wir die aus den früheren Entwieklungen hervorgehende Identität 


I. 
be \ y S,(o s;) (36 
F(=?2F(tT Beorhl.-n „S u: 
S; S.(s)+4g s;D8, (4 s;) 
Br 
und gehen außerdem mit dem Ansatz 
FI 
,: IB; tn) S,(o s;) e si 
et 


unter der Voraussetzung, daß die Konstanten B; von der Zeit r abhängen, in die Gleichu 
hinein‘). Es entsteht dann für das “te Glied 


S.los;) d B; Be 
) + ‘ v “u I 0) ) r 0+» Y | ! _— T 8;8 Q 
B;e-"si".48,(os;) 2 FT —- s? B; e-"8i"8,los;)+ e-"83i7.,(0 
u s;8, (8) 4 S;j Sig S;) a dr 


! 


s, Vgl. aneh Wilh. Müller: Einige instationäre Bewegungen der zähen Flüssigkeit unter zylindrischer 
erenzung, Annalen der Physik, 5. F., Bd. 25 (1935), S. 185 bis 204, insbes. S. 200; ferner ZAMM Bd. 15, 1959, S.372 I 


Vorträge der Hauptversammlung, Stuttgart). 
I 
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| daraus wegen 


LA 5” 
sender Ausdruck für DB; 
T 
‘) . 
B, = - a - Fir) BAR: 4 2.4. “re (37) 
ö S; Ss, (Ss) rg s;9S,(4 s;). 
\ 
| mithin die gesuchte Lösung’) 
L T SH T 
, S (os;) e="8i? ae E ‚7 (s)8S,(os;)e-tsi? fr __ z 
> \ o\E a. Ir (r) e’ si’ d 7 7 \ 0 t TE ı |Fners#a: (35). 
s;»D,(8;) +4S:ı98,(48i). PS J/,(s;) r/.tg ur 
i ) 0 i | 0 


onn die Gesehwindiekeit am Anfang über den Querschnitt konstant und gleich U gesetzt 
ir. so ergibt sieh die vollständige Lösung durch Verbindung von (35) mit (23). Man erhält 








| 
lo FL jr 2 . 
) u Hy > ” O \ ri (S;) Do (0 s;) e : a I S N [ Y 0) 
0" 3 al|r: x I Fir) e’si? 30 
$ | g (1 1 l mn Js) + J.(g s;) u ES” ’ . g ar% dr “ ) 
lg ’ 2 0 


Die mittlere Geschwindigkeit wird dagegen 


I Vrlhlasd -FAlsdle®(, Su, f | 
“ 8 y' ‚T8j? 
" u | ee sg s;) +, (s;)] \ € Ss; | / u Bo don) ca 
i=] 1 


"ür den Fall des Kreiszylinders gehen diese Ausdrücke über in 


F. 7 


J (os;)) Su 4 
. .) oa) -—- 92 \ I U n—r5j2 T L- y yrsj? 
v—=2ull— 0°’)+2 dg" L PR. - [Fin a a di . . (4), 
=) 0 
3 In du R 2 
v=u+4 \ — e-tTs?|l . + | Fin ers”) ee et A 
Ss; Ss; 2 
ij l 0 


wo ım ersten Falle über den Nullstellen der Funktion (8), im zweiten Falle über den Null- 
stellen der Funktion .J,(s;) zu summieren ist. Damit die mittlere Geschwindiekeit während 
ler ganzen Zeitdauer der Bewegung unverändert bleibt, muß die Funktion Fr) der Inteeral- 
sleichung 
L T 
 J tqs;))- +J,(s;) S f 
\ 0 / ı ) q rsj? | - —- | F(t)e’si? dr () i 5 - (43) 


r e 
——— s(J. lg S;) n_ J,(s;)) \ ES; 
0 


nügen, deren Lösung für den besonderen Fall 5=0 von Vogelpohl angegeben ist. Indem 
" «lie weitere Behandlung der einschlägigen Fragen auf eine spätere Gelegenheit verschieben, 
sen wir nur kurz auf die beiden wichtigsten Funktionengruppen hin, die einem exponentiell 
ehmendem und einem periodisch veränderlichem Druckgefälle entsprechen. 


a) Wenn F die Form F(ir)=m e-“" hat, so ergibt sich für den maßgebenden Teil der 
ung der Ausdruck 


u. == .) m \ > J,(0 s;) (e ‚ si” e r a?) | | & 
3 _ = = . ’ i ’ i a j yon 
ns J (s;)s; (af s;’) 
=} 


Vel. auch G. Vogelpohl: Über die Ermittlung der Rohreinlaufströmung aus den Stokes-Navierschen 
ıngen, ZAMM, 6 (1933), S. 25. (Vorträge der Hauptversammlung in Würzburg d. Ges. f. a. M. u. M.) 
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Wenn a einer der Nullstellen, z. B. s; gleich wird, so ist, wie ein einfacher Grenzübere: 
zeigt, das entsprechende Reihenglied dureh 
>) ‚Fe a. sk? 
In re i 
Sk; J, (S7;) 
zu ersetzen. 


b) Haben wir F(r) = wsinpr, so kommt im Falle des Kreisrohres'") 


\ ” n} Io 8;)e "3° 
5) £ N 1} > ! 


\ Y s;./,(o s;) \ ’ J,los;) 
J a n J COS PT 
—— j 8; 1J,18;) u Sl -8;)J,(8;) 





Dieser Ausdruck entsprieht der Longitudinalschwingung einer Flüssigkeit im Kreisrohr bei 
periodischer Druckänderung, wenn für #==0 Ruhe herrscht. Es ergibt sich daraus u. a., «dal; 
die Geschwindiekeitsänderung gerenüber der Druckänderung eine Phasenverschiebung zeigt. 
die sich mit wachsender Frequenz 5 vergrößert''). 


6. Strömung in einem zylindrischen Rohr bei gegebener Anfangsverteilung der Geschwindig- 
keit und zeitlich veränderlichem Druck. Schließlich wollen wir noch den allgemeinsten Ausdruck 
für die Axialbewegung in einem Hohlzylinder angeben, wenn der Druck mit der Zeit sieh 
ändert und die Anfangsverteilung der Geschwindigkeit (r), „=P (o) gegeben ist. Wenn wir 
die Gleichung in der Gestalt (35) annehmen, so läßt sich die Geschwindigkeit aus dreı Teilen 
zusammensetzen: 


| 
lo . 
Inu ’ 5 0 ; Ge 2 
1 | 0” (1 4") | r A;S,(os;)e-"s tv, : i . ‚ ® (16) 
le ii 
W 


und die Koeffizienten ‚1; des zweiten Ausdrucks (als Lösung der homogenen, vom Druck 
unabhängigen Gleichung) können dadureh bestimmt werden, daß man das oben eingeführte 
Integrationsverfahren verwendet und die Orthogonalitätsbedingungen (11) der Besselschen 
"unktionen, sowie die Formeln (14) bis (16) berücksichtigt. Man erhält dann nach ent- 
sprechender Vereinfachung folgenden Ausdruck für die Geschwindigkeit 


| i | L l 
u °o0 1” v I (s;) N, (Oo S;) P si? f 
/ | 0” ( | ) - r-2 \ e Bw = «D (0) S, (oO s;) odo 
: i Ä 2 Je(gs)— Jr (8) - we . Ba 
4 N l 4 
1 (47). 
"J (s;)S,(os;)e- "SI f_, ” Sul 
7 \ u — j (r)e’si’ dr „| 
dm J,(s) + I, (y Ss; Es; | 
! | () 





der für 5.0 übereeht in 


FL l L T 
Yy J » . y . | . ll | 
d /,(oSs;) J.(os;) i , Su 
)‘ I (| 0°) —+ ) \ - e- : | P(o).J, (os;))odo | >) \ en | F (r) e’si” dr — | (4 
| selrl) Tag I, (sh. s; 





! | u ! | u 


Dabei ist vorausgesetzt, daß die Funktion P(o) nach den Besselschen Funktionen S; (0 : 
bzw. ./,(os;) entwickelbar ist. Die Anwendungen auf besondere Fragestellungen sind 
manniefaltie, daß sie einer besonderen Untersuchung vorbehalten bleiben müssen. Ich h 
im besonderen bei einer späteren Gelegenheit zu dem Fall der Schwingungsbewegung ein 


Beitrae liefern zu können. d 


- 


10) Vgl. aueh P. Szymanski: a. a. ©. NS. 94, wo der Fall nach einer anderen Methode behandelt ist. 
11) Die Lösung ist inzwischen vom Verf. in der unter ®) zitierten Arbeit zahlenmäßig ausgewertet und graph 
veranschaulicht worden, 





ang 


AM: 


r 


indig- 
druck 
bt sich 
In wir 
leilen 


(416) 


Druck 
führte 
schen 
ı ent- 


(47). 


aph 





| 16, Heft 4 


ıgrııst 1036 


Collatz, Über das Differenzenverfahren partieller Differentialzleichungen 239 








‘ber das Differenzenverfahren bei Anfangswertproblemen 
partieller Differentialgleichungen. 


Von L. Collatz in Karlsruhe. 


rotz der Erfolge, die das Differenzenverfahren bei Randwertproblemen elliptischer Differential- 
eleichungen zu verzeichnen hat, ist es bisher nur wenig zur Behandlung von Anfangs- 
tproblemen partieller Differentialgleichungen verwendet worden. 

Der Grund hierfür hegt wohl an einer namentlich bei Auftreten höherer partieller 
leitungen in der Differentialgleichung ungünstigen Fehlerfortpflanzung, die, wie das Beispiel 
ner Schwingungsdifferentialgleichung vierter Ordnung in I zeizt, das Verfahren vollkommen 
hrauechbar machen kann, wenn es in der üblichen Weise verwendet wird, indem die Differential- 
notienten dureh Differenzenquotienten ersetzt werden. Die vorliegende Arbeit behandelt nun 

‚Ilezemein die Fehlerfortpflanzung und bringt insbesondere zwei neue Gesichtspunkte, die der 
\sıfstellune von „brauchbaren“ finiten Ausdrücken dienen sollen. Ein ausführliches numerisches 
‚eispiel in 6 zeigt, daß das Differenzenverfahren bei Verwendung der auf diese Weise gefundenen 
niten Ausdrücke an Stelle der Differenzenquotienten zufriedenstellende Resultate liefert. 


1. Beispiel, daß das Differenzenverfahren in der gewöhnlich verwendeten Form unbrauchbar 
sein kann'!\. Es liege das Problem vor, die Differentialgeleichung 


er af 


te 3 a re 


unter den vorgegebenen Randbedingungen 


of,  ofl2. m) 


V == fi, Bo ee rede 

IV, = f(2,Yy) pr ar (=) 
oO) bc TEN . 
XL, S x ei de Bi ; . 

dy fi. 0) In 4 ar i ir) (>) 


iin Gebiet OT =2, 0=y zu lösen?) Physikalisch kann man f(r,y) als die Ausbiegung 
der Stelle = eines transversal schwingenden Balkens der Länge 2 zur Zeit fies ist F propor- 
tional y) deuten. Der Balken ist an beiden Enden eingespannt, zur Zeit 79° 0 gemäß der 
vorgegebenen ungefähr sinusförmig aussehenden Funktion f(.r,0) ausgebogen und wird dann 
sich selbst überlassen. 


Wir gehen genau wie auch sonst beim gewöhnlichen Differenzenverfahren vor, legen 
ein reehteckiges Gitter zugrunde, das aus den Schnittpunkten der Geraden 


weh Er ia ER 3. len Er 


besteht, und suchen Näherungswerte F 


„> für die Funktionswerte f(h,=zl. Dazu ersetzen 
wir üblicherweise 


"fh, j ats Fa tfı.s-: 
25 durch - 
y I 
BER 
"fühl lure] F\ t 2» X F', t „+6 F..x | P, 1, X r I, 2X 
pr durch ji 


| 


' entsprieht dann die Differenzengleichung 


TR ER. tr Pus-ı ı Bra IF ı.,2+6F, AN 7 RR 0 (5) 
l? h' El re 


!) In den bisher bekannt gewordenen Fällen, in denen die Lösungen des Differenzenverfahrens nieht gegen die 
x des Differentialgleichungs-Problems konvergieren vgl. die Arbeit von R. Conrant, K. Friedriehs und 
ewyv: „Über die partiellen Differenzengleichungen der mathematischen Physik“, Math. Ann. 100 (1928), 32 bis 74, 
.. Collatz: „Das Differenzenverfahren mit höherer Approximation für lineare Differentialgleichungen“, Schriften 
lath. Sem. u. d. Inst. f. ange. Math. Univ. Berlin, Bd. 3, I bis 34), lag die Niehtkonvergenz stets daran, daß die 
\laschenweite, d. h. die Maschenweite in Fortpflanzungsriehtung, zu groß gewählt war, so daß man schneller 
ärtseine, als es dem Fortschreiten der Charakteristiken entsprach. Bei Verkleinern dieser Masche nweite trat 
" Konvergenz ein. In dem vorliegenden Beispiel jedoch liegt das Versagen des Differenzenverfahrens nieht daran, 
te Zeit-Maschenweite zu groß ist, sondern ist schon im Ersetzen der Differentialquotienten dureh Differenzenquotienten 
det. 


dia) dig?) 
dx dx“ 


°) Die Differentialgleichung (1) ist wegen 0 aueh in den beiden Endpunkten erfüllt. 
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of N 


Wir setzen nun, genau so wie man es bei der Wärmeleitungsgleichung Nr” nr - zu tun pflegt 


h? BL ’ . . : r s ; ; . i . . . i FL 


Bei dem Ansatz (6) nimmt die Differenzengleichung (5) die Form an: 


77 7 1 Y 77 7’ 
Fiat: Fi+2.% IF +1, 2 ur ser F-.% Fux-: 4 
Wählen wir als Maschenweite etwa h = 50 entspricht den Anfangsbedingungen (2) 
= er: 2 EG 


(aus Symmetriegründen können wir uns auf das Gebiet == 1 oder 1: 4 beschränken) und (3) 
entspricht: 


F,..=g(:h) (;). Pu GL 238. : : 3:3... 


Kür die erste Zeile (Y I) besagt (7) 


.) F r2,0 7 F, 71,0 i Ft F, 10 N PR) (10), 


für die weiteren Zeilen ergeben sich nach (7) die Zahlen des folgenden Schemas: 


Nach dem gewöhnlichen Differenzenverfahren berechnete Näherungswerte F, 











«X 
! | l () { 1 l 2 ( 3 l | 
/ 1 1 1 ; 
A / 1 1 ( l 1 1 am 
| | 2 | 
() () V.OZ3O0T1TDD (0) 0.030071735 V.I11I805S158 V,107568754 0,214601837 
| 0.04465260>? () 0.044652602 0.109863 744 0.176862238 0.203013 1841 
2 0.125 VUOO25TIOB () 0.000257993 0.125533042 0.134975166 0.168515 
3 0321214435 0) 0.321214435 | — 0.239579119 0,324199479 0.0483413 1? 
| 0.25 > SSSS46125 () 2, S88846123 3.462780402 3.?2?8868720 3.0976034175 
iM) 31.203006508 0 + 31.203006508 | — 41,.180005336 | —+- 44.950525752 | — 45.098583152 


Schon an diesen wenizen Zahlen erkennt man die völlige Unbrauehbarkeit des Verfahrens 
in der hier angewandten Form, da derart springende Zahlen nicht als geeignete Näherungen 
für die Ausbieeungsfunktion angesehen werden können. Es wurde mit so vielen Dezimal- 
stellen gerechnet, damit nieht etwa das Anwachsen der Abrundungsfehler als Grund für die 
Unbrauchbarkeit angesehen werden kann, sondern nur die Formel (7), nach der gerechnet wurde 


In 6 soll gezeigt werden, daß dasselbe Beispiel bei Verwendung „besserer“ Formeln als 
(7) eine brauehbare Durchreehnung zuläßt. 


2. Allgemeinere Fälle, in denen das bisher verwendete Differenzenverfahren unbrauchbar is! 
Die Formel (7) hat folgende zwei unzünstiren Eigenschaften: 


I. Die F,,, haben dem Betrage nach große Koeffizienten im Vergleich zu dem Koeffizient: 
I des Gliedes F,,,; ,, das berechnet werden soll. 

2. „Benachbarte“ F,,,- Werte werden mit verschiedenen Vorzeichen genommen. 

Diese beiden Eigenschaften wirken zusammen dahin, daß die zu berechnende Größe F',,, 
als Differenz großer Zahlen erscheint. Wie stark auf diese Weise ein vorhandener Fehle: 
anwachsen kann, zeigt sich deutlich aus der folgenden Fehlertabelle, in der zur Illustratı 
einmal angenommen sei, es bestünde kein Verfahrenfehler und die Funktionswerte F,,,, F 
seien alle bis auf einen einzigen F'„,, exakt richtig. 


s) BR. Courant, K. Friedrichs, H. Lewy a.a. O.t) S. 68. 











(10). 
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Fehler der F,.z-Werte 








f a f 7 1 { 7 2 f a 3) 
% | & ( () 0 
% 2 hi 1: 0) 
% 3 19: 10: 248 NS: 
% | D60 & 196; IT 12 6 
% #) 68939 ı 6200: 


den Zahlen der Tabelle in 1 ist nicht die Fortpflanzung der Abrundungsfehler, sondern 
der Verfahrenfehler entscheidend. Die Abrundungsfehler könnte man durch Mitschleppen 
»r zu Anfang großen Stellenzahl klein halten, die Verfahrenfehler jedoch nur durch Ver- 
«lung besserer Formeln als (7). 


Die gleichen Schwierigkeiten treten allgemein auf, wenn man Anfangswertprobleme bei 
‚tiellen Differerentialgleichungen höherer Ordnung nach dem Differenzenverfahren in der 
ise behandelt, daß man Differentialquotienten durch Differenzenquotienten ersetzt, z. B.: 


wi w\._ 
‚Pr | l ae Were 2 et I ,o 


on f l 
»@ dureh 
( y"' ( I’ 


und die so entstehende Differenzengleichung nach dem Gliede F,,o...„ auflöst. (o ist an sich 
n 
7 .) 


Differentialquotienten einzeln durch höhere annähernde finite Ausdrücke ersetzt, also z. B.: 


I» 
1) 


' gesetzt.) Auch wenn man die 


D|° 


villkürlich, wird aber gewöhnlich gleich 2 oder 


Ofı.% f ‚x+2 BPR.S3,; Be; f 
z (dureh : 
07, 121 


hleiben dieselben Schwieriekeiten bestehen. 


3. Zwei neue Gesichtspunkte bei der Aufstellung finiter Ausdrücke Zur Aufstellung 
biauechbarer finiter Ausdrücke stehen zwei Auswege zur Verfügung. 


A. Zur deutlieheren Fixierune der Vorstellung beschränken wır uns im foleenden auf 
ineare partielle Differentialgleichungen m-ter Ordnung bei 2 unabhängigen Variablen ., y: 


\ all ; 


LKif) A) ARTE 


Be Fi RO 


s seien etwa auf der a-Achse die zur Eindeutigmachung von f(x, y) nötigen Anfangsdaten 
rzegeben, und es werden Näherungswerte F,,, der Lösungsfunktion f Gr, y) für y > 0 gesucht. 


\lan pflegte bisher bei Zugrundelegung des rechteckigen Gitters (4) finite Gleichungen zu (11) 


Sn n a. 6 

PA SAG OPee 7 "oe re ee ; © : 
‚ustellen, indem man 

EEE AO a a ee 


’ 


dem Taylorschen Satz nach f und den Ableitungen von f an der Stelle «,,z, entwickelte 
die €,,,„ so bestimmte, daß (13) mit L (fi...) bis auf Glieder möglichst hoher Ordnung 
einstimmte. Die Ausdrücke (12) waren von der Form, daß ein C,,, mit dem größten z, 
,r auftrat, und (12) wurde benutzt, um F;,, aus den F,,, mit z<k zu berechnen. 
"ı wurde der Punkt :,,2, möglichst in die Mitte des „Sternes“ gelegt, der von den in 
auftretenden (von O verschiedenen) €,,, gebildet wurde: war L(f) symmetrisch, so wurde 
auch der Stern der €',, „ symmetrisch gewählt. Zur Erlangung brauchbarer finiter Ausdrücke 
“ jedoch oft wesentlich, daß der Nullpunkt der Taylorentwicklung :,.z2, von der Mitte 


"0 
S1 


ernes verlagert wird nach der Fortpflanzungsstelle hin, d.h. in Richtung wachsender y. 


16 
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Von wie entscheidendem Einfluß diese Verlagerung sein kann, zeigt 5. 


B. Man nehme außer den zur Übereinstimmung von L(fi,,x,) mit den entsprechen 
(liedern der Taylorentwieklung von (13) benötigten (,,„ noch weitere Punkte, also weit 
zur Verfügung stehende Unbekannte €,,, hinzu, die nicht einer höheren Approximation die: 
sollen, sondern einer Verkleinerung des „Index“ des finiten Ausdruckes. Dabeı verstehen 
unter dem Index ./ eines finiten Ausdruckes (;,; Fi; +2 0,.z F 


l _— 6; 


,,z. bei dem nach dem Gli: 


“<kK 
F;, , aufgelöst werden soll, den Ausdruck 
l y 
J ' \ ( / , (] 
Girl ww 
uch 


Der Index eines Ausdruckes der Gestalt (15) wird also erhalten, indem man (12) nach di 
(Gliede mit dem größten 2 auflöst: 


und «die Summe der Beträge der Koeffizienten aller anderen F,,, bildet. Wır können die 
Ergebnisse von I dahin zusammenfassen: Das Differenzenverfahren versagte bei dem Beispiel 
ın 1, weil der Index des dort verwendeten finiten Ausdruckes (7) zu groß war. In 6 win 
vezeigt werden, wie stark Zusatzglieder C,,, den Index herabdrücken können. 

4. Fehlerbetrachtung. Stimmen in der Entwicklung von (13) nach dem Taylorschen Satz 
die Glieder mit denen von L(fi,.,) bis auf Glieder r-ter Ordnung (es sei r = m) überein, so 
lassen sich die Restglieder unter der Voraussetzung, daß sämtliche (r + 1)-ten partiellen 
Ableitungen von fin dem betrachteten Gebiet existieren und dem Betrage nach durch f" 
abgeschätzt werden können, zusammenfassen zu 9:D. 7 "'f7#», Dabei ist !9|=1, D ein 
vom Gitter und vom Differentialausdruck Z, aber nieht von f abhängende abschätzbare Größe, 
hängt also ım allgemeinen von den Maschenweiten A und ! ab, und zwar ist D in! ein Polynom 
von höchstens m-tem Grade: also 


P 0 HRRpee pre I, Fri -U DIN u ERTER, 


.. 


Nach (12) gilt daher für den Fehler /,,,=F,.,-fı.z 
Base BDETM Ei, DE e dia rin uk 


Da man auf der Anfangsgeraden y=0 alle gewünschten Ableitungen von f aus den 
gegebenen Anfangsdaten berechnen kann, kann man auf allen „Anfangsgeraden* y-%, 
yl....9y=(N DI (in dem finiten Ausdruck (12) mögen höchstens N +1 aufeinander- 
folgende 2 vorkommen) Näherungen F,,, berechnen nach 


r 
177 _Vvıv f 
F., £ 7 I \ 2 ı,o . . . . . . ° . " . ’ ” o (16). 
". ( „ 
v=] 
so daß 
!’ + 1 „'! 4] 2 
En YTzr Aa BEN EN<=T 3 un 0 
S il)! 


Mit den Formeln (15) und (17) ist die Frage nach der Fehlerabschätzung im Prinzip erledir! 


Ist der Index / 1°), so kann man die Fehlerfunktion Y,,, im Falle r > m explizit so a 


schätzen, daß die Konvergenz der Lösungen des Differenzenverfahrens gegen die Lösung der .\ 
fangeswertaufgabe evident ist. Wir geben eine Schrankenfolge O< Y,<Y,<Y,<...< Y,: 
an mit 7. Y.: (man kann (12) dureh ©, , durehdividieren, also C; „== 1 erreichen) 


u Max ( L; .; Er ER Ei v)—+ Dir-m-+ 1 ger 4 2 


- a 0° of 4 x ne 
4) Das ist z.B. der Fall, wenn man der Wärmeleitungsgleichung I die finite Gleichung 
8x2 NR 7] 
F, # IF, # l 4 F # F, A 
l a) I k l 
hh? I 


mit den Maschenweiten A, I „ahh? gegenüberstellt, wo @ eine beliebige Zahl aus dem Intervall 0 <a<] 
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lieı 


Gl 
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en die 
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6b wirt 


n Satz 
in, so 
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ya 
D ein 
(röße. 
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(1Dı. 


us den 


B=%W%, 
nander- 


(16). 


(1%) 


rledıiyt 


t so dl!) 


der A 


n) 
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ser Forderung genügen wir in bekannter Weise durch’) 
r —— 4 ur m 1 »«T — 1) 
I; u. pe "Ye zu Di / 


Er Y, 1, D I’ m 1 gır4 


‚hlen wir als Y, das Maximum aller Ausdrücke (17), so folgt 


(N j)7+ ı ır4 1 u 2 rm Ifır 1) 


einer festen Bereichstelle #,y, (y=z:1=const) geht die Fehlersehranke mit ! gegen 0, 
un nur r > m ist oder wenn im Falle r m das Polynom D die Maschenweite I als Faktor 
thält. (Das ist z. B. stets der Fall, wenn in der betrachteten Differentialgleichung (11) die 
ınktion f undifferenziert nicht vorkommt, d. h. A, ®). 


5. Anwendung des Gesichtspunktes A. Im folgenden will ich für einen einfachen Fall, 
dem der Index / >1 ist, den Weg zu einer expliziten Fehlerabschätzung andeuten. Zugleich 
|| dieser Fall den Einfluß des ersten Gesichtspunktes aus 3 illustrieren. Es werde die Gleichung 


"fo „of 


I; (/ Ir, Y) ) Q r? ‘ A) y 


ge, y)f !(x,y) 
hei konstantem e betrachtet. Als Maschenweite wählen wir h und !=-eh. Wir suchen 
Iinite Ausdrücke höherer Annäherung, die in ihrer Tavlorentwieklung bis auf Glieder 4. Ordnung 
einschließlieh mit Z(f) übereinstimmen. 

Je nachdem, ob man die Taylorentwieklung auf die Mitte des Sternes (den Punkt ı,2 — 2 
oder auf den Punkt ı,2 1 oder auf die Spitze des Sternes, den Punkt ı,% bezieht, ergeben 
sich (zwangsläufig) die folgenden finiten Ausdrücke: 


FF, =16F,,x-ı Eh’ g.x- fun, +16 F,.x-5 Firz-ı | (15) 
Pr IP 4,0 BE Rt, » 1Eh’f,,x-: 
1, = - Erg aa ya Bl. an Fun st Fs- 
FRE ae Fre (FF 41.% Eee (19) 
(F+.% ? -F, 1, 4-2) ISh’t,,x 1 
85 +12’, ,.) FM. =32F .x-: Br. 2 
Sb(H, 1, X ‚+#, 1, X ı) 32H, H1,X-2 t F, 1, X%—2 (F, 2,%-2 T N, 2,%-2) (Er). 
IF in m FF, ; 3) 12h’t,,. 





Diese drei Ausdrücke sind ohne die #-Glieder übersichtlicher in folgender Tabelle, deren 
Darstellungsweise wohl klar ist, zusammengestellt: 

















Ausdruck (18) Ausdruck (19) Ausdruck (20) 
ade Ze ie (35 + 124? g,.x%) 
16 2 -(4+12Rgq,.x-,) 12 36 32 36 
16 12h gu,u.2— 161 1 4 12 4 -dH—1 8 a8 2 —ı 
+16 ) 12 32 12 
1 1 11 








of 0° 


Ausdruck (18) ist derselbe, den man erhält, wenn man 3 und er einzeln dureh 
Or” Of 


»r approximierende, nur die F,,, einer Geraden benutzenden Ausdrücke ersetzt und ist 
den in 2 angeführten Gründen unbrauehbar. Beim Übergang von (18) zu (19) wächst der 
tor des Gliedes F,, im Vergleich zu den übrigen Faktoren stark an, der Index wird 
entlich kleiner. Allerdings muß man das „Brauchbarwerden“ des finiten Ausdrucks mit 
ı Komplizierter-Werden, d. h. einer Vergrößerung der Rechenarbeit in Kauf nehmen. 


») Ähnlich wie bei R. v. Mises: Zur numerischen Integration von Differentialgleichungen, ZAMM 1030, 8.86 bis SS. 


16* 
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Wir skizzieren die Fehlerabsehätzung ım Falle 9,.,„=0 für das Rechnen mit der finit 
(il. (20), die wir in der Gestalt schreiben 


t 
\ Y \ Y 12 
n, () 7" 3 i 
Fir ax Fi-uk-x 3, i,k; 
—— — en 
f % | 


für die Lösungsfunktion f(x, y) gilt 


fir an. 35h Lf,d—Aur; 
%=] 


{ 
v Y 12 
\ \ a, 


das Restglied /,,, kann dabei abgeschätzt werden durch 


f . WM, n = 

rl >ZA= ng "14 + 105 ce +580 ce + S10 ce! +0 ce! + 1759 ce?) , 
wo M, der Maximalbetrag aller fünften partiellen Ableitungen von f im betrachteten Gebiet ist. 
Also gilt für den Fehler 


N 
y » y Y a 
P, h Fi; h liı.k 6. = a Y; uk-x TAi,k; 
I 4 l 


drückt man in gleicher Weise die Y,,,_, durch 7,, mit k—-5=z2= k— 2 aus, so folgt 


Is 
i 
pP; h > 2a, a Y nn - 8-1 r hi, k 7 Da, - M-uk-ı 
! “=] l 
So fortfahrend erhält man 
! k-4 
pP; l, a K,, „ih u pP, NER WERE a 2 ki, l; u 220, Nude 7; u k-vy ä (21 ), 
et %=] t 2m] 


Während sich nun die «a,,„'” in sehr unübersichtlicher Weise durch die Ausgangsgrößen, die 
.. . » . zer . y ’ Y ’ 
a,,,” ausdrücken, bestehen einfache Relationen für die Summen a,” = N'a,,"”’, und es 


I 
kommt gerade auf diese Summen an, da sich bei gleichem Vorzeichen der «,,,'” (bei festen 
x,» für verschiedene ı) nach (21) sofort die Abschätzung ergibt: 


1 k—4 


V, P ® Mi N ur ’) 
p .\< > la, Pia _,+al +, ai \ 
v—] 


(22) 
wobei Y, = Max /,,,| sei. 


Entsprechend der Entstehung der «a,,"” gilt 


I» 
(+1 s eo’) ıv) (DD 
(dl, = d,} N na, dl, 
Ks genügen alle a,” (d.h. fürz = 1,2,3,4, es ist a,” 0 zu setzen) der totalen Differenzen- 
' 
gleichung 2, 5) a,’ 2,00, für die man den Lösungsansatz z,= {," durchführen kann, wo 
“, die Wurzeln der Gleichung «* — a,” x! 2-0 sind. Im vorliegenden Fall hat 35." 101. 
a Iv+4:iy>3 Ä 
+ 114 0° — 562 +11=0 die Wurzeln 1;1; — 5. I” so daß man alle a,” dureh a, 
.)») 
Const.» abschätzen kann. Aus (22) folgt daher 
PU, > const k- + const k?- 2, 
wo für den Anfangsfehler Y=Max(V/,. VW, Y,) nach (16) AD) V= const- 1? erreie 
werden kann. An einer festen Bereichstelle (= kl= const) ist mithin |; | = eonst - ". 


Für einige weitere Typen von Differentialgleichungen stellen wir in der folgenden Tabe 
finite Ausdrücke höherer Approximation zusammen, die nach den in 3 aufgestellten Gesich' 
punkten gefunden wurden (der Übersichtlichkeit halber sind an Stelle der Indizes i,k uı 
i-ı,k- x die Indizes 0,0 und 1, x geschrieben: A, B,C,t sind Funktionen von X, 








ıingr 


M« 


nit 


et ist. 
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Tabelle: Einige finite Ausdrücke Y(F) höherer Annäherung. 
0” f of 
‚(f ) :—-B-—; rUCf: 
KN=Ay a Byytref=t 
2 (P) = Con Fan } rl BF, + 104 F,_,— 114, ,+56F,,—1UF,- 
A, 0 - Y j p) Y I) / 5) y 
opel FF, +36 (F, ‚+/# 1» 4 66 F,,, v S2(F,, ‚+/ I» >) 
(F', 2 t F_ 2» >) 24 F. 3 I 12( A 3 1 F 1» »)] AR 
telied: 
)—-tuolf): - (14 + 1051? 14380 hP + 810 P + 970 e - 1002 +17%0|B-P 
/) 0,0 % . | | h fiat! 1? . 


mit |9|S1. 


"ir 1 B ein einfacherer Ausdruck mit einem um I niedrigeren Annäherungsgrad®) (h=1): 


A 
Ku ( F) v ER F 0.0 + 34 3F 0,0 + 3(F, 0 + F I», - ı) 2 F %,—2 (F = + l 1» 


"8h 
4 F', n san Fa 
v0 L(f) + Restglied 4. Ordnung. 


0? f 


020 Y 





L(f)= A +Bf=t 


2, MB, #,r : auf > Fin —- (Fu. ıtF-.0) +7T2F_,-ı+21(F.->+F_.0) 
24(F_,-.+F_, +) rs. 
Er, „+. DBek 
Yu) = I:(f) + Restglied 5. Ordnung. | 


x 
— 
(u 
Du 
> 
=) 
. 





of „of 


L(f)= Fr b\ y +Cf=t, Maschenweiten I, k mit 2-1 bh’: b sei eine Konstante. 
1: ah E To 0 (F) u 30 + 12 h? SER Fo + 16 F,. 1 r F 1» 1) E. 4 r F 2, 1) 124° Ins 


=L(f)+Ph'R mit |dISI1;: 


\ f 
ol 


[? setzt sich aus Maximalbeträgen partieller Ableitungen von f zusammen (ist von h unabhängig). 


m 





Aufstellung eines „brauchbaren“ finiten Ausdrucks für die Balkenschwingungsgleichung 
(\nwendung der Gesichtspunkte A und B). Zum Schluß wollen wir noch zeigen, daß unsere 
\lethoden ausreichen, um das in 1 genannte Beispiel in Ordnung zu bringen. Zunächst brauchen 
wir einen geeigneten finiten Ausdruck 2C,,F,, für den Differentialausdruck (1). Wir 
bezeichnen die C,,„ mit a,b,e,d,e,f,g wie in der Darstellung: 


xıax 
XXX 
GXEXG 


N 
> u _ 
Ar e IN e ZN 
ar WW. 
4 xTx9g, 


ei wir ausnutzen, daß die beiden Gitter (das Gitter, das die Punkte mit a,b,e,d,e,f,g 

ält und das Gitter, dessen Punkte durch Kreuze angedeutet sind) sieh nieht beeinflussen, 

wir so bei den Maschenweiten 2-1=- A? zum Vordringen bis zu einer Stelle y 9, die 
sche Anzahl von Funktionswerten berechnen müssen wie in 1. Dann ergeben sich für die 
bei einer Taylorentwicklung an der Stelle a die Gleichungen: 


%) Der Ausdruck verwendet nur F,,., bei denen © +% gerade ist. Das Gitter aller ı,x-Punkte zerfällt dann in 


sieh nieht beeinflussende Gitter, und man kann die Rechnung für jedes Gitter getrennt durehführen; jedes Gitter 
eh liefert einen Überbliek über den Funktionsverlauf. 


Ztsehr. f. angı 
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Koeffizient von f: a+2b+ c+ 2d+ 2e+ + 29=0 
If, b+2c+ 4d+ 6be+ 4f+ Sy=O 
hi? | | | 
la 2b +8d+2e + 8g=0 
P 2 (> 
> Tun: 2b +4c+ Sd+1Se+16f+32g Ts 
hl ; ® . .) 
.) E X y: Ih T I6d-+ he 921 =—() 
h" ; .) | *).) .) *»).) 24 
oylx vr: 2b +932d- ae 920 MR 





Die Koeffizienten von fr: fr fr xx fr y„» verschwinden wegen der symmetrischen Annahmen 
über a,b,e,d,e,}.y von selbst. Es sind hier alle Glieder bis h* einschließlich aufgeschrieben 
Diesen 6 Gleichungen könnte man auch ohne die Größe g genügen, es wird sich aber zeigen 
(als Anwendung des Gesichtspunktes B), daß g einen starken Einfluß auf den Index hat. Die 
Auflösung des Gleichungssystems (23) liefert als finite Gleichung für (1) 


(2 IG)(F,, ; E Pr, X -1 F\ 2-1) r2+2G)(F,x +F,.x J+(ü GH, 9» x -2 


F_.z2-.))-2+46G)(F, 41.x- + f),-.%-.)+ Golf, 42.2, + F 
We 2 ! Is; 3 I,» X 3 {I 2, % 


{ l 2, X 


(24). 
4 „=V 
wo @=4gF° noch willkürlich ist. Den Index J als Funktion von @ aufgetragen, veranschan- 
lieht die Abb. 1. Für G=0 ıst der Index /J=7, während der kleinste Wert von J. der 








Lan mul an mer ne u 
. 


Abb. 1. 


... | 
Wert 3,5 für @: „ angenommen wird. Der finite Ausdruck wird also viel besser, wenn 


wir die Punkte mit g hinzunehmen, Das Beispiel von 1 ist im folgenden noch einmal mit 


| 
dem finiten Ausdruck (24) für G = > 
SF S(F, 67 E 1, 2-1) -Z FF, tur, u 4 2X +2P,% | 
(2 
(FF + en | 


lureheereehnet. Auf der Zeile yeıd sind die Werte B: F' 


{0 


nach (3) vorgegeben, (8) und () 


werden ohne Anderune übernommen. Bezeichnen wır die Werte von F auf den Geraden 


y=0,4,21 mit 4,B,C,D,a,ß: 
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ind die Werte €, D auf den Geraden y= + I! nach dem gewöhnlichen Differenzenausdruck (5) 


‚ h? gereehnet”); für die beiden zunächst noch unbekannten Werte a,p auf den Ge- 


n y= +21 ergeben sich nach (25) zwei lineare Gleichungen: 
Sa-SsC+24A+3(B+A)+2a- lat pP) 0 
Ss —S(C+D)+2b+5(4+5b)+2P -la+P)=d. 


it kennt man die F/,,, auf den Zeilen I, y-0,y=I, y- 21 und der Anwendung 
(25) zur Berechnung von F,,, auf den weiteren Zeilen y=>3T1, y=4l,.... steht nichts 


rim Wege. So ergab sich das Schema der Näherungswerte F°): 


Näherungen F,.z. berechnet nach (25) für das Problem (1) (2) (3). 














> > 
| L— | UV ! | ( ( - m 4 AIR | 
t=— 2=0 = r - r 2 1 summen Sy, 
4 4 2 | 

( 0) 0,03907 0.03907 0,11806 0.18757 0.271460 0,45328 
0 0,11600 0.21170 0,44571 
> 0.03870 0.035870 0.17837 0.413413 
> 0 0.10614 0.109288 0.410516 
| 0.125 0.033683 0.053363 0,15446 V.BT61S 
ß 0) V.OS669 V.I5HT1S8 0,33056 
0.02521 0,02521 0,11726 0.284093 

() 0.05943 0.105852 0.227838 

0.25 0.01481 001481 007011 0,16982 

() 0.027383 0,0508 0.100555 

10 0.000330 0.003330 0.017535 0.041?9 
| N 0.00657 0.OTOSO 0.023983 

0.375 0.0082? 0.0082 0,03636 0.080915 

Ö 0,03936 0,0TONS 0,14969 

V.OISSB 0.0188 0.086290 0,21025 

15 ( 006837 0.12425 0.26099 
0.5 0.027776 002776 0.125811 0.31175 

0 0.09122 0.165098 034842 

0.03426 — 1.053426 0,15829 0.358510 

0) 0.10594 0.19269 0.140457 

3) 0.625 0.03771 0,.03771 — 0,17491 0.412405 
0 — 0,1112? 0.20231 0,412475 

0.03779 0.039779 0,17494 0.412545 

() 0.110666 0.191419 0,40751 

0.755 0.034951 0,0341 0,.16027 0.580957 
27 0 VVHIDTS | 0.16023 0,3478 


Die ersten Maxima der Ausbiegunz ergeben sich somit zu 


Mn = 0,2023, m.# 


3 


0.100. 


.) = m | ” ’ 
m, 0.21 16, Hl, Ds Fun 1,42 


\\ie aueh sonst beim Differenzenverfahren erster Annäherung wird man an die hier er- 
haltenen Werte keine zu hohen Genauigkeitsansprüche stellen dürfen (man bedenke, daß in 
ler Spalte der letzten Tabelle durchschnittlich nur 5 Zahlen zwischen Maximalwert und 
\ullage berechnet wurden, für genauere Werte müßte man kleinere Maschenweiten benutzen, 


beonm Reehnen mit h fr wird m, > 0,212, das ist etwa 1°/, Fehler), sondern der Wert des 
2 

(ahrens wird in der Mögliehkeit gesehen, mit recht geringer Rechenarbeit einen Überblick 

" den Verlauf der Lösungsfunktion zu gewinnen, aYB 


’) Dieser Ausdruck ist auf der Anfangsgeraden y=T brauchbar wegen F,,ı Fı. 1, Der Faktor des zu 
hnenden Gliedes F,, ı verdoppelt sieh mithin (vgl. (10) ), während diese Stärkung des Faktors des zu bereehnenden 
\os bei den folgenden Zeilen fortfallen würde und das gewöhnliche Differenzenverfahren also weiterhin nicht 
brauchbar ist. 

3 Mi 
; Fs, vn 


“, Das Bilden der Zeilensummen 5, +1, für gerades „, S, =D Fs,,, für ungerades # liefert eine 


ı 0 } 1 
tige Kontrolle: Sy + Say +: =2-Sıxryı für alle 
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Tabellen für Besselsche Funktionen erster und zweiter A 
ö > | 3 
mit den Parametern » +35. +7: +, 


Von K. Karas ın Prag. 


nläßlich einer demnächst erscheinenden Untersuchung über die Eigenschwingeun: 
A ınhomogzener Saiten zeigte es sich, daß die strenzen Lösungen, die sich in Besselse] 
| | 2 | 3 
Funktionen mit den Parametern vr ta, +5; + 14:27 darstellen ließen, mit and». 
weitie erhaltenen Näherungslösungen nicht befriedigend übereinstimmten und insbeson: 
bei einieen Fällen höhere Schwingungeszahlen ergaben, als sie sich mittels des Ritzsel 
Verfahrens ermitteln hießen. Dies heß sich darauf zurückführen, daß die bisher vorhande: 
Tabellen 1. zu großen Argumentsehritt aufweisen, 2. z. T. nieht genau genug, 3. z. T. fehle 
haft sind. 


Während die Besselsechen Funktionen erster und zweiter Art mit den Parametern 


"+, in dem Tafelwerk von G. N. Watson') mit dem Argzumentsehritt 0,01 bis zum 


.) 
Areumentwert 16 bis zu 7 Stellen genau berechnet sind, hegen für die Besselschen Fink 
2 | 3 
tionen erster Art mit den Parametern vr = +2. + ‚+ nur die 4-stelliven Tafeln von 
.) { | 
A. N. Dinnik°) mit dem Argumentsechritt 0,2 bis zum Argumentwert S vor, während di: 


Besselsechen Funktionen zweiter Art überhaupt nicht berechnet sind. 


In Zahlentafel 1 wurden daher die Besselschen Funktionen erster und zweiter Art 
.) 

mit den Parametern vr = +5 für den Argzumentschritt 0,1 bis zum Argumentwert 10 auf 5 
+.) 


Dezimalstellen genau durehwer neu berechnet und die Dinnikschen Werte bloß zur Kontrolle 
benutzt. Hierbei wurden die Reihen’): 


I 


\ Y bc | / | r r 24 yr I, ? ‘ gi u \ 
m Uliy)lliiv +-i)\2 2 /v) \ 2.(22v+2) 2.4.2 r+2),(2v +4) 


u 


pP (,r) 


mit //(2)= /'(z--1), wobei /' die Gammafunktion bedeutet. benutzt. Die Besselschen 
Funktionen zweiter Art wurden in üblieher Weise *) dureh 


| 


r.(2Y:=eotle u) I, 1) 
sın (vr) 


„() 


berechnet’). Die Ergebnisse wurden bezüglich Jı» (“) durch Bildung der Differenzenstaffe! 


3 
bıs zur 3. Ordnung, bezüelieh Yı»(.r) mittels der Relation 


) 2 le) Pn ir). 5 } (a) 
[2 
kontrolliert). Dabei zeigte es sich, daß die Dinnikschen Werte F>(w)-- 115 Je (a) | 
3 +) 3 
.) 
»1S und Fe (a) Il ,|-J 2 dr) für ce 48 unrichtig sind und dureh 0,5451 (statt 0,599 


bzw. dureh 0.3529 (statt 0.3560) ersetzt werden müssen. Zur Berechnung der Reihen wurd 


1), G.N. Watson: A treatise on the Theory of Bessel Funetions, Cambridge 1922, insbesondere S. 714 bis 
°), A.N. Dinnik: Tafeln der Besselsehen Funktionen Jıı und Jı2, Archiv der Math. und Physik | 
3 | 
18. Bd. (1911), S. 337 bis 338. \. N. Dinnik: Über die Stabilität der flachen Form der Biegung (Über die K 
erscheinungen), Annalen des Donschen Polytechn. Institutes zu Nowotscherkask (1913). - A.N. Dinnik: Tafeln 
Besselschen Funktionen J,ı und Jı2, Archiv der Math. und Physik 111, 21. Bd. (1913), S. 324 bis 326. 


)), P,Schafheitlin: Die Theorie der Besselschen Funktionen, S. 11 Gl. (7), Leipzig und Berlin 1908. Es sei | 
auch auf das demnächst im Verlage von B. G. Teubner erscheinende Bueh: „Die Zylinderfunktionen und ihre 
wendungen* von Prof. R. Wevrich, Brünn, verwiesen. 

\ E. Jahnke und F. Emde: Funktionentafeln mit Formeln und Kurven, S. 93, Leipzig und Berlin 1928. 

5), Es möge bemerkt werden, daß die Besselschen Funktionen zweiter Art z.B. bei P. Schafheitlin a.a 
S. 44, Gl. (4) mit entgegesetztem Vorzeichen eingeführt werden. 

6) E. Jahnke und F. Emde: a.a. O„,NS. ». 
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Zahlentafel 1. 


Fünfstellire Besselsche Funktionen erster und zweiter Art für die 
> 














Parameter + 9, mit dem Argumentschritt 0.1 bis zum Ärgzumentwert 7 10, 

» (r) .(r (X) _e(r) r J»(r) e(.r) Yelr) . (r) 

3 3 3 B 3 3 3 
0 X Sc x 5.0 0.359712 0.192416 001605 0.530125 
0.15011 2,72981 3.235785 1.714939 5.1 0.35025 0.210928 V.ODOOS 0,2778 
0.23723 1.658082 2.077851 1,24435 2.2 034010 0.243544 0,U8475 0.25216 
0.308552 1.235374 1.650273 1.016855 5.3 0.326853 0,26477 011704 0.22452 
0.360981 0,96247 1.32487 0.08270 5.4 0.310653 0.28312 014758 0.109522 
0.412331 0,7683) 1.153161 0.93240 2.) 0.29169 29838 0.17613 0.164555 
0.47006 0.615955 0.98217 0.809817 2.6 0.27025 0,3104 0.202455 0,13282 
0.51063 0.ASTS6 0.855815 0,5712 5.7 0.2465 0,31928 0.22632 010036 
0,5453 0.370691 0.75007 0.84731 2.8 0.220856 0.532485 01,24757 V.OGTAS 
0,57440 0,27799 0,65262 0.82376 3.9 0,19344 0,32711 0,266083 0,03451 
0,59795 0.188734 0.6270 0.799019 6,0 016460 0.532616 0.28158 V.OO1T6 
0.611609 0.10626 0,ATSAO 0.171275 6.1 0.131462 0.532202 0.29412 V.OBOAS 
(1.628092 0.03070 0.309855 0.743094 6.2 0,10381 0.531481 030358 0,06189 
0,63656 0.0390] 0.532247 0.71251 6.3 0.007247 0,30462 — 0.300990 0,00219 
0.653912 0.105327 0.249795 0.607837 6.4 0.040902 0.2016] 0.531309 012111 
! 0.063673 0,16233 0.180158 064151 6.5 0.009046 0.275094 0.31316 0.148390 
1.6 0.620955 0.216300 0.1137 0.650206 6.6 - 0.02162 0.2780 0.31016 0.173850 
I. 0.651777 0.26527 0.05036 0.56019 6.7 0,05202 0.253740 0,30416 0,19713 
® 0.650159 0.30926 0.VO9TS 0,D1611 6.8 0.085145 0.214985 0.209526 0,21817 
5 0,8124 0,34828 V.UHHHS V.ATOOS 6,0 0.100966 0,19077 0.283509 0.23676 
0.235697 0.38232 0,1190 0,42240 70 0.136535 0,16903 0,26930 0.2276 
0.529305 0,41158 0.16957 0.373539 re 0.161385 0.138504 0.25257 0.266504 
0,49750 0.4348 0.215495 0.32338 if 0.184495 0.11006 0.235358 0.2765: 
0.463551 0.415464 0.257537 0,27273 1.3 0.205383 0,08139 0.21256 0,28414 
0.426593 0,468092 0.29520 0.22179 7.4 0.2240] 0.0523 0.180973 0.8886 
R 0,3872] 0.47837 0.32882 (.17092 1. 0.240158 0.025311 0.169539 0.29067 
6 0.345991 V.AS311 0,3581 0.12050 1.6 0.25377 0.00594 0.139966 0.285960 
| 0.30300 0.48323 0.38305 V.OTOSS 1.7 0.264609 0.0345 0.11292 0,28569 
„O 0,25886 0,17890 0.408354 + 0,02241 7,8 0,27285 0.06245 0,08542 0,27901 
0.215387 0.417029 0,41957 0.02457 1.9 0.27822 | — 0.089357 0,05743 0.26966 
0.116841 0.415760 0.43115 0.016973 Su 0.280799 0.111508 0.020283 VIHTTS 
0.122858 - 0.44104 0.438833 0.11275 8,1 0.28055 0.139534 000108 0.214350 
0,0776 0.412087 0.441116 0,1533 8,2 0.2775 0.161092 0.026753 0,22700 
0.033305 0.39737 0.43976 0.109126 8.3 0,27184 0.182763 0.059394 0.20845 
0,01051. — 0,37081 0.143424 0.22622 8.4 0.206352 0.20129 V.U8029 0.18807 
0,05271 0,34151 0.142478 0.25804 8.5 0.253270 0.21775 0.105553 0,16608 
0.009322 0.30980 0.411154 0.28650 5.6 0.230951 0.231856 0.120944 0.14270 
0.131795 0.27600 0.309476 0.31146 8.7 0.224111 0.248351 0.151790 0.118519 
0.16796 024048 0.37465 0.33279 8.8 0.200667 0.252603 0.172539 0.009279 
0.201659 020358 0.235149 0.35037 8.9 0.18739 0.250916 0,19106 006676 
01,23254 0,16566 0,3255 0,36416 3,0 V,16647 0.26504 0,20769 (,01035 
0.26045 0.127099 0.209712 0.3741] 91 0,14414 0.264299 0.221096 0.013855 
0,28517 (LUOSS?; 0.266952 0.538023 92 0.12062 0.262091 0.233094 | — 0,01251 
030658 0,04943 0.23408 0.385254 33 V.UHG1T 0.2589 0.243145 0,03846 
0.524583 0.01105 0,20013 038111 9.4 0.071038 0,25246 0.2501 0,00374 
1.338095 0.026585 0.156500 0.376004 35 0.045495 0.243955 0.254099 0.088813 
01,34977 006312 0,12906 0.336744 3,6 0,0969 0.23233 0.256090 011140 
0.356598 0.0826 0,09264 0.535547 7 V,UOGOO 0,21892 0.25625 0.13332 
- 0.360357 0.13171 0.056509 0.340531 3,8 0.03136 0.205348 0,2307 0,15369 
0,36060 0.165320. -+- 0,01975 0.32216 997 — 0,0616 018619 0,24742 0.172934 
10.0 G.O8015 0.16723 0.230537 0.185910 
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Zahlentafel 2. 
Vierstellize Besselsche Funktionen erster und zweiter Art für die 
Parameter + mit dem Arzumentsehritt 0,1 bis zum Arzumentwert 2 8, 
r Jı(z +8) Yılz) Dane: r Jı(#) 3:1) il) y  % 
4 i 1 4 % 4 4 4 
() () . sc IC 1.0 - V,3748 0,3595 0,1336 - 0,171 
0.1 0,5207 1.7200 1.9117  — 0.9836 1.1 0.3816 0,3364 0,0941 0.20% 
(,> 1,6155 1.4310 1.408? 0.606 % 0,3845 0,3106 V.OD4S 0.23 
0.3 0.6743 1.2722 1.1249 0.3186 1.3 0.3836 0,2823 0,0156 0,260 
0.4 0.7145 1.1560 0,0205. — (1.1458 1.4 0.3789 0,2519 - ),0227 0.2894 
0,D (1.7417 1.0596 V,TD6N 0.0107 4.5 0,3705 0.2197 V.OSOS 0.3043 
0.6 V.T7D89 0.0737 V.6181 0.0995 1.6 0,3587 0.1861 0.0955 0.321 
(1.7 0.7677 V.NIB6 0.4960 0,1921 1.7 0.3436 0.1513 0,1296 0.3546 
O.S 0.7690 0.S170 0.3864 01.2705 {.8 0,3255 0,1158 0.1617 0.3445 
0,0 0,7637 0.7425 1,2864 0.3375 1.9 0.3046 0,0799 0.1916 0.3500 
1.0 0.7522 0.6604 0,1945 0.5944 2.0 0,2810 0,0459 0,2189 0.3535 
1.1 0.1392 0.972 0,1004 0.4425 2,1 0.2551 0,0082 0.2439 0.3526 
1.2 0,7120 0.5260 0,0310 0.4822 3.2 0.2272 0.0269 0.2652 0,5482 
1.3 V.HS59 0,4556 0,0416 0,5144 33 1.1976 0,0611 0.2840 0,3405 
1.1 0,6545 0,386 0.1083 0.2394 2.4 0.1666 0,0940 0.2995 0,5206 
1.5 V.H10> V.B180 (0,1695 0.5577 DD 0.1345 0.1254 0.3118 0,3156 
1.6 1.80 0.2512 0.2251 0,5695 2,6 0.1017 0.1550 0.3209 V.2O8S 
1.7 0.5383 0.1861 0,2751 0.792 9.7 V.OH84 0.1825 0.3269 0.2792 
1. 0.409937 0,1229 0.3109 0.5753 DS — 0.0351 0,2075 0,3290 0,2574 
1.0 0.4467 V.OHP>U 0,3590 ().5697 2,9 0.0020 0,2305 0.3280 0.2333 
>0 VIOTS 0.0036 0.3927 0,290 6,0 m 0.0306 0.2506 0.3238 0.2073 
2.3 0.3475 0,0520 0.4210 0.5434 6.1 0,0623 0.2679. — 0.3166 0.1798 
A 1,2962 0,1045 0.4440 0.5234 6.2 0.0929 0.2823 0.3063 0.1504 
2,3 01.2119 0,1936 0.4615 0,4991 6.3 0.1221 0.2936 0,2931 0.1200 
2.4 0.1923 0,1992 0.4740 0.4712 6.4 0.1497 3015 0,2771 0.0901 
3.5 0.1406 0.2410 AStA 0.4308 6.2 0,1755 0,3069 0.2587 0,0590 
> VUNUT (LITSS V.4ND3N 0.4054 6.6 0.1988 0.3090 0.2382 0.0270 
A | 0,0304 0.3124 (0.481? 0.3681 6.7 0.2200 0.3079 0.2154 + 0,0052 
2.8 VUOO> VSFIS 0.4742 0.3288 6,0 (1,2387 V.BOB8 0,1909 0.039> 
2,9 1,0560 IZHON 0,4627 1.2876 6,0 0.2548 0.2967 V.1648 0.009 
3.0 0.1006 1.3875 0.4474 0.2452 7.0 0.2681 V,2868 0.1375 0.092 
3.1 0.1429 0.4037 0,4280 0.2016 7.1 0.2785 0.2742 0,1093 0,119 
3.2 0.1824 0.4154 0.4051 0.1574 7.2 0.2861 0.2591 0.0803 0,14) 
3.9 1,2180 0.4227 V.BZTSI 0.1131 1.9 (1.2907 0.2415 V.OSOS 0.10% 
3.1 0.2521 0.4256 1.3498 0,0691 71.4 0.2925 0.2219 0,0213 0.191 
3.5 2810 0.4243 0.3182 0.0256 v5) 0.2913 0.2003 - V.OOSO 0.211 
36 BON OLÄISS 1.2842 0.0169 1.6 0.2873 0.1770 0.0370 0.22% 
3.4 0.3306 1.4004 V.2484 0.0581 1.1 0,2805 0.1522 0.0653 0,244 
3.8 0,3493 USE T ENG, 0.2110 VODTS TS 0.2711 0,1263 0,0925 0.257 
3 0,5640 0.3790 0,1727 0.1359 7.9 0,2591 0.0094 0,1185 0.267 
S.0 0.2449 0.0719 0.1432 0.274 
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Zahlentafel 
Vierstellire Besselsche Funktionen erster und zweiter Art für die 
Parameter + , mit dem Ärgumentschritt 0,1 bis zum ÄArgumentwert 8. 

Js (x) J s (Tr Ya (r) 4 8 (r) E Js (x) J s\E) Ys,(r) a i%) 

4 4 4 4 N 4 N 4 
() x X x 4,0 0.1934 0.1114 0,3509 0.3849 
0,1149 2,5824 3.7670 2,7449 4,1 0.2243 0.0719 0.3260 0,5801 
0,1924 1.4802 2984 1.7615 4,2 0.2524 | — 0.0326 0.2085 0,3895 
> 0.2589 1.0423 1.7329 1.4084 4.3 0.2774 0.0061 VHS 0,386? 
| 0.3180 TTV 1.4168 1.2267 4,4 0,2992 0.0438 0.2373 0.3793 
5 0,5711 0,5899 1.2053 1.1147 1,5 0,3176 0.0803 0.2040 0,3689 
6 0,4187 0.4442 1.0469 1.0363 4.6 0.3325 0.1152 0.1696 0.3550 
(1.7 0,4612 0.3277. — 0.9246 0.9799 4,7 0.3459 0,1482 0.1343 0,3381 
ILS 0.4956 0.2194 0,8089 0.0245 1.8 0.3D1S 1791 0.0985 1.3184 
0.0 0.5311 0.1275 0,7111 0,8784 4,9 0.3561 0,2076 0.0625 0,2960 
1.0 1,987 0.0447 0.6219 0,8348 2,0 0.3569 0,2336 0,0265 0.2711 
14 0,5812 0,0302 0.5385 0.7917 1 0,3542 0.2567 V,LOOSS 0.2442 
1.2 0.DOSS 0.0985 0.4595 0.7483 5.2 0,3482 01.2770 0.0435 0.2154 
1.3 0.6115 0.1607 — (), 3842 0.7041 3.3 V.B3S8 0,2941 0771 0.1850 
1.4 0,6104 - 0,2172 0,3122 0.HDSS nd 0.3264 0.308? 0.1095 0.1534 
1.D 0,6224 0,2684 0.2428 V.H118 D.D 0,3110 0,3189 0.1400 0.1209 
1.6 0,6208 - 0,3143 0.1763 0,5636 2.6 0,2929 0,3264 0.1687 VUSTS 
1.7 0.6145 0.3950 0.1125 0,5140 5.7 0,2723 0.3306 0.1952 0,0545 
1.8 V,HOBS 0,3006 00914 0.4633 D.S 0.2493 0.3314 0.2104 0,0212 
1.0 U.DSSS 0.4?11 +- 0,0067 0,4116 5,9 0.2243 0.3290 0,2410 VOL1S 
2.0 0.5608 - (1,4467 0.0619 0.3591 6.0 0,1976 0,3234 0.2598 0.0440 
2,1 0.5469 0,4672 0.1135 0.3062 6,1 0.1603 0,3148 0.2759 0.0754 
2,2 0.5204 0.4829 0.1625 0.2531 6. 0.1399 0.3033 0.2890 0.1055 
2.3 0.4906 0.4936 0.2075 0.2002 6.3 0.1096 0.2889 0.2990 0.1339 
2.4 VADTS 0.4996 0,2487 0,1478 6.4 0.OTST 0.27>0 0.3060 0.1607 
2.5 0.4223 0.5008 0,2859 0.0964 6.5 0.0475 0.2526 0.3097 0.1854 
6 0,3844 0,4977 0.3195 + 0,0459 6.6 0.0163 0.2312 0.3107 0.2081 
2,4 0,3445 0.4901 0.3486 0.0029 6.1 0.0145 V.>OTS 0.3084 0.2283 
2.8 0.3029 VATSA 0.3737 0.0500 6.8 0.0447 0.1827 0,3031 0,2459 
A 0.2600 0.4628 0,3945 0.0951 6.0 0.0741 0.1561 0.2949 0.2609 
v 0,2162 — 0,4434 0.4109 0.1376 7.0 0,1024 0.1284 0.2840 0,2732 
| 0,1719 0,4204 0.4226 0.1773 7.1 0.1293 V.OIO8 0.2704 0.2827 
4 0.1274 | — 0.3943 0.4302 | — 0.2141 712 0.1545 0.0707 0.2545 (1,2892 
. 0,0830 0,3656 0,4340 0.2482 173 V.1778 0.0413 0.2362 0.2027 
4 0.0391 0.3343 0.4337 0.2790 7.4 0.1990 V.O118 0,2157 0,2932 
I V.OOZ8 0.3007 0.4291 0,3061 1. 0.2180 0.0173 0,1935 0,2910 
6 0,0459 0,2651 0.4204 0.3294 7.6 0.2347 0.0459 0.1608 0.2860 
Ad VOS5D . — 0.2281 0.4081 0.3490 27 NV.DASS 0.0737 0.1446 0,2782 
B 0.1238 0,1899 0.3924 0.3650 7.8 0.2602 0.1004 0.1182 0,2676 
) 0,1598 0,1509 0.3732 0.3769 79 0,2689 0.1257 - 0,0911 0.2546 
8.0 0,2749 0,0633 0,2392 
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die Potenztafeln in K. Havyashi’) und die fünfstellige Tafel für & von J. Kuno°) benu 
Bezüglich weiterer eventueller Unterteilune des Tafelschrittes sei auf die wertvollen Arbei 
von F. Emde*) verwiesen. 


In der Zahlentafel 2 wurden die Besselschen Funktionen erster und zweiter Art 
Pr «+ i 3 
den Parametern v = + m Zahlentafel 3 mit den Parametern v = + 4 T. aus der oh: 
angegebenen Reihe, z. T. aus den Dinnikschen Werten durch Interpolation *’) für « 
Arzumentsehritt 0,1 bis zum Arzumentwert S auf 4 Dezimalen berechnet. Die Besselsehon 
Funktionen erster Art wurden dann wieder durch Bildung der Differenzenstaffel bis zur 


> 


3. Ordnune, die zweiter Art dureh die Relationen 


.) 


s |y2 | 5 ] 
Y oil) Jılae) + Yıla)ı - und Y_s{2)= | Jıle) — Yılz) 
4 4 4 


4 | 1 ) a 


2 f 5 2 
kontrolliert. Hierbei zeigte es sich, daß neben anderen kleinern Abweichungen um eine 
Kinheit der letzten Stelle sämtliche Dinniksche Werte für 2=5,4 und insbesondere für 
“—=4,S, letztere ziemlich beträchtlich, von den wahren Werten abweichen. Eine weitere 
Kontrolle, wenigstens bis zur 3. Dezimale, kann bezüglich Jı1(8) aus den Reihen D (.r) 





bzw. (G (ie) der Dissertation von L. Prandtl'') errechnet werden '"?). 


Bei Berechnung der Zahlentafeln wurde der Verfasser durch Frau J. Karas weitgeheni 


unterstützt. DIL 


'), K. Havashi: Fünfstelliee Funktionentafeln, Berlin 1930, insbesondere Tafel XXXVII, S. 166 bis 169. Di 
Tafel sollte bis einsehließlieh x# unzeekürzt wiedergegeben werden, da sonst der relative Fehler zu groß ist, was « 
Brauehbarkeit der Tafel herabsetzt. Bei x* sollte außerdem für x — 4,1! richtig 79849,25229 (statt 78949,25229) stehe 
s, .J. Kuno: Teeh. Rep. of the Kyushu Imp. Univ. I, Fukuoka 1926, abgedruckt bei K. Hayashia.a. O.,>.] 
bis 110 als Tafel NANVII. FKüra 1,3 sollte es dort riehtig 1,19114 (statt 1,19108), für x=5,7 richtig 3,19092 (statt 3,1909 


für \ 6,1 riehtig 3,3381 (statt 3,32861) heißen. 
», F. Emde: Unterteilung des Tafelschritts („Untertafelung“ „Brigegs-Moutonsche Aufgabe"), ZAM! 


Bd. 14 (1934), S. 333 bis 339. F. Emde: Fortlaufende Reehnungen auf der Rechenmaschine. In dieser Arbeit find 
sieh unter F wertvolle Hinweise zur Bereehnung der Reihen und unter G solehe zur Bildung der Tafeldifferenzen un! 
riinstirer Anusnutzune der Reehenmasechine. Klektroteehn. Institut, Stuttgart. Diese Arbeit wurde dem Verfasser freu! 
lieherweise von Herrn Prof. P. Funk, Prag, zur Verfügung gestellt. Sie ist derzeit noch nieht im Druck erschien: 


10), B. Sporer: Niedere Analysis, S. 74 bis 76, Leipzig 1903. 


11) 1, Prandtl: Kipp Erscheinungen, Inaug.-Diss. München (1899) S. 25. Man vgl. auch E. Jahnke undE. Emd 


Funktionentafeln a. a. O., S. DL, 92, 104 und Tafel I, S. 106. 
1?) Man vel. auch das „Verzeichnis bereehneter Funktionentafeln*, herausgegeben vom Institut für angewanı 


Mathematik nnd Mechanik, Berlin 1928, insbesondere 8. 13, Nr. 46 und 48. 
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en BUCHBESPRECHUNGEN 
Be Ine. PAUL WERNERS, Energie-Über- eine verhältnismäßig leicht verständliche Einfüh 


‚une und -Umwandlung mit Wech- 
trom. 204 8. Leipzig und Berlin 1935, Ver 
(6. Teubner. Preis geb. ISM. 


A rt ı 


ist ein uraltes Bestreben der Wissenschaft, die 
‚erscheinungen und Gesetze möglichst alle, 
mindesten aber die Gesetze gleicher oder ähn- 

(Gruppen auf eine gemeinsame Grundformel 
kzuführen. Solange die Entwicklung und Er- 
ung noch rasch fortschreitet, wird man sich 
»iner mehr geschichtlich gegebenen Darstellung 
iven. Sobald aber einmal ein gewisser Ab 
h zefunden ist, werden Wissenschaftler und 
niker von einer höheren Warte aus rück- 
nd versuchen, das Grundproblem eines 
eren zusammenhängenden Gebietes der Wissen 
t und Technik aus der Fülle der Einzelerschei- 
en herauszuarbeiten. 


er ol N} 
ür ( 
| sc] I 


m eine 
ere für 
weitere 
n Dir) 


J 


‚ hat es der Verfasser unternommen, die Ge- 
» der Erzeurung, Umwandlung und Übertragung 
triseher Energie mittels Wechselstrom auf eine 
1 gemeinsame Schaltungsart und Grundformel 
bringen, die Betriebseigenschaften daraus ab 
‚leiten und «damit eine einheitliche Theorie der 
tunzeen. Transformatoren und Maschinen zu 
schaffen. Fernleitung und Maschinen erhalten auf 
aYR irse Weise im wesentlichen die aus Wechselstrom 
erständen gebildete T-Schaltung als Ersatz- 

schaltbild zugeordnet. 


gehend 


Die Lösung dieser Aufgabe stellt hohe Anforde- 
neen an das mathematische Rüstzeug und wird 
eicht dureh die symbolische Betrachtungsweise, 
. I. durch das Rechnen mit komplexen Größen, 
ie von Fall zu Fall dureh das Vektor-Diagramm 
ınterstützt werden. Die Darstellung ist klar und 
ibersichtlich. Das Verständnis des schwierigen 


Stoffes wird erleichtert durch die am Anfang ge- 
schene Darstellung der einschlägigen mathemati 


en Hilfsmittel. 
as Buch wendet sich wohl in erster Linie an 
en wissenschaftlich geschulten Pionier der Elektro- 
teehnik. Es ist aber auch geeignet. den jungen 
Stiulierenden in die wissenschaftlichen Aufgaben 
elektroteehnischen Fachgebietes einzuführen 
nl ihm das nötige wissenschaftliche Rüstzeug zu 
vermitteln. Als Beispiele sind Probleme der Stark- 
omteehnik gewählt. Doch werden auch der 
Schwaechstrom- und Hochspannungstechniker an 
en Buch nieht vorübergehen können, da be 
sonders die ersten Kapitel schon tief in die Theorie 
Leitungen hineinführen. Auch der Mathematiker 
ıl manche Anrerung daraus entnehmen können, 
wie seine Vorlesungen über das Rechnen mit kom 
biexen Größen durch praktische Beispiele für den 
jungen Studierenden nieht nur anrerend. sondern 
überzeugend für die Notwendigkeit einer 
sieen Beschäftigung mit diesem Zweig der Ma- 


= 1} s 
il > 


natik gestaltet werden können. 
9. Dir Druck und Ausstattung verdienen hervorgehoben 
‚„ was («d verden. 
1) st hr I . \ r r r * u 
init "reibere/Sa. W. Krug VDE. 574 
0.8.1 


tt 3,1908 a6 ER RR 
phil. LUDWIG HANERT, Marine-Studienrat in 

‚„ ZAMN "orwik, Geschützund Schuß. Einführung in 

eit find 0 Geschützmeehanik und Ballistik. 2. verb. Aufl. 

zen un! 370 S. m. 161 Textabb. Berlin 1935, Verlag 

er freu s Springer. Preis geb. 27 M. 

‚schien . ’ r 
» Inhalt der vorlierenden Auflage stimmt im 


en ganzen mit der in dieser Zeitschrift 9 (1920) 


.Emd \ 233 von P’opoff besprochenen, ersten 


| Auflare 
ein. hat aber in 


Kinzelheiten eine starke Über- 


gewanı tung erfahren. Dabei ist natürlich der Cha- 
r des Buches der gleiche geblieben. Es ist 





rune in das Gebiet der Geschützmeehanik und der 
Ballistik, die an mathematischen und physikalischen 
Kenntnissen etwa «das voraussetzt, was eine gzrute 
höhere Schule vermittelt. Heute ist ein derartiges 
Buch unentbehrlich: vor allem in der Hand des jun 
ven Oftiziers wird es gute Dienste leisten. Darüber 
hinaus möchte man es aber auch in den Händen 
unserer heranwachsenden Ingenieure sehen, die so 
verhältnismäßig bequem mit den wissenschaftlichen 
Grundlagen der Wehrteehnik bekannt werden könn 
ten. Auswahl «(des Stoffes und Darstellung sind 
sehr zeschiekt, nur wäre vor allem bei der Dar 
bietung «der mathematischen und physikalischen 
(Grundlagen eine etwas schärfere Fassung erwünscht. 
Doch das sind Kleiniekeiten, die der Verfasser bei 
einer Neuauflage, die hoffentlich im Interesse einer 
möglichst weiten  Verbreitune  wehrwissensehaft 
licher Kenntnisse recht bald nötie ist. sicher leicht 
beseitiren wird. 

Freiberr (Sa.). Fr. 
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Prof. Dr. HEINRICH REICHEL, Graz, Die 
wiehtiesten mathematischen Metho 
denbeiderbBbearbeitunevonVersuchs 
ererebnissen und Beobachtungen. 88 N. 
m. 30 Abb. Berlin-Wien 1935, Verlae Urban & 
Schwarzenberg. 


Die vorliegende Einführung in die statistischen 
Methoden, ein Sonderdruck aus Abierhaldens Hanıd 
buch der biologischen Arbeitsmethoden, Abt. V, 
Teil 10, ist von einem Biologen für Biologen ze 
schrieben. Er beabsiehtiet damit nicht nur die rich 
tigen Verfahrungsweisen an der Hand von Formeln 
uud Rechenanweisunzen zu lehren, sondern anen 


„tie Ableitung dieser Rereln aus verständlichen 
Voraussetzunzen klar oder «doch glaubhaft zu 
machen“, Behandelt werden: die normale Häufig 


keitsverteilung, einiges aus der Fehlerreehnung, «lie 
Annäherung von Messungen, die von einer Ver 
änderlichen abhängen. dureh die nach der Methorle 
der kleinsten Quadrate beste Gerade und (der Korre- 
lationskoeftizient. An einer größeren Anzahl von 
Beispielen wird die Anwendung der Formeln und 
kechenmethoden erläutert. Nieht überall steht leider 
der Verf. so über der Sache, dab er sie einwandfrei 
darstellen könnte. Neben schwerfällizen Ableitun 
gen finden sich auch hier und da Unriehtirkeiten: 
so, wenn Verf, behauptet, daß durch drei Punkte 
eine Kurve „zweiten Grades (Kreis, Ellipse. Hyp®r 
bel, Parabel oder ähnliche Kurve (!))* bestimmt sei. 
oder wenn er sagt: „Die Neirunge tra der Geraden 
im log-arithmetisch geteilten Feld bedeutet nichts 
anderes, als die Krümmung der Kurve im mumer 
schen Feld” usw. Auch die Ableitung «es Fehlers 
des Korrelationskoeffizienten ist nieht einwandfrei. 
Sie führt auch auf einen anderen Wert als auf den, 
den man zewöhnlieh unter Annahme einer zwei 
dimensionalen Gaußschen Verteilung ableitet. 

Der Verfasser dürfte also mit dieser Darstellung 
nieht durchaus seinen Zweck erreichen, doch wirt 
der Praktiker daraus an Hand der gegebenen Bei 
spiele die riehtire Verwendung der wiehtiesten 
Methoden zur Bearbeitung von Versuchserzebnissen 
erlernen können. 

Freiberz (Sa.). Fr. 


E55 
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0, SCHREIER und E. SPERNER, Kinführung 
in die analytische Geometrie und 
Aleebra. I. Bd. 238 S. (1931). Geb. 9,60 M.; 


Il. Bd. 308 S. (1955). Geb. 9,60 M. (Hamburger Ein- 
Verlag B. G. Teubner. 


zelschriften Bd. 10 und 19.) 
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In klarer lesbarer Form lieet hier die denkbar 
solideste Grundlage moderner analytischer Geo- 
metrie vor, die, der heutieen Entwieklune dieser 
Disziplin entsprechend. enz mit der Algebra ver- 
knüpft ist. Freilich werden nieht, wie in manchen 
beliebten Lehrbüchern vorwierend einschmeichelnde 
Beispiele eboten. hinter denen sich dann bisweilen 
unzweekmäßiee Beeriffe und fadenscheinire Beweise 


oder nur schwer strene beweisbare anschauliche 
Theoreme verbereen. Bewährte und vernünftiee 
Berriffe werden an die Spitze gestellt. und mit 


Ihnen wird die Theorie systematisch aufge- 
baut. Beispiele sind vielfach in den Übungsaufgaben 
untereebracht. Bei aller Strenge ist die Theorie im 
Hinbliek auf (die mannigfachen Anwendungen ge 
arbeitet. Wer das Bueh Lehrbuch dureh- 
arbeitet. gewinnt nieht nur. was die Kapitelüber 
schriften ankündigen: affiner Raum. lineare 
(Gleiehuneen: euklidischer Raum. Determinanten 
theorie: Körpertheorie. Fundamentalsatz der Alge- 


» 
«ıS 


bra: Gruppentheorie: lineare  Transformationen, 
Matrizen: projektive Geometrie sondern besitzt 
darüber hinaus z. B. die orthogonalen Transfor 
mationen und die metrische Geometrie, die pro- 


jektive, affine und metrische Klassifizierung 


der Flächen zweiten Grades, die Elementar- 
teilertheorie in besonders einfacher unıl 
durehsiehtieer Darstellung). und er hat nur einen 
leiehten Schritt zu tun, um die nichteuklidische 
Geometrie und die geometrischen Grundlagen der 


analvtischen Mechanik und der Relativitätstheorie 
u. a. verstehen zu lernen. Beeriffsbilduneen und 
Sätze der analytischen Geometrie gewinnen bekannt- 
h dann ihren höchsten Grad von Einfachheit 

in einer der Verallee- 


I: . 
IC erst 


umd Übersiechtliehkeit. wenn sie 
meinerun® auf z Dimensionen fähileen Weise formu 


liert werden. Dieser Höchsterad wird überall im 
Buche erreicht. wie überhaupt die Beweisführung 
mit Vorteil so eestaltet ist. daß sie auf die allee- 


meinsten Fälle anwendbar bleibt. in denen die Sätze 
noch welten.  Nireends aber wird einem Veralloe 
meineruneswahnsinn eefrönt. der die alleemeinsten 


Buchbesprechungen 


Fälle auch da ausspinnt. wo sie überhaupt kein oder 


doeh für den Leserkreis (des Buches kein Interesse 
mehr haben: der dareebotene Stoff ist im Geerenteil 
weise beerenzt. 


Dresden. Threlfall. 


607 


Handbuch der Vermessungskunde von weil. Dr. 


phil.h.e. W.Jordan,. Prof. a. d. Techn. Hochsech. 
zu Hannover. fortzeesetzt von weil. Dr. ©. Rein- 
hertz. Prof. a. d. Techn. Hochsch. zu Hannover. 
l. Bl. Auszleiecehunges-Rechnung nach 
ler Methode der kleinsten Quadrate. 
S., erweit. Aufl. bearb. von Dr.-Ing. e. h. O0. Eeg- 
eert, Prof. a. d. Techn. Hochseh. zu Berlin. VII 

632124 S. m. zahlreichen Abb. Stutteart 1935. 


J. B. Metzlersche Verlagsbuchh. Preis eb, 30 M. 


Neben dem bekannten Buch von F. R. Helmert 
über die Auseleiehungesreehnune kann man das von 
W, Jordan beeründete, von Ü. Reinhertz fort 
und von O0, Erzert weitereeführte Werk 
beachtenswerteste bezeichnen: eine Neuauf 
beilentet deshalb immer in den in Frare kom 
menden Kreisen ein Ereienis. Das vorliegende Buch 
ist zunächst für den Geodäten bestimmt: bietet 


vesetzte 
als (las 
laoe 


. 


e 


aber aueh anderen Fachleuten. die sieh mit der Aus- 
eleiehun®e von Beobachtungen beschäftieen. viel 
Wertvolles. Dies eilt insbesondere von dem ersten 
und dem letzten Kapitel. 

Das erste, mit „alleemeine Theorie der kleinsten 
Fehlerguadratsumme"  überschriebene Kapitel ent 
hält alles, was theoretisch und praktisch notwendig 
ist zur Auseleiehunge nach vermittelnden Beobach 


tuneen. nach 
vermittelnden 


und nach 
Bedineungsglei- 


beiineten Beobachtunren 
Beobachtunzen mit 


Ztsehr. f. an: 
Math. und M 


chungen. Mit Rücksicht auf die Bestimmung 
Buches wird die „Theorie der Fehlerwahrschein! 
keit“ in dem für die Ausgleichungsreehnung eı 
(terlicehen Umfang erst im letzten Kapitel behan 
Das zweite Kapitel enthält die Ausgleichung 
„T[riangulierungsnetze“ oder die Ausgleichung 
Punktbestimmungen nach bedineten Beobachtun: 
Das zerenüber der letzten Auflage erweiterte dı 


Kapitel umfaßt die „Punktbestimmung du 
Koordinatenausgleichung“ oder die Ausgleich 


von Punktbestimmungen nach vermittelnden Bi 
achtungen: die Erweiterungen beziehen sich ins 


sondere auf die zur Beurteilung der Genauie) 
einer Punktbestimmung wichtige mittlere Fel 
ellipse. Das vierte Kapitel enthält neben Angal 


über die Bestimmung der Genauigkeit von Trian: 
lierungen insbesondere geschichtliche Mitteilung 
über ausgeführte Triangulierungen. 

In einem Anhang sind einige bei der Durehfü 
rung von bestimmten Aufgaben der Auseleichung: 
rechnung verwendbare Tafeln beigefügt. Ein alpha 
betisches Sachverzeiehnis erleichtert die Benutzun 


les Buches. 

Dresden. P.Werkmeister. 507 

Dr.-Ine. M. EWEIS, Reibunes- und Un 
diehtiekeitsverluste an Kolbenrin 
een. (Forschungsheft 371. Beilage zu „Forschung 
auf dem Gebiete des Inzenieurwesens“, Ause, B. 
Bd. 6. März/April 1935.) 23 8. m. 37 Abb. uni 
3 Zahlentafeln. Berlin 1935. VDI-Verlag m. b. H 
Preis 5 M. 

Der Verfasser hat es sich zur Aufgabe gesetzt, 
über diesen für alle Kolbenmaschinen so wesen! 


lichen Maschinenteil dureh Reehnunge und Versuch 
Klarheit zu schaffen. 

Die Forderune zeringer Undiehtiekeitsverlust: 
brinet für die Kolbenringe zumeist zwangsläufig 
las Auftreten größerer Reibungsverluste mit sich. 
Seit Jahren kann man daher im Kolbenmaschinen 
bau ein Suchen nach dem Bestwert beider Verlust: 
dureh Wahl verschiedener Formen und Kolbenring 
zahlen beobachten. 

Ausgehend von einer zusammenfassenden matlr 
matischen Behandlung der Reibung und der Druck 
verteilune unter den Ringen sowie des Druck ve: 
lustes hinter den Ringen in Abhängiekeit von deı 
Drehzahl werden Beziehungen abeeleitet für veı 
schiedene Strömungsarten der Gase dureh die ein 


zelnen Ringstufen. Aus diesen Gleichungen ergibt 
sich dann ein mathematischer Ansatz für die Be 


rechnung der Undiehtigrkeitsverluste bei pulsierender 
Strömung. 

In einer besonderen Reibungsmeßmaschine, bei 
der der Kolben mit den an den Ringen eingebauten 
Meßzeräten zerenüber einem beweeten Zylinder 
ruhend aufgehängt ist, werden die Rechnung®! 


Adnreh Messungen bei wechselnden äußeren Beil 
eunzeen der Ringezahlen,. der Drehzanlen und (« 


Drucke nachgeprüft. Neben einer Druckmessu 
hinter den einzelnen Rinzen mit Hilfe des Farnboı 
Indikators zelanget zur Bestimmung der Reibun: 
kräfte ein piezo-elektrischer Quarz-Indikator 
Verwendung. 

Als Ergebnis der Messungen wurden Werte | 
die Reibungeskräfte und für den Druckverlauf au 
bei Einstellung verschiedener Schmierzustände 
wonnen. Wenn auch diese Ergebnisse wegen © 
eerinzen Höchstdrehzahl der Versuchsanlage stre' 
nur für das Gebiet der halbflüssigen Reibung zelt: 
sind doch als Anhalt für die Vorausbere« 
nung und Planung bei Neuentwürfen von Kolbi 
maschinen bereits eine sehr wertvolle Unterlage. 


Holfelder VDI. 


SV sıe 


Berlin. 0, 6 











f.an 16, Heft 4 

ind M gust 167 EZ _ 

UNe erner sind bei der Schriftleitung folgende 

hein! er eingegangen (ausführliche Besprechung 

ee t vorbehalten): 

ae HARRY SCHMIDT, Prof. a. d. Staatl. Hoch- 

no \e f. angew. Technik Köthen, nicht-planm. a. 0. 

htun, ‚a. d. Univ. Leipzig. Kinführungein die 

te dr (tor- und Tensorreehnung unter 

An üceksiehtigung ihre u physikali- 

leich en Beide utung. \ 125 S. mit 20 Abh. 

B n je 1955. Verlagsbuchh. Dr. Max „Jänecke. 

h E s broseh. 5.80 M. 

aulel Grimsehls Lehrbuch der Physik. Neubearbeitet von 
Fel ‚Dr. R. Tomaschek., Direktor (des Physika- 

\ngal en Instituts der Techn. Hochschule Dresden. 


Trian: l 


i ter Band. Zweiter Teil: Materie und Äther. 
eilungen 


\utl. VI1/430 S. mit 314 Abb. Leipzig und Berlin 
Verlaxe B. G. Teubner. Preis geb. 14 M. 


r, W, LIETZMANN, Oberstudiendirektor in 
{inzen. Dr. J. JAROSCH, Direktor (ler Bunides- 


schule in Wien, Geometrie und zgeome- 


urehfiü 
Ichun: \ 
n alpha 
nutzun 


sches Zeichnen für die Ill. und IV. 
ur lasse der Realevymnasien,. Real 
er. hulen und Frauenoberschulen. 94. 
\utl. Leipzie und Wien 1935. Verla@ Franz 
ıtieke. Preis ’ M. 
d Un ‘ieh. Rat Prof. Dr. PHILIPP LENARD, Heide! 
enrin Grobe Naturforscher in Lebens- 
rschung schreibunzen. 324 S. mit 67 Bildn. Mün- 


\use, B. 
bb. un 


Nachriehten 





Preis oeh, 


chen 1929, J. F. Lehmanns Verlag. | 


10.80 M. 


Dr. phil. L. C. DUE, Die Brückenverbin 
dunestheorieundihre Anwendung zur 
Nlasseneinteilune und Klassenzu 
sammensetzung quadratischer Irra 
tionalzahlen und binärer quadrati 
scher Formen. 39 S. Kopenhagen 1936, Verlag 


levin & Munksgaard. Preis 4.50 Kr. 
M. VON LAUE und R. VON MISES, Ste- 
reoskopbilder von Kristallgittern. 


Teil IH. 56 S. mit 24 Tafeln. Berlin 1936, Verlaz 
Julius Springer. Preis in Mappe 18 RM. 


Prof, Dr. X. THUM VDI und Dr.-Ing. K. OESER, 
Gummifederuneen für ortfeste Ma 
schinen. Mitteilungen «der Materialprüfungs 
anstalt an der Technischen Hochsehule Darmstadt, 
herauseereb. von Prof. Dr. A. Thum VDI) V/7T2 S. 
mit 52 Abb. und 7 Zahlentaf. Berlin 1955. VDI 
Verlax G.m.b.H. Preis 6.50 M. 


P’rof. Dr. A. THUM V\DI und Dr.-Ing. F. DEBUS 
\DI. Vorspannung und Dauerhaltbar 
keit von Schraubenverbindungen. (Mit 
teilunzen (der Materialprüfungsanstalt an der Tech 
nischen Hochschule Darmstadt, herauszegeeb. von 
Prof. Dr. A. Thum VDI) VIl1/72 S. mit 90 Abb. 
und 8 Zahlentafeln. Berlin 1936. VDI-Verlae G.m. 
b.H. Preis broseh. 6.50 M. 


NACHRICHTEN 


oPSse! X. 
wesent! 
Versuch 


74. Hauptversammlung und 80-Jahrfeier des Vereins 
Deutscher Ingenieure. 

In Verbindung mit der Hundertjahrfeier der Teeh- 
nischen Hochschule Darmstadt beringe der Verein 
Deutscher Ingenieure auf seiner 74. Hauptversamm- 
iz die Feier seines S0Ojährigen Bestehens. Über 
00 Ingenieure und frühere Schüler der Teech- 


verluste 
Yesläufie 
nit sich. 
(schinen 


Verluste schen Hochschule Darmstadt nahmen an der Feier 
benring teil. die am 26. Mai mit der Hauptversammlung (des 


1} 


\DI begann. am 27. mit den Fachvorträgen aus 
en Einzeleebieten der Technik fortgesetzt wurde, 


| mathe et s 
nl «lie ihren Höhepunkt in einer gemeinsamen 


* Druck 


ruckver "eierstunde fand, die der Feier des 100jährizen Be- 
von dei stehens der Technischen Hochschule und (des 80- 
für veı ührigen Bestehens des Vereins gewidmet war, 


die ein "in Sommerfest im Hofe des Schlosses ließ die 


n ereibt | iuehleutsche Gastlichkeit gebührend zur Geltung 
die Be i sommen. Mit Ausflügen und Besichtieunesfahrten 
ierenier >), Mai fand die Tagung ihren Abschluß. 


90 Jahre Akademischer Verein Hütte. 


reihe \m 16. Mai feierte (ler Akademische Verein 
ulinder 'iutte sein 90jähriges Bestehen. Er wurde im ‚Jahre 
nn sb von Studierenden des damaligen ‚Königlichen 
 Madin 5 verbe-Instituts, aus der die Fechnische Hoch- 
ee schule Berlin hervorgegangen ist. reeründet. Aus 
can sem Kreise stammte die Anregung zur Grün- 

des Vereins Deutscher Inzenieure, (ler auf 


arnhoı 


a zehnten Stiftungsfefst der „Hütte“ Ins leben 
Kar en wurde, Der Verein hat sich in den 00 
en seines Bestehens manniefachen Aufeaben 
| lem Interessenkreis der Ingenieure gewidmet: 
erte | bekanntesten ist die im Jahre 1856 erfolgte 
uf an “ffung des Taschenbuches, das heute in 26. Auf- 
nde 2 vorliegt. Dieses dem Ingenieur unentbehr 
gen Werk ist ein schönes Zeuenis für die un- 
e strei (lliche Energie, mit der der Verein jederzeit 
> ie \lleemeinheit zu dienen bestrebt ist. 
sherei 
Kolbi Preisausschreiben für den Max-Eyth-Preis. 


Jare. 


'äblich des 100. Geburtstages von Max Eytl 
"DI. 6 lie 


Max-Eyth-Gesellschaft in Verbindung mit 





dem Verein Deutscher Ingenieure ein Preisaus 
schreiben erlassen. Der Max-Eyth-Preis in Höhe 
von 1000 RM. wird alljährlich am 6. Mai. dem Ge 
burtstage Max Evths, verliehen, erstmalige am 6. Mai 
137. 


Die Bewerber haben jeweils bis zum 31. Dezember 


an die Geschäftsstelle des Vereins Deutscher In- 
venieure, Berlin NW 7. Hermann-Göring-Str. 27T. 
eine Arbeit einzureichen. «deren Umfang bis zu 


» Scehreibmaschinenseiten (Format DIN A4. 30 Zei 
len) betraxzen soll. Der Aufsatz soll in allzemein 
verständlicher, fesselnder Form und in Prosa einen 
Stoff aus dem Gebiet der Technik behandeln. wobei 
Arbeiten. die sieh mit dem Wirken der Teehnik in 


der Landwirtschaft befassen. besonders erwünscht 
sind. 
Das Preisgerieht setzt sich zusammen aus je 


einem Vertreter des Vereins Deutscher Ingenieure, 
ler Max-Eyth-Gesellschaft. der Techniseh-Literari 
schen Gesellschaft und der Reichssehrifttumskammer. 

Näheres siehe in der Zeitschrift des Vereins 
Deutscher Ingenieure Bd. 80 (1936) Nr. 20. S. 612. 


Lilienthal-Gesellschaft für Luftfahrtforschung. 


\uf Veranlassung des Reichsministers für Luft 
fahrt. Generaloberst Göring, ist die Lilienthal-Ge 
sellschaft für Luftfahrtforsehung zeeründet worden. 
Die neue Gesellschaft wird sieh den Forschungs 
aufeaben widmen. «die bisher der Wissenschaft 
lichen Gesellschaft für Luftfahrt (WGL) und der 
Vereinieunge für Luftfahrtforsehunz oblaren. Die 
Wissenschaftliche Gesellschaft für Luftfahrt. die in 
diesem Jahr auf eine 25jährige erfolgreiche Arbeit 
‚„urückblieken konnte, wird sich auflösen. Ihr bis 
herizer Präsident. Prof. Dr. Schütte. und Herr 
Oberstleutnant a. D. Wagenführ sind zu Ehrenmit- 
eliedern der neuen Gesellschaft ernannt wortien. 
/u Präsidenten der Lilienthal-Gesellsehaft sind Prof. 
Carl Boseh-Heidelberz und Prof. L. Prandtl-Göttin 
een. und als geeschäftsführender Präsident Mini 
sterialrat Baeumker ernannt worden. 


Ztschr. f. ang 








256 Nachrichten Math. und M 
Der 12. Deutsche Physiker- und Mathematikertag Gesellschaft für angewandte Mathematik 
findet vom 13. bis 19. September d. J. in Bad Salz und Mechanik. 

brunn statt. Die Hauptversammlung der GAÄMM findet 


Donnerstag, 24. 9, bis Sonnabend, 26. 9., im 
schluß an die Tagrunz der Gesellschaft Deuts 
Naturforscher und Ärzte in Dresden statt. 

eungsort: Technische Hochschule, Bismarekplatz 


Veranstaltungen 
des Mathematischen Reichsverbandes 1936. 


I.Im Anschluß an (die Tarunz der Mathematiker Zur Tagung sind bis jetzt die folgenden Vortr: 
und Physiker in Bad Salzbrunn am Freitar. anzemeldet worden: 
len 18. 9. l.R. Beyer, Zwiekau: Schmiegungskegelsehn 
\. 9 Uhr: Geschäftssitzune. und ihre Anwendung in der praktischen Analı 
B. 0% Uhr: Allgemeine Sitzung. Thema: Ausbil und Kinematik. = a | 
dun®e «der zukünftieen Lehrer der Mathematik >,R. Gran Olsson. Frondheim: [Thema vorbehalt 
| ı (den höheren Schnlen. ınter besonderer Berück > Henschke, Berlin: Über Strahlgebläse m ho 
siehtirunz der Anrewandten Mathematik. Be Wirkungsgrad, insbesondere zur "Tunnellüftu 
riehter die Herren Dreetz und Zabel. ann Göttingen: Die tragende Fläche \ 
Lee, xreis[orm. 
G8ı ne abends: ‚ „usammenkunft und freie Aus s.K. Ludwie, Hannover: Über Träger mit Kin 
sprache der Beiräte des MR. eewieht und gleichem Widerstand zeeen B 
II.Vor «der Naturforscherversammlung in Dres une, 
den am Sonnabend, den 19. 9. zusammen mit 6. K. Marguerre, Berlin-Adlershof: Über die Tr: 
der Dresdener Ortsgruppe des Fördervereins fähiekeit eines Plattenstreifens, «der über 
und der \bteilune \V für den mathematischen A\usbeuluneserenze auf Druck belastet ist. 
und naturwissenschaftlichen Unterricht. 1.v. Ohnesorge,. Berlin: Untersuchungen ii 
D. Nachmittaes 5 Uhr: Herr Salkowski: „Was kann Tropfenbildung, >Strahlzerfall und Zerstäubung 
(ler eeometrische Interrieht für (lie A\usbillune mit Z/eitllehneraufnahmen). 
ler Raumansehanune tun?" s. Th. Pöschl. Karlsruhe: Die Blasiusschen Fo: 
E. Abends 8 Uhr: Zusammenkunft und freie Aus mein und der Mittelpunkt eines Profils. 


R %.L. Prandtl, Göttingen: Bemerkungen zur Theoı 
der trarenden Fläche. 
I0.E. Schapitz, Berlin-Adlershof: Über einige Pro 


sprache ler Beiräte des M. 
\nallen Veranstaltuneen können Gäste teilnehmen. 


Der Vorsitzende «des Arbeitsausschusses des M. R. bleme der versteiften Schale im Flugzeugrbanu. 
Prof. Dr. Hamel. 11.H. Schliehting, Friedrichshafen: Tragtlügeltheori 
bei Überschalleeschwindigekeit. 
r y -) SN.ol 1 @; } . EN ’" vo ar (vn 1; 
Die 94. Versammlung 12.H. Schmidt. Köthen: Strenge Lösungen 


Prandtischen Theorie der trarenden Linie. 
15. E. Seidl. Berlin: Betrachtunz bleibender Forn 
änderungen auf Grund der Arbeitshypothesen 


der Gesellschaft Deutscher Naturforscher und Ärzte 


timdet vom 20, bis 23. September in Dresden statt. 


Vorläutiges Programm der Tagung: „Formungsprinzip" und „Individualprinzip" 
Sonnabend. 19. September: ») Uhr: Abenıdvortra®z. 14.1. Sörensen,. Dresden: Krzwungzene und freie 
Sonnlaur U) September: 130 Uhr: Sitzune des Schwingungen elastischer Systeme, 

wissenschaftlichen Ausschusses. 12.15 Uhr: Ent- 15.W. Tollmien,. Göttingen: Zum Übergang von 

' \ ie se \ a ER a s gs x B ö I » sel ‘ 1 7 ırgehhe u »fir on 

hüllun® einer Carus-Gedenktafel am Carushaus. Unter: hall in Ub rschall Strömung N. 

16.30 Uhr: Eröffnunessitzune. 19 Uhr: Abend 16. E. Trefftz, Dresden: Über die Berechnung den 

vortrae. 20.15 Uhr: Beerüßunesabend /Zirkulationsverteilung bei der geraden, tragen 

Re ERS ad (len Linie. 
Montag. 21. September: 9 Uhr: Erste allgemeine - 


Pam) 


17. Umpfenbach, Berlin: Kalorimetrisches Verfahren 


Sitzung, > Uhr: Sitzunz der meidizinischen vs 
2 1 x N spin zur Wirkungseradbestimmung® bei Wassertn 


Hauptgruppe. 20 Uhr: Abendvortrag. 





binen. 

Dienstag, 22. September: 9 Uhr: Sitzung (der natur- IS.6G. Vogelpohl,. Berlin: Fortschritte in der >Strö 
wissenschaftlichen Hauptgruppe. 15 Uhr: Kom- mungsuntersuchung mit dem Ultramikroskop. 
binierte Sitzung der Abteilungen. 20 Uhr: Vor 19.6. Vogelpohl, Berlin: Hydrodynamische Lagen 
stellune in ler Staatsoper und Aben Ivortra®e. theorie und halbtlüssige Reibung. 

Mittwoch, 23. September: 8.15 Uhr: Geschäftssitzung 20.0. Weber, Dresden: Streifen mit Eigenlast. 
der Gesellschaft. 9 Uhr: Kombinierte Sitzung der 21.G. Weinblum, Berlin: Theorie der aktiv: 
\bteilungen. 15 Uhr: Zweite allzemeine Sitzung, Schlingertanks. 

Schlußansprache. 20 Uhr: Vorstellung im Staat 22.F. Weinig, Berlin-Adlershof: Betrachtung 

liehen Sehauspielhaus und Abendvortrag. räumlichen Strömune um den Tragtlürel. 

Der Preis der Teilnehmerkarte beträgt: für Alt- Ktwaige weitere Vortragsanmellungen wert: 
mitglieder der Gesellschaft 10 M.. für neubeitretende umgehend an den-Unterzeichneten erbeten. 

Mitglieder 15 M., für Niehtmitglieder 20 M., für An- lin besonderer Versammlungsbeitrag wird v' 

eehörige 10 M. Teilnehmer an den Tarunren (der der GAMM nieht erhoben, doch wird auch «dı 

Kachvereinizungzen, die an der Tagung der Gesell- Teilnehmern der GAMM-Tarıung, «die am Nat 

schaft nieht teilnehmen. können eine Sondertarunges forscher- und Ärztetag nicht teilnehmen, empfohl 

Teilnehmerkarte zum Preise von 3 M. lösen. (lie die eine Sondertagungs-Teilnehmerkarte der Ges 

Vergünstirungen gewährt, die die Teilnehmer am sehaft Deutscher Naturforscher und Ärzte 

Naturforscher- und AÄrztetar wenießen (z. B. er- Preise von 3 M. zu lösen, wofür, von anderen \: 

näßiete Eisenbahnfahrt zwischen Do. 17. 9. und rünstieungen (verbilliete Eisenbahnfahrt usw. 

Shyd. 26. 9.. gesehen, die Unterkunftsvermittlune übernomm: 
Die Geschäftsstelle der Gesellschaft Deutscher wird. Diese Sondertagungs-Teilnehmerkarte ist | 

Naturforscher und Ärzte, Dresden-X\. 16. Elisenstr. 4. ler Geschäftsstelle der Gesellschaft Deutscher N 

versendet auf Wunsch kostenlos Einladungen und fwrforscher und Arzte, Dresden-A. 16, Elisenstr. 

’rorsramm der Tazune, aus dem Näheres zu er Prof. Dr. R. Zauniek) anzufordern. 

sehen Ist. Prof. Dr. E. Trefftz. Dresden-A. 20, Cäeilienstr 
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